"reinit"
   "set display color white"; "c"
   Inp_Dir="/home/tanpv/Progs/S2S_4Giang/Res/"
   Out_Dir="/home/tanpv/Progs/S2S_4Giang/PDF/"

   ;* Veri_Score_03-05.nc
   ;* Veri_Score_09-11_corrected.nc
   Nseason=4
   Seas.1="12-02"; Seas.2="03-05"; Seas.3="06-08"; Seas.4="09-11"
   Ncol=3; Nrow=1      ;* model in row, char in col 
   V.1="rme"; V.2="rmae"; V.3="cor"
   T.1="RME"; T.2="RMAE"; T.3="CORR"

   y00=8.0; y01=1.0
   dy=(y00-y01)/Nrow
   x00=0.7; x01=11
   dx=(x01-x00)/Ncol-0.3

 isea=0
 While (isea<Nseason)
   isea=isea+1
  ;* "sdfopen " Inp_Dir "Veri_Score_" Seas.isea ".nc"   ;* Origin
   "sdfopen " Inp_Dir "Veri_Score_" Seas.isea "_corrected.nc"     ;* BC
   "set mpdraw on"
   "set mpdset hires"
   ;*"set lat 11 15.7"
   ;*"set lon 106.3 109.2"
   Nt=44; delt=5
   it=delt+2
   While (it<Nt-delt)
     it=it+delt
     irow=0    ;* model
     valid_plot = 1
     count_valid = 0
     While (irow<Nrow)
       irow=irow+1
       if (irow=1); y2=y00; y1=y2-dy; y11=y1; endif
       if (irow>1); y2=y11-0.4; y1=y2-dy; y11=y1; endif

       icol=0; 
       While (icol<Ncol)
         icol=icol+1; 
         if (icol=1); x1=x00; x2=x1+dx; x11=x2; endif
         if (icol>1); x1=x11+0.2; x2=x1+dx; x11=x2; endif

         "set parea "x1" "x2" "y1" "y2
         "set ylab off"
         "set xlab on"; "set xlint 0.5"
         if (icol=1); "set ylab on"; "set ylint 0.5"; endif
         Colo(icol) 
         "d "V.icol"(t="it")"
         draw_tit((x1+x2)/2,y2-0.3, T.icol " for Lead time="it" days")
         Sub_R6
         "xcbar "x1+0.0" "x2-0.2" "y1-0.1" "y1+0.1" -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on "
         
         count_valid = count_valid + 1
       Endwhile
     Endwhile
     draw_tit((x00+x01)/2, y00, "Season: "Seas.isea) 
     if (count_valid = Ncol)
       ;*"printim "Out_Dir"Veri_Org_rm_t"it".png png white"
     endif
     
     pull dd; "c"
   Endwhile
   "close 1"
 Endwhile    
******************************************************
function draw_tit(x, y, title)
  "set string 1 c 5.5 0"; "set strsiz 0.1 0.13"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return

********************************************
function Colo(cha)
    if (cha=1); levels = "  -80 -60 -40 -20 -10 0 10 20 40 60 80 "; endif 
    if (cha=2); levels = " 30 40 50 60 80 100 120 140 "; endif
    if (cha=3); levels = " 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 "; endif 
    if (cha=1); "run color.gs -gxout grfill -kind green->lightgreen-(0)->yellow->red -levs " levels ; endif
    if (cha=2 | cha=3); "run color.gs -gxout grfill -kind green->lightgreen->yellow->red -levs " levels ; endif
return
********************************************
function d_Stn(var)
    "run color.gs -gxout shaded -kind blue->blue -levs " 1 500
    "set gxout stnmark"
    "set digsiz 0.2"
    "set cmark 3"
    "set cthick 7"
    "d "var
    "set grads off"
return