"reinit"
   "set display color white"; "c"
   Varn.1="EOF1"; Varn.2="EOF2"; Varn.3="EOF3"; Varn.4="EOF4"; Varn.5="EOF5"; Varn.6="EOF6"

   Dat_out="/home/tanpv/Progs/Dr_SEA/PDF/"
*  Trend_spei_12.nc
   Dr_var="spei_12"

   Dat_src="APH"
   Dat_nm="observation"

*   Dat_src="Ens"
*   Dat_nm="model ensemble"

   y00=7.8; y01=0.8
   x00=0.5; x01=10.5

   "sdfopen /home/tanpv/Progs/Dr_SEA/MCA/Trend_" Dat_src "_" Dr_var ".nc"

          x1=x00; x2=x01; y1=y01; y2=y00
          "set parea "x1" "x2" "y1" "y2
          "set mpdraw on"; "set mpdset mres"; "set frame on"; "set xlab on";"set ylab on"
          "set grads off"
          "set ylint 5"; "set xlint 5"
          Cor_Clor()
          "d sen*120"
          "hatch pval 0 0.05 -angle 45,135 -int 0.1 -color 1 -thickness 0.01 -type 1 -density 0.01"
          "xcbar "x00+1" "x01-1.2" 0.3 0.5 -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on -edge triangle "
   draw_tit((x00+x01)/2, y00+0.4, "Trend (per decade) of "Dat_nm" drought ("Dr_var") during the baseline period")
   draw_tit((x00+x01)/2, y00+0.2, "(Hatching: p-value<0.05)")
   "gxprint " Dat_out "Trend_drou_" Dat_src "_" Dr_var ".pdf pdf"
   "close 1"
*******************************************************
function draw_tit(x, y, title)
  "set string 1 c 5.5 0"; "set strsiz 0.11 0.11"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return
*******************************************************
function Line_pl(mark,colo,thick,styl)
"set gxout line"
*"set gxout stnmark"
"set cmark "mark
"set ccolor "colo
"set cthick "thick
"set cstyle "styl
return
********************************************
function Bar_pl(clor,gap)
"set gxout bar"
"set ccolor "clor
"set bargap "gap
"set barbase 0"
"set baropts outline"
"set cthick 2"
"set baropts filled"
*"set ccolor 1"
return
********************************************
   function Cor_Clor()
   levels = "-1 -0.8 -0.6 -0.4 -0.2 -0.1 -0.05  0  0.05 0.1 0.2 0.4 0.6 0.8 1"
   "run color.gs -gxout grfill -kind red->orange->white-(0)->white->cyan->blue -levs " levels
   "set grads off"
   return
********************************************
function Truc_0(a0)
"set gxout line"
"set ccolor 1"
"set cstyle 1"
"set cthick 2"
"set cmark 0"
"d "a0
return
********************************************