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

   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"

   Dir_Inp="/home/tanpv/Progs/Dr_Daily/Res/"
   Dat_out="//home/tanpv/Progs/Dr_Daily/PDF/"

   Nvar=4; Nyear=60

   V.1="frq"; V.2="ndy"; V.3="dur"; V.4="dmx"
   Var.1="Annual Number of CDH"; Var.2="Annual Number of the CDH days" 
   Var.3="Annual mean duration of CDH"; Var.4="Annual Maximum duration of CDH"
   Ncol=2; Nrow=2 
   y00=8.0; y01=0.2
   dy=(y00-y01)/Nrow-0.4
   x00=0.5; x01=10.5
   dx=(x01-x00)/Ncol-0.0

   "sdfopen "Dir_Inp "CDH_Chars_1.nc"
   iyr=0
   While (iyr<Nyear)
     iyr=iyr+1; 
       irow=0; ivar=0 
       While (irow<Nrow)
         irow=irow+1; icol=0 
         if (irow=1); y2=y00; y1=y2-dy; y11=y1; endif
         if (irow>1); y2=y11-0.3; y1=y2-dy; y11=y1; endif
         While (icol<Ncol)
            icol=icol+1; ivar=ivar+1 
            if (icol=1); x1=x00; x2=x1+dx; x11=x2; endif
            if (icol>1); x1=x11+0.02; x2=x1+dx; x11=x2; endif
            "set parea "x1" "x2" "y1" "y2
            "set mpdraw on"; "set frame on"; "set xlab off"; "set grid off"; "set ylab off";"set xlint 5"; "set ylint 5"
           ;* "set lon 101 111"; "set lat 7.5 24.0"
            if (icol=1); "set ylab on"; "set grid off";  endif
            if (irow=Nrow); "set xlab on"; "set grid off";  endif

            "set t " iyr
            Char_Clor(ivar)
            "d "V.ivar"(t="iyr")"
            draw_tit((x1+x2)/2, y2+0.1, Var.ivar);
            "xcbar "x2-0.25" "x2-0.05" "y1+0.4" "y2-0.4" -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on "
            ;* Sub_Reg2
            ;* "Title_lonlat -lon " 109 " -lat "  17.0 " -siz " 0.10 " -ang "300 " -tit "Bien_Dong""
         Endwhile
         if (irow=1)
            Main_tit((x00+x01)/2,y2+0.25, "CDH characteristics of year "iyr+1960)
         endif
       Endwhile
       ;* "gxprint " Dat_out V.ivar"_forecast.png png"
       
     pull dd
     "c"

   Endwhile
   "close 1"
   "quit"
*******************************************************
function Main_tit(x, y, title)
  "set string 1 c 5.5 0"; "set strsiz 0.11 0.13"
  "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 Char_Clor(var)
   if (var=1); levels = " 1 2 3 4 5 6 7 8 9 10 "; endif              ;* frq
   if (var=2); levels = " 5 10 15 20 25 30 35 40 50 "; endif  ;* ndy
   if (var=3); levels = " 2 3 5 7 10 15 20 "; endif    ;* dur
   if (var=4); levels = " 5 10 15 20 25 30 "; endif ;* dmx

   "run color.gs -gxout grfill -kind lightyellow->orange->red -levs " levels
   "set grads off"
   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 Con_pr()
   levels = " 5 10 20 30 50 70 100 150 200 250 300 350 400 450 500 "
   "set clevs " levels
   "set ccolor 4"
   "set clab off"
   "set gxout contour"
   "set grads off"
   return
********************************************
   function Pr_Clor()
*   "run color60_4rainfall.gs"
*    levels = " 5 10 20 30 40 50 60 70 80 90 " 
   levels = "  5 10 20 30 50 70 100 120 140 160 180 200 240 280 320 350 400 450 500 "
*    "run color.gs -gxout grfill -kind green->lightgreen->yellow->orange->red -levs " levels
*   "run color.gs -gxout grfill -kind lightblue->blue->darkblue -levs " levels

   colors = " 32 36 40 48 50 52  56  58  60  62  64  66  68  70  74  76  80  84  88 90 "
   "set clevs " levels
   "set ccols " colors
   "set gxout grfill"
   "run color60_4rainfall.gs"
   "set grads off"
   return
********************************************
   function Con_pr_ano()
   levels = " -100 -80 -60 -40 -30 -20 -10 -5 0 5 10 20 30 40 60 80 100 " 
   "set clevs " levels
   "set ccolor 4"
   "set clab off"
   "set gxout contour"
   "set grads off"
   return
********************************************
   function Pr_ano()
   levels = "-100 -80 -60 -40 -30 -20 -10 -5 0 5 10 20 30 40 60 80 100 "
   "run color.gs -gxout grfill -kind red->orange->yellow->white-(0)->white->lightblue->blue->darkblue -levs " levels
   return
********************************************
   function Con_t2()
   levels = " 10 12 14 16 18 20 22 24 26 28 30 32 "
   "set clevs " levels
   "set ccolor 4"
   "set clab off"
   "set gxout contour"
   "set grads off"
   return
********************************************
   function t2_Clor()
   levels = " 10 12 14 16 18 20 22 24 26 28 30 32 "
   "run color.gs -gxout grfill -kind darkgreen->green->lightgreen->yellow->orange->red->darkred -levs " levels
   "set grads off"
   return
********************************************
   function Con_t2_ano()
   levels = " -2 -1.5 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.5 2 " 
   "set clevs " levels
   "set ccolor 4"
   "set clab off"
   "set gxout contour"
   "set grads off"
   return
********************************************
   function t2_ano()
   levels = " -2 -1.5 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.5 2 " 
   "run color.gs -gxout grfill -kind darkgreen->green->lightgreen-(0)->yellow->orange->red -levs " levels
   "set grads off"
   return
********************************************