"reinit"

  mark.1=0; mark.2=0; mark.3=0; mark.4=0; mark.5=0; mark.6=9; mark.7=3
  colo.1=4; colo.2=1; colo.3=8; colo.4=4; colo.5=1; colo.6=1; colo.7=2
  thic.1=5; thic.2=5; thic.3=5; thic.4=8; thic.5=5; thic.6=5; thic.7=7
  styl.1=1; styl.2=1; styl.3=1; styl.4=1; styl.5=1; styl.6=1; styl.7=1;

 Var.1="prec"; Var.2="spi12";Var.3="spe12"; 
 Varn.1="PDSI"; Varn.2="SPI_12"; Varn.3="SPEI_12"
 Mod.1="PC1"; Mod.2="PC2"; Mod.3="PC3"; Mod.4="PC4"; Mod.5="PC5"

 Out_Dir="/home/tanpv/Progs/Dr_VN_2020/VN_PP/PDF1/"
 In_Dir="/home/tanpv/Progs/Dr_VN_2020/VN_PP/NC/"
 X_Axis="1980 2018"
 idx=1

*  File.1="Areal_mean_48.nc"
 File.1="/home/tanpv/Progs/Dr_VN_2020/VN_PP/Results/spi_spei_Reg_1980-2018_48timescales.nc"
 EOVar.1="spi"

 Dat_src="observation"

* Spatial PC patterns of Drought over SEA
 x00=0.7; x01=8.0;  
 Ncol=3; 
* dx=(x01-x00)/Ncol-0.2

 y00=10.6; y01=0.5
 Nmode=3

   "reset"; "set display white"; "c"; "set grads off"; "set grid off"
* EOF Timeseries of Drought over SEA:
   mark.1=0; mark.2=0; mark.3=0;
   colo.1=2; colo.2=4; colo.3=1;
   thic.1=5; thic.2=5; thic.3=5; 
   styl.1=1; styl.2=1; styl.3=1;
   
   Nvar=2; Nrow=8;
   x1=x00+0.0; x2=x01-0.0
*   "set x 1 "Ntime
   dy=(y00-y01)/(Nrow+0)-0.3
   
   Nreg=8; Nlen=48; Ntime=468
   R.1="R1";R.2="R2";R.3="R3";R.4="R4";R.5="R5";R.6="R6";R.7="R7"; R.8="VN"
   Nrow=Nreg
   "sdfopen " File.1
   Len=Nlen-36
  irow=0
  While (irow<Nrow)
    irow=irow+1; ireg=irow
    "set x "ireg
    if (irow=1); y2=y00-(irow-1)*dy; y1=y2-dy; y11=y1; endif
    if (irow>1); y2=y11-0.3; y1=y2-dy; y11=y1; endif
    "set ylab on"; "set ylint 3"; "set yaxis 1 "Len; "set grads off"
    "set parea "x1" "x2" "y1" "y2

*    if (irow=Nrow); "set xlab on"; "set xlab %.0f"; "set xlint 5"; "set xaxis "X_Axis; endif
    "set xlab on"; "set xlab %.0f"; "set xlint 2"; "set xaxis "X_Axis;
    "set grid on"
*    if (irow<Nrow-1)
    "set t 1 "Len; "set y 1 "Ntime
    T_clor()
    "d spi"
    draw_tit_ver(x1-0.3, (y1+y2)/2, R.ireg)
  Endwhile
  "xcbar "x00" "x01-0.2" 0.3 0.45 -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on "
  draw_tit((x00+x01)/2, y00+0.2,"Spatial drought extent (Black: Moderate, Yellow: Severe, Red: Extreme)")
*  draw_tit((x00+x01)/2, y00+0.1,"and areal mean (blue) of drought index")
  "gxprint " Out_Dir "Time_evolution_spi.pdf pdf"
  "set parea off";
  "close 1"

*******************************************************
function draw_tit_ver(x, y, title)
  "set string 1 c 5.5 90"; "set strsiz 0.12 0.12"
  "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
return
*******************************************************
function T_clor()
levels = "  -2.0 -1.5 -1 -0.5  0  0.5  1 1.5 2 "
"run color.gs -gxout shaded -kind darkgreen->green->lightgreen->white-(0)->white->orange->red->darkred -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 color.gs -gxout grfill -kind red->yellow-(0)->aqua->blue -levs" levels
"set grads off"
return
*******************************************************
function Line_pl(mark,colo,thick,styl)
"set gxout line"
"set cmark "mark
"set ccolor "colo
"set cthick "thick
"set cstyle "styl
return
********************************************
function Mark_leg(colo,type,x,y,size)
"set line "colo
"draw mark "type" "x" "y" "size
return
********************************************
function Truc_0(a0)
"set gxout line"
"set ccolor 1"
"set cstyle 1"
"set cthick 2"
"set cmark 0"
"d "a0
return

********************************************
function ENSO()
*    "sdfopen /home/tanpv/Progs/Dr_SEA/EOF_New/ONI_ENSO.nc"
*     "sdfopen /home/tanpv/Progs/ENSO_Years/ONI_ENSO_1951-2020.nc"
    "sdfopen /home/tanpv/Progs/ENSO_Years/ONI_ENSO_1961-2005_.nc"

    "set x 13 540"
*    "set x 13 540"
    "set gxout linefill"; "set ylpos 0 l"; "set vrange -2 2"; "set ylint 1";
    "set y 1"
    "define a1=oni-oni-0.5"; "define a0=oni-oni"
    "set lfcols 15 15"; "d oni; a0"
    "set lfcols 0 4"; "d oni; a1";
    "set y 3"
    "define a1=oni-oni+0.5";
    "set lfcols 15 15"; "d oni; a0"
    "set lfcols 2 0"; "d oni; a1";
    "set y 2";
    "define a1=oni-oni-0.5"; "define a2=oni-oni+0.5";
    "set lfcols 15 15";
    "d oni; a0";
    "set y 4"
    "set gxout line"; "set line 4 1 0.1"; "set ccolor 1"; "set cmark 0"; 
    "set ylpos 0 l"; "d oni";
return
*******************************************************