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

   Ncol=2; Nrow=1      ;* model in row, char in col 
   V.1="me"; V.2="mae"
   T.1="ME"; T.2="MAE"

   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_Org_Absolute_Errors_R6.nc"

   "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
     
     if (count_valid = Ncol)
       "printim "Out_Dir"Veri_Org_me_mae_t"it".png png white"
     endif
     
     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 = " 30 40 50 60 80 100 120 140 "; 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 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