"reinit"
   "set display color white"; "c"
   Inp_Dir="/home/tanpv/Progs/S2S_Vn_Ope/Res/"

   Ncol=3; Nrow=1      ;* model in row, char in col 
   V.1="rme"; V.2="rmae"; V.3="cor"
   T.1="RME(BC)-RME(ORG)"; T.2="RMAE(BC)-RMAE(ORG)"; T.3="CORR(BC)-CORR(ORG)"

   Nseason=4
   Season.1="12-02"; Season.2="03-05"; Season.3="06-08"; Season.4="09-11"
   y00=8.0; y01=1.0
   dy=(y00-y01)/Nrow
   x00=0.7; x01=11
   dx=(x01-x00)/Ncol-0.3

   ;* "sdfopen " Inp_Dir "Veri_Score_R6_BC_test.nc"
   ;* "sdfopen " Inp_Dir "Veri_Score_12-02_corrected.nc"
   "set mpdraw on"
   "set mpdset hires"
   "set grid off"
   ;* "set lat 11 15.7"
   ;* "set lon 106.3 109.2"
   Nt=44; delt=10
 isea=0
 While (isea<Nseason)
   isea=isea+1
   "sdfopen " Inp_Dir "Veri_Score_" Season.isea "_corrected.nc"
   "sdfopen " Inp_Dir "Veri_Score_" Season.isea "_org.nc"
   it=0
   While (it<Nt-delt)
     it=it+delt
     irow=0    ;* model
     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; ;* char
       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 2"
         if (icol=1); "set ylab on"; "set ylint 2"; endif
         Colo(icol) 
         "d "V.icol".1(t="it")-" V.icol".2(t="it")"
         draw_tit((x1+x2)/2,y2-0.3, T.icol " for LT="it" days")
         Sub_Reg2
         "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 "

       Endwhile
     Endwhile
     draw_tit((x00+x01)/2, y00+0.1, "Differences between BC and ORG for season "Season.isea)
     pull dd; "c"
   Endwhile
   "close 2"
   "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 = " -80 -60 -40 -20 -10 0 10 20 40 60 80 "; endif
*     if (cha=3); levels = " -0.8 -0.6 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.6 0.8"; endif 
     if (cha=3); levels = " -0.5 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.5 "; 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
     "run color.gs -gxout grfill -kind green->lightgreen-(0)->yellow->red -levs " levels
   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"
*    "set ccolor 4"
    "d "var
    "set grads off"
return
********************************************