"reinit"
DIR="/home/tanpv/SEACLID_New/RES_Post_Grd_2019/"
Fname="_T_P_Changes_Reg.nc"
IDX.1="PDSI"; IDX.2="SPI"; IDX.3="SPEI"
Nidx=3; Nreg=8; Nsce=4; Nmod=7
SCE.1="BL";SCE.2="MF45";SCE.3="MF85";SCE.4="FF45"; SCE.5="FF85"; SCE.6="OBS"
MODL.1="CNRM";MODL.2="CSIR";MODL.3="ECEA";MODL.4="GFDL";MODL.5="HadG";MODL.6="MPI";MODL.7="Ens";
**********************************************************
** Color <and> style of marks for each model
cmdl.1 = 2 ; smdl.1 = 1
cmdl.2 = 3 ; smdl.2 = 2
cmdl.3 = 4 ; smdl.3 = 6
cmdl.4 = 5 ; smdl.4 = 6
cmdl.5 = 6 ; smdl.5 = 7
cmdl.6 = 7 ; smdl.6 = 1
cmdl.7 = 2 ; smdl.7 = 5

**********************************************************
Var.1="MF45"; Var.2="MF85";Var.3="FF45"; Var.4="FF85";

V.1=dt_mf45; V.2=dt_mf85; V.3=dt_ff45; V.4=dt_ff85;
V.5=dp_mf45; V.6=dp_mf85; V.7=dp_ff45; V.8=dp_ff85;
*V.9=dgeo3_mf45; V.10=dgeo3_mf85; V.11=dgeo3_ff45; V.12=dgeo3_ff85;

*	float T_BL_OBS(Subreg, Model, Mon) ;
*	float P_BL_OBS(Subreg, Model, Mon) ;

*Seve.1="Moderate_Drought"; Seve.2="Severe_Drought"; Seve.3="Extreme_Drought";
Seve.1="Regional mean of monthly Temperature (C)"; Seve.2="Regional mean of monthly Precipitation (%)";
Fnam.1="Regional_mean_of_monthly_Temperature"; Fnam.2="Regional_mean_of_monthly_Precipitation"

R.1="R1";R.2="R2";R.3="R3";R.4="R4";R.5="R5";R.6="R6";R.7="R7"; R.8="VN"
Mon.1="Jan";Mon.2="Feb";Mon.3="Mar";Mon.4="Apr";Mon.5="May";Mon.6="Jun";
Mon.7="Jul";Mon.8="Aug";Mon.9="Sep";Mon.10="Oct";Mon.11="Nov";Mon.12="Dec"
X_labs="|J|F|M|A|M|J|J|A|S|O|N|D|"
*xt Y_labs="CNRM|CSIR|ECEA|GFDL|HadG| MPI| Ens"

*x00=0.5; x01=8.0; dx=1.5; y00=10.5; dy=1.25
 x00=0.9; x01=8.0;  y00=10.3; y000=10 
 Ncol=Nsce; Nrow=Nreg; Nsev=2
*xt dx=(x01-x00)/Ncol; dy=1.1

ifil=0;Nidx=1
While (ifil<Nidx)
   ifil=ifil+1
   "sdfopen "DIR IDX.ifil Fname
   "set display white"; "c"; "set mpdraw off"
   ivar=0; isev=0
   While (isev<Nsev)
      isev=isev+1; icol=0; 
*xt      "set ylab off"
      "set grid off"
      While (icol<Ncol)
         icol=icol+1; ireg=0; irow=0; ivar=ivar+1
*xt         x1=x00+(icol-1)*dx; x2=x1+dx; x11=x2; if (icol>1); x1=x1+0.01*(icol-1); x2=x1+dx; endif
*xt         "set xlab off"
*xt         if (icol=1); "set ylab on"; "set ylabs "Y_labs; endif; if (icol>1); "set ylab off"; endif
         ylb = "off" ; if (icol=1) ; ylb = "on" ; endif
         While (irow<Nrow)
            irow=irow+1;ireg=ireg+1
