"reinit"
   "set display color white"; "c"
   Inp_Dir="/home/cordex_sea/Chi/Res/"

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

   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.nc"
   "sdfopen " Inp_Dir "Veri_Score_R6_BC.nc"

   Va.1="rme.1"; Va.2="rmae.1"; Va.3="cor.1"
   Vb.1="rme.2"; Vb.2="rmae.2"; Vb.3="cor.2"

   "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
     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 0.5"
         if (icol=1); "set ylab on"; "set ylint 0.5"; endif
         Colo(icol) 
         if (icol=1); "d "Vb.icol"(t="it")"; endif                    ;* RME of BC case
         if (icol>1); "d "Vb.icol"(t="it")-"Va.icol"(t="it")"; endif  ;* Differences between BC and ORG
         draw_tit((x1+x2)/2,y2-0.3, T.icol " for LT="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 "

       Endwhile
     Endwhile
     pull dd; "c"
   Endwhile
   "close 1"
     
*******************************************************
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 = " -40 -30 -20 -10 0 10 20 30 40 "; endif
     if (cha=3); levels = " -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 "; 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); "run color.gs -gxout grfill -kind red->yellow-(0)->lightgreen->green -levs " levels ; endif
     if (cha=3); "run color.gs -gxout grfill -kind green->lightgreen-(0)->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"
*    "set ccolor 4"
    "d "var
    "set grads off"
return
********************************************