"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"; T.2="RMAE"; T.3="CORR"

   

   y00=8.0; y01=0.5
   dy=0.3
   x00=0.5; x01=9.5
   dx=(x01-x00)/Ncol-0.3

   "sdfopen " Inp_Dir "Veri_Org_Absolute_Errors_R6.nc" ;*"Veri_Score_R6_ME_MAE_BC_Ann_new.nc" ;*"Veri_Score_R6_ME_MAE.nc"
   "set t 3 44"
   Nst=12
   RMAE_1()
   "set vrange 5 35"; "set ylint 5"
   "set parea "x00" "x01" "y01" "y00
   x1=x01+0.1; x2=x1+0.5
   y1=y00-(y00-y01)/4
   "set xaxis 3 44 2"; "set grads off"
  ;* RMAE
   ist=0
   While (ist<Nst)
     ist=ist+1
     Line_pl(0,ist,5,ist)   ;* mark,colo,thick,styl
     "d S"ist
     y1=y1-dy; y2=y1
     Legend(ist,ist,5,x1,y1,x2,y2)  ;* colo,type,thick
     draw_tit_left(x2+0.05, y1, _Stn.ist)
   Endwhile
   pull dd
   "c"
   "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 draw_tit_left(x, y, title)
  "set string 1 l 5.5 0"; "set strsiz 0.1 0.13"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return
********************************************
function RMAE_1()
   "set lon 107.833"; "set lat 14.650"; _Stn.1="DAKTO"; "define S1=mae"
   "set lon 108.000"; "set lat 14.333"; _Stn.2="KONTUM";"define S2=mae"
   "set lon 108.017"; "set lat 13.967"; _Stn.3="PLEIKU";"define S3=mae"
   "set lon 108.650"; "set lat 13.950"; _Stn.4="ANKHE" ;"define S4=mae"
   "set lon 108.450"; "set lat 13.383"; _Stn.5="AYUNPA" ;"define S5=mae"
   "set lon 108.267"; "set lat 12.917"; _Stn.6="BUONHO";"define S6=mae"
   "set lon 108.767"; "set lat 12.733"; _Stn.7="MDRAK"; "define S7=mae"
   "set lon 108.050"; "set lat 12.667"; _Stn.8="BMTHUOT";"define S8=mae"
   "set lon 107.617"; "set lat 12.450"; _Stn.9="DAKMIL";"define S9=mae"
   "set lon 107.683"; "set lat 12.000"; _Stn.10="DAKNONG";"define S10=mae"
   "set lon 108.450"; "set lat 11.950"; _Stn.11="DALAT";"define S11=mae"
   "set lon 107.817"; "set lat 11.533"; _Stn.12="BAOLOC";"define S12=mae"
return 
********************************************
function COR_1()
   "set lon 107.833"; "set lat 14.650"; _Stn.1="DAKTO"; "define S1=cor"
   "set lon 108.000"; "set lat 14.333"; _Stn.2="KONTUM";"define S2=cor"
   "set lon 108.017"; "set lat 13.967"; _Stn.3="PLEIKU";"define S3=cor"
   "set lon 108.650"; "set lat 13.950"; _Stn.4="ANKHE" ;"define S4=cor"
   "set lon 108.450"; "set lat 13.383"; _Stn.5="AYUNPA" ;"define S5=cor"
   "set lon 108.267"; "set lat 12.917"; _Stn.6="BUONHO";"define S6=cor"
   "set lon 108.767"; "set lat 12.733"; _Stn.7="MDRAK"; "define S7=cor"
   "set lon 108.050"; "set lat 12.667"; _Stn.8="BMTHUOT";"define S8=cor"
   "set lon 107.617"; "set lat 12.450"; _Stn.9="DAKMIL";"define S9=cor"
   "set lon 107.683"; "set lat 12.000"; _Stn.10="DAKNONG";"define S10=cor"
   "set lon 108.450"; "set lat 11.950"; _Stn.11="DALAT";"define S11=cor"
   "set lon 107.817"; "set lat 11.533"; _Stn.12="BAOLOC";"define S12=cor"
return
*******************************************************
function Line_pl(mark,colo,thick,styl)
"set gxout line"
"set cmark "mark
"set ccolor "colo
"set cthick "thick
"set cstyle "styl
return
********************************************
function Legend(colo,type,thick,x1,y1,x2,y2)
"set line "colo " "type" "thick 
"draw line "x1" "y1" "x2" "y2
return
********************************************
   function Colo(cha)
     if (cha=1); levels = "  -20 -10 -5 -3 -2 -1 0 1 2 3 5 10 20 "; endif 
     if (cha=2); levels = " 1 2 3 5 7 10 15 20 30 "; 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.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"
*    "set ccolor 4"
    "d "var
    "set grads off"
return
********************************************