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

   Ncol=2; Nrow=2      ;* model in row, char in col 
   V.1.1="sshf"; V.1.2="slhf"
   V.2.1="ttr";  V.2.2="tp"

   y00=10.5; y01=1.0
   dy=(y00-y01)/Nrow
   x00=0.3; x01=7.8
   dx=(x01-x00)/Ncol

   "sdfopen " Inp_Dir "LH_SH_OLR_TP_cf_2008-04-29.nc"

   "set mpdraw on"
   "set mpdset hires"
   Nt=46
   it=0
   While (it<Nt)
     it=it+1
     ivar=0
     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; ivar=ivar+1
         if (icol=1); x1=x00; x2=x1+dx; x11=x2; endif
         if (icol>1); x1=x11+0.3; x2=x1+dx; x11=x2; endif

         "set parea "x1" "x2" "y1" "y2
         Colo(ivar) 
         if (ivar<=3); "d "V.irow.icol"(t="it")";endif
         if (ivar=4) ; "d "V.irow.icol"(t="it")";endif
         draw_tit((x1+x2)/2,y2+0.1, V.irow.icol " - t="it)
         Sub_Reg2
         "xcbar "x2-0.4" "x2-0.2" "y1+0.1" "y2-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.09 0.10"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return

********************************************
   function Colo(cha)
     if (cha=1); levels = " 10 20 30 40 50 60 70 80 90 100 "; endif 
     if (cha=2); levels = " 30 40 50 60 70 80 90 100 110 120 130 "; endif
     if (cha=3); levels = " 140 160 180 200 220 240 260 280 300 320 "; endif 
     if (cha=4); levels = " 2 3 4 5 6 8 10 15 20 25 30 35 40 50 "; endif 
     "run color.gs -gxout grfill -kind green->lightgreen->yellow->red -levs " levels
   return
********************************************