"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_SEA/Post_APH/PDF1/"
 In_Dir="/home/tanpv/Progs/Dr_SEA/Post_APH/NC1/"
 X_Axis="1962 2005"
 idx=3

  File.1="eof_APH.nc"
  Geo.1="APH_geo_reg_Mod.nc"
  Geo.2="APH_geo_reg_Sev.nc"
  Geo.3="APH_geo_reg_Ext.nc"
  Geo.4="Areal_mean_and_Trend_APH_spei_12.nc"

  GeoSEA.1="APH_geo_Mod.nc"
  GeoSEA.2="APH_geo_Sev.nc"
  GeoSEA.3="APH_geo_Ext.nc"
  GeoSEA.4="Areal_mean_and_Trend_APH_spei_12.nc"

 EOVar.1="spe12"

 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=1.0
 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=11;
   x1=x00+0.0; x2=x01-0.0
*   "set x 1 "Ntime
   dy=(y00-y01)/(Nrow+1)-0.05
   
   "sdfopen " In_Dir File.1
   "query dims"; str=sublin(result,2);  Nmons=subwrd(str,8)
   Ntime=Nmons; "set x 1 "Ntime
   "set gxout linefill"
   irow=0; it=0; iplt=0
   While (iplt<Nmode)
         irow=irow+1; iplt=iplt+1
         if (irow=1); y2=y00; y1=y2-dy; y11=y1; y12=y2; endif
         if (irow>1); y2=y11-0.09; y1=y2-dy; y11=y1;y12=y2; endif
         "set parea "x1" "x2" "y1" "y2
         "set string 1 c 5.5 0"; "set strsiz 0.02 0.02"
         "set vrange -0.12 0.12"; 
         "set ylint 0.05";
         "set xlab off"; "set ylab on";
         if (iplt=Nmode); "set xlab on"; "set xlint 5"; "set strsiz 0.11 0.11"; "set xaxis "X_Axis ; endif
         it=it+1
         "set y "it
         "define aa=" EOVar.1
         "define a0=aa-aa"; 
         "set lfcols 5 0"; "d aa; a0";
         "set lfcols 5 8"; "d aa; a0";
         Line_pl(0,1,3,1)
         "d aa"
         "set gxout linefill"
         draw_tit(x00+0.3,y1+0.1, "EOF" it)
         if (irow=1); draw_tit(x00+(x01-x00)/2, y2+0.15, "a) Time series of EOFs of SPEI_12"); endif
   Endwhile
   "close 1"

   Nreg=5
   R.1="R1";R.2="R2";R.3="R3";R.4="R4";R.5="R5";R.6="SEA"
   y00=y1-0.7; y01=1
   Nrow=Nreg+2; NDr_type=3+1; 

  irow=0
  While (irow<Nrow)
    irow=irow+1; ireg=irow
    if (irow=1); y2=y00-(irow-1)*dy; y1=y2-dy; y11=y1; endif
    if (irow>1); y2=y11-0.1; y1=y2-dy; y11=y1; endif
    "set xlab off"; "set ylab on"; "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
    if (irow<Nrow-1)
*    if (irow<Nrow)
    ityp=0
    While (ityp<NDr_type)
      ityp=ityp+1
      "sdfopen " In_Dir Geo.ityp
      if (ityp<NDr_type)
         "set gxout linefill"
         "set t "idx
         "set y "irow;
         "query dims"; str=sublin(result,2);  Nmons=subwrd(str,8)
         "set x 1 "Nmons
         "set ylpos 0 l"; "set vrange 0 100"; "set ylint 20"
*         Line_pl(mark.ityp,colo.ityp,thic.ityp,styl.ityp)
*         "set lfcols "colo.ityp" 0"
         if (ityp=1); "set lfcols 1 0"; endif
         if (ityp=2); "set lfcols 7 0"; endif
         if (ityp=3); "set lfcols 2 0"; endif
         "define a0=georeg-georeg"
         "d georeg; a0"
         "close 1"
       endif
       if (ityp=NDr_type)
         "set x "irow;
         "query dims"; str=sublin(result,3);  Nmons=subwrd(str,8)
         "set y 1 "Nmons
         "set ylpos 0 r"; "set vrange -3 3"; "set ylint 1"
         Line_pl(mark.ityp,colo.ityp,thic.ityp,styl.ityp)
         "d dridx"
         "close 1"
       endif
    Endwhile
    draw_tit((x1+x2)/2, y2-0.1, R.ireg)
    endif
