"reinit"
   "set display color white"; "c"

   Mon.1="1961-1965"; Mon.2="1966-1970"; Mon.3="1971-1975"; Mon.4="1976-1980"; Mon.5="1981-1985"; Mon.6="1986-1990"
   Mon.7="1991-1995"; Mon.8="1996-2000"; Mon.9="2001-2005"; Mon.10="2006-2010"; Mon.11="2011-2015"; Mon.12="2016-2020"

   L.1="OBS"; L.2="LT=1 mon"; L.3="LT=2 mon"; L.4="LT=3 mon"; L.5="LT=4 mon"; L.6="LT=5 mon";

   Var.0="Multivariable Flash Drought Index"
   Var.1="Num of FDs per year"; Var.2="Mean duration (day)"; Var.3="Mean severity (%)"; Var.4="Onset Decline Speed (%/pentad)"
   Var.5="Rate of Depletion (%)"; Var.6="Intensity (%/day)"
   Varn.1="Pr"; Varn.2="T2m"; Varn.3="Tx"; Varn.4="Tm"
   Dat_out="/home/cordex_sea/Minh/Out/"

   Fn.1="NumFD";  Mod.1="pr_bs_"
   Fn.2="Duration"; Mod.2="t2m_bs_"
   Fn.3="Severity";  Mod.3="tx_bs_"
   Fn.4="Speed";  Mod.4="tm_bs_"
   Fn.5="RateDepletion"
   Fn.6="Intensity"

   Nvar=6
   Ncol=3; Nrow=2 
   y00=10.5; y01=0.5
   dy=(y00-y01)/Nrow+0.25
   x00=0.5; x01=8.3
   dx=(x01-x00)/Ncol

   "sdfopen /home/cordex_sea/Minh/Out/MFDI_and_Trend_1yrs_Daily_New.nc"
   "set display color white"; "c";
   "set gxout grfill"

    x1=x00; x2=x01; y1=y01; y2=y00
       "set parea "x1" "x2" "y1" "y2
       "set mpdraw off"; "set frame on"; "set xlab on";"set ylab on"
*       "define a=ave(mpi, t=1, t=12)"
       "set t 1 "
       "set lon 101 116"
       cor(0); "d sen"
       "xcbar "x2-1.0" "x2-0.75" "y1+3" "y2-1.2" -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on "
*       draw_ver_tit(x1+1, (y1+y2)/2-0.5," Mean "Var.0" for the period 1961-2020")
*       endif
*       if (icol=1); draw_ver_tit(x00-0.1, (y1+y2)/2, L.irow); endif
*       if (irow=1); draw_tit((x1+x2)/2-0.05, y00+0.05, Mon.icol); endif
       Sub_Reg2

       HS_TS
       "Title_lonlat -lon " 112.1 " -lat " 15.500 " -siz " 0.15 " -ang "0 " -tit "Hoang_Sa""  
       "Title_lonlat -lon " 113.5 " -lat "  8.250 " -siz " 0.15 " -ang "0 " -tit "Truong_Sa""         
       "Title_lonlat -lon " 109.0 " -lat " 17.250 " -siz " 0.15 " -ang "310 " -tit "Bien_Dong""

*       "hatch pv 0 0.05 -angle 135 -int 0.05 -color 1 -thickness 0.01 -type 1 -density 0.01"
*       "hatch prob 0 0.10 -angle 135 -int 0.05 -color 1 -thickness 0.01 -type 1 -density 0.01"

   Main_tit((x00+x01)/2, y00+0.2, "Trend of "Var.0" for the period 1961-2020")
*   Main_tit((x00+x01)/2, y00+0.15, "(Hatching: p-value<0.5)")
   "gxprint " Dat_out "MFDI_Trend.pdf pdf"
   pull dd
   "c"
*******************************************************
function Main_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 draw_tit(x, y, title)
  "set string 1 c 5.5 0"; "set strsiz 0.09 0.10"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return
*******************************************************
function draw_ver_tit(x, y, title)
  "set string 1 c 5.5 90"; "set strsiz 0.09 0.10"
  "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(icas)
      if (icas=0); levels = " -0.4 -0.3 -0.2 -0.1 -0.01 0 0.01 0.1 0.2 0.3 0.4 "; endif
      "run color.gs -gxout grfill -kind blue->lightblue->white-(0)->white->orange->red -levs " levels
      "set grads off"
   return
********************************************
   function Colors(var)
      if (var=0); levels = " 90 92 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 "; endif
      if (var=1); levels = " 0.5 1 1.5 2 2.5 3 3.5 4 "; endif
      if (var=2); levels = " 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 "; endif
      if (var=3); levels = " 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 "; endif
      if (var=4); levels = " 5 10 15 20 25 30 35 "; endif
      if (var=5); levels = " 4 6 8 10 12 14 16 18 20 "; endif
      if (var=6); levels = " 10 12 14 16 18 20 22 24 26 28 30 "; endif
      "run color.gs -gxout grfill -kind lightgreen->yellow->orange->red -levs " levels
      "set grads off"
   return
********************************************