"reinit"
Fname="../RES/PDSI_Res/PDSI_Stn_No_Rain_Center.nc"
IDX.1="PDSI"; IDX.2="SPI"
Nidx=2; 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";

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"
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.7; x01=10.5;  y00=7.5; y01=1 
 Nst=131; Ntime=419
 
   "sdfopen "Fname
   "set display white"; "c"; "set mpdraw off"; "set grid on"
   "set xlab off"; "set ylab off"
    x1=x00; x2=x01; y1=y01; y2=y00
   "set parea "x1" "x2" "y1" "y2
   "set grads off"; 
   "set xlab on"; "set xaxis 1980 2014 "; "set ylab on"; "set yaxis 1 131"
   "set xyrev on"; "set yflip on" 
   ist=132
*   While (ist<Nst)
*     ist=ist-1; "set y "ist 
     T1_clor()
*     "set xyrev on"
     "d ori_pdsi"
*   Endwhile 
*   draw_tit(x1+0.2, y2-0.2, R.ireg)
*   draw_tit(x00+dx*icol-dx/2, y000+0.15, Var.icol)
    "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 sub-regions (R1-R7) and entire country (VN)"
*    "draw string "x00+(x01-x00)/2" "y00+0.1" for different scenarios"
*    "gxprint Pic_Apr_30/"IDX.ifil"_Changes_in_"Fnam.isev".pdf pdf"
    "set parea off";
pull dd
"close 1"
*******************************************************

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 T1_clor()
levels = " -6 -5 -4 -3 -2 -1  0 1  2 3 4 5 6 "
*"run My_Color.gs -gxout grfill -kind blue->lightblue->aqua->white-(0)->white->green->yellow->red -levs" levels
"run My_Color.gs -gxout grfill -kind darkred->red->white-(0)->white->aqua->green -levs" levels
"set grads off"
return
*******************************************************