*xt            if (irow=1); y2=y000-(irow-1)*dy; y1=y2-dy; y11=y1; endif
*            if (irow>1); y2=y000-(irow-1)*dy-0.05; y1=y2-dy; y11=y1; endif
*xt            if (irow>1); y2=y11-0.04; y1=y2-dy; y11=y1; endif
*xt            "set parea "x1" "x2" "y1" "y2
            "set t "ireg; "set grads off"; "set gxout grfill"
*xt            if (irow=Nrow); "set xlab on"; "set xlabs "X_labs; endif
            "set xlabs "X_labs
            xlb = "off" ; if (irow=Nrow); xlb = "on" ; endif
*xt            if (isev=1); T1_clor();endif; if (isev=2); P1_clor();endif; 
            "AIO_plot -ir "irow" -ic "icol" -nr "Nrow" -nc "Ncol" -xlb "xlb" -ylb "ylb" -omr 0.3 0.4 0.7 0.1 -space 0.05"
            im = 1 ; ccols = ""
            if (isev = 1)
              'set vrange 0 6' ; 'set ylint 2'
            else
              'set vrange -100 100' ; 'set ylint 50'
            endif
            while (im <= 7)
               'set x 0 13'
               'set y 'im
               'set cstyle 0'
               'set ccolor 'cmdl.im ; ccols = ccols" "cmdl.im
               'set cmark 'smdl.im
               "d "V.ivar
*xt            draw_tit(x1+0.2, y2-0.2, R.ireg)
              'nx_title -main 'R.ireg' -pos top -siz 0.12 -space -0.1'
              addline(V.ivar)
               im = im + 1
            endwhile
         Endwhile
*xt      draw_tit(x00+dx*icol-dx/2, y000+0.15, Var.icol)
      Endwhile
*xt      "cbarn"
    "set string 1 c 5.5"
    "set strsiz 0.12 0.12"
    "draw string "x00+(x01-x00)/2" "y00+0.5" Changes in "Seve.isev
    "draw string "x00+(x01-x00)/2" "y00+0.3" over subregions (R1-R7) and entire country (VN)"
    "draw string "x00+(x01-x00)/2" "y00+0.1" for different scenarios (left to right: MF45, MF85, FF45, FF85)"
    "gxprint Pics_2019/"IDX.ifil"_Changes_in_"Fnam.isev".pdf"
    "set parea off";
pull dd
    "c" 
   Endwhile
"close 1"
Endwhile
*******************************************************

function draw_tit(x, y, title)
  "set string 1 c 5.5"; "set strsiz 0.12 0.12"
  "draw string "x" "y" "title
return
*******************************************************
function T_clor()
levels = " 0 0.2 0.5 1 1.5 2 2.5 3 3.5 5"
"run My_Color.gs -gxout grfill -kind white-(0)->green->yellow->red -levs" levels
"set grads off"
return
*******************************************************
function P_clor()
levels = " -100 -70 -50 -30 -20 -10 -5 0 5 10 20 30 50 70 100 "
"run My_Color.gs -gxout grfill -kind red->yellow-(0)->aqua->blue -levs" levels
"set grads off"
return
*******************************************************
function Bar_pl()
"set gxout bar"
"set ccolor 15"
"set bargap 15"
"set barbase 0"
"set baropts outline"
"set cthick 2"
"set baropts filled"
*"set ccolor 1"
return
********************************************

*******************************************************
function T1_clor()
levels = " -6 -5 -4 -3 -2 -1 -0.5 0 0.5 1 2 3 4 5 6 "
*"run My_Color.gs -gxout grfill -kind blue->lightblue->aqua->white-(0)->white->green->yellow->red -levs" levels
levels = " 0 0.5 1 1.5 2 2.5 3 3.5 4 5 6 "
"run My_Color.gs -gxout grfill -kind white-(0)->green->yellow->red -levs" levels
"set grads off"
return
*******************************************************
function P1_clor()
levels = " -200 -100 -50 -30 -20 -10 -5 0 5 10 20 30 50 100 200 "
*"run My_Color.gs -gxout grfill -kind red->yellow->green->white-(0)->white->aqua->lightblue->blue -levs" levels
"run My_Color.gs -gxout grfill -kind red->yellow-(0)->aqua->green -levs" levels
"set grads off"
return

function addline(var)
  'a = 0*'var
  'set cstyle 1'
  'set ccolor 1'
  'set cmark 0'
  'set cthick 3'
  'd a'
return