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

   Mon.1="01"; Mon.2="02"; Mon.3="03";  Mon.4="04"; Mon.5="05";  Mon.6="06"
   Mon.7="07"; Mon.8="08"; Mon.9="09"; Mon.10="10"; Mon.11="11"; Mon.12="12"

   Dir_Inp="/home/tanpv/Data/562_VNU/0.25_VN/Dr_IDXs_Ope/"
   Dat_out="/home/tanpv/Progs/562_VNU/2024/PDF/"

   Nidx=2        ;* spi, spei
   Nlen=3        ;* 1,3,6 months
   NDr_lev=3     ;* mild, moderate, extreme
   Nlead=6
   Nvar=4        ;* pc, pod, ts, hss
   Nmon=12;

   Var.1="pc"; Var.2="pod"; Var.3="ts"; Var.4="hss"
   Didx.1="spi"; Didx.2="spei"
   Len.1="01"; Len.2="03"; Len.3="06"
   Lev.1="Mild"; Lev.2="Moderate"; Lev.3="Extreme"

   Ncol=Nmon; Nrow=Nvar
   
   y00=8.0; y01=0.3
   dy=(y00-y01)/Nrow+0.09
   x00=0.3; x01=11
   dx=(x01-x00)/Ncol+0.12
   yy2=(y00+y01)/2+0.05

   idx=0
   While (idx<Nidx)    ;* =2
     idx=idx+1; ilen=0
     While (ilen<Nlen) ;* =3
       ilen=ilen+1; 
       "sdfopen "Dir_Inp "Veri_Dr_IDX_1983-2020_" Didx.idx "_" Len.ilen ".nc"
       ilead=0
       While (ilead<Nlead)
         ilead=ilead+1; ilev=0
         While (ilev<NDr_lev)    ;* mild, moderate, extreme
           ilev=ilev+1; icol=0; 
           While (icol<Ncol)
             icol=icol+1; irow=0
             if (icol=1); x1=x00; x2=x1+dx; x11=x2; endif
             if (icol>1); x1=x11-0.2; x2=x1+dx; x11=x2; endif
             While (irow<Nrow)
               irow=irow+1
               if (irow=1); y2=y00; y1=y2-dy; y11=y1; endif
               if (irow>1); y2=y11+0.2; y1=y2-dy; y11=y1; endif
               "set parea "x1" "x2" "y1" "y2
               "set mpdraw off"; "set frame off"; "set xlab off";"set ylab off"

               "set t " ilead; "set z " ilev
*               "d "Var.irow"_"Mon.icol
               if (irow<Nrow); prob_Clor(); "d "Var.irow"_"Mon.icol; endif
               if (irow=Nrow); ss_Clor(); "d "Var.irow"_"Mon.icol; endif
*               "d "Var.irow"_"Mon.icol
               if (irow=1); draw_tit((x1+x2)/2-0.1, y2-0.05, "Month "Mon.icol); endif
               if (icol=1); draw_ver_tit(x1+0.1, (y1+y2)/2, Var.irow); endif
               if (icol=Ncol); "xcbar "x2-0.1" "x2+0.1" "y1+0.1" "y2-0.2" -fs 2 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on "; endif
             Endwhile
           Endwhile
           Main_tit((x00+x01)/2, y00+0.2, "Drought IDX: "Didx.idx" - Timescale="Len.ilen" month(s) - Lead="ilead" month(s)- Drought level: "Lev.ilev)
           pull dd
           "c"
         Endwhile
       Endwhile
       "close 1"
     Endwhile
   Endwhile
*******************************************************
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 ss_Clor()
*   levels = "  -2 -1.5 -1  0  1 1.5 2  " 
*   "run color.gs -gxout grfill -kind darkred->red->yellow->white-(0)->white->cyan->blue->darkblue -levs " levels
   levels = " -0.3 -0.2 -0.1 -0.05 0 0.05 0.1 0.2 0.3 " 
   "run color.gs -gxout grfill -kind green->lightgreen-(0)->yellow->red -levs " levels

   "set grads off"
   return

********************************************
   function prob_Clor()
*   levels = "  -2 -1.5 -1  0  1 1.5 2  " 
*   "run color.gs -gxout grfill -kind darkred->red->yellow->white-(0)->white->cyan->blue->darkblue -levs " levels
*   levels = " 0.01 0.02 0.03 0.05 0.1 0.2 0.3 0.5 0.6 0.7 0.8 0.9 " 
   levels = " 0.05 0.1 0.2 0.3 0.5 0.6 0.7 0.8 0.9 " 
   ;*levels = "   -1.0  0  1.0  "
   "run color.gs -gxout grfill -kind yellow->orange->red -levs " levels

   "set grads off"
   return
********************************************
   function T_Clor()
   levels = " 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 " 
   "run color.gs -gxout grfill -kind lightgreen->yellow->orange->red -levs " levels
   "set grads off"
   return
********************************************
   function Pr_Clor()
*   levels = " 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 " 
   levels = "  0.8 0.85 0.9 0.95 "
   "run color.gs -gxout grfill -kind white->yellow->orange->red->darkred -levs " levels
   "set grads off"
   return
********************************************
   function R_Clor()
*   levels = " 10 30 50 70 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 "
*   "run color.gs -gxout grfill -kind lightgreen->yellow->orange->red -levs " levels
    levels = " 5 10 20 30 50 70 100 150 200 250 300 350 400 450 500 550 600 650 700 "
    "run color.gs -gxout grfill -kind green->lightgreen->yellow->orange->red -levs " levels
*    "run color.gs -gxout grfill -kind red->orange->yellow->lightgreen->green -levs " levels
   "set grads off"
   return