*   ------------- Entire SEA:
    if (irow=Nrow-1)
*    if (irow=Nrow)
    ityp=0
    While (ityp<NDr_type)
      ityp=ityp+1
      "sdfopen " In_Dir GeoSEA.ityp
      if (ityp<NDr_type)
         "set gxout linefill"
*         "set t "idx
         "set y "idx;
         "query dims"; str=sublin(result,2);  Nmons=subwrd(str,8)
         "set x 1 "Nmons
         "set ylpos 0 l"; "set vrange 0 100"; "set ylint 20"
*         Line_pl(mark.ityp,colo.ityp,thic.ityp,styl.ityp)
*         "set lfcols "colo.ityp" 0"
         if (ityp=1); "set lfcols 1 0"; endif
         if (ityp=2); "set lfcols 7 0"; endif
         if (ityp=3); "set lfcols 2 0"; endif
         "define a0=geo-geo"
         "d geo; a0"
         "close 1"
       endif
       if (ityp=NDr_type)
         "set x "Nrow-1;
         "query dims"; str=sublin(result,3);  Nmons=subwrd(str,8)
         "set y 1 "Nmons
         "set ylpos 0 r"; "set vrange -3 3"; "set ylint 1"
         Line_pl(mark.ityp,colo.ityp,thic.ityp,styl.ityp)
         "d dridx"
         "close 1"
       endif
    Endwhile
    draw_tit((x1+x2)/2, y2-0.1, R.ireg)
    endif
*   ---------------- Moisture:
    if (irow=Nrow)
       "sdfopen " In_Dir "Moisture_flux.nc"
       "set gxout linefill"
       "set y 3"
       "query dims"; str=sublin(result,2);  Nmons=subwrd(str,13)
       "define a0=qflux-qflux"
       "define aa=-qflux/1.0E+11"
       "set ylpos 0 l"; "set vrange -1.5 1.5"; "set ylint 0.5"
       "set lfcols 4 0"
       "d aa; a0"
       "set lfcols 4 2"
       "d aa; a0"
       "close 1"
    endif
*    "set ylab off"
  Endwhile
  draw_tit((x00+x01)/2, y00+0.3,"b) Spatial drought extent (Blue: Moderate, Yellow: Severe, Red: Extreme)")
  draw_tit((x00+x01)/2, y00+0.1,"and areal mean (blue) of drought index")
      y2=y1-0.5; y1=y2-dy 
      "set parea "x1" "x2" "y1" "y2
      ENSO()
      draw_tit(x00+(x01-x00)/2,y2+0.1, "c) Time serie of ONI")
   "gxprint " Out_Dir "Timeseries_eof_Dr_Geo.pdf pdf"
   "set parea off";
   "close 1"

*******************************************************
function draw_tit_ver(x, y, title)
  "set string 1 c 5.5 270"; "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 = " -0.6 -0.4 -0.2 -0.1 0 0.1 0.2 0.4 0.6 "

*levels = " -1 -0.8 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.8 1"
levels = " -1 -0.8 -0.6 -0.4 -0.2 -0.1 0 0.1 0.2 0.4 0.6 0.8 1"
*levels = " -10 -8 -6 -5 -4 -3 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 3 4 5 6 8 10"
*levels = " -10 -8 -6 -5 -4 -3 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 3 4 5 6 8 10"
*levels = " -30 -20 -10 -8 -6 -5 -4 -3 -2 -1.5 -1 0 1 1.5 2 3 4 5 6 8 10 20 30"
*"run color.gs -gxout grfill -kind blue->cyan->white-(0)->white->red->darkred -levs" levels
"run color.gs -gxout grfill -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
*******************************************************