"reinit"
   "set display color white"; "c"

   Mon.1="Jan"; Mon.2="Feb"; Mon.3="Mar"; Mon.4="Apr"; Mon.5="May"; Mon.6="Jun"
   Mon.7="Jul"; Mon.8="Aug"; Mon.9="Sep"; Mon.10="Oct"; Mon.11="Nov"; Mon.12="Dec"

   L.1="LT=1 mon"; L.2="LT=2 mon"; L.3="LT=3 mon"; L.4="LT=4 mon"; L.5="LT=5 mon"; L.6="LT=6 mon";

   Var.1="Relative Mean Absolute Error (%) - ANN correction"; Var.2="Correlation between ANN correction and observation"; Var.3="Added value of ANN correction over the model output"; 
   Dat_out="/home/tanpv/Progs/562_VNU/PDF/ANN/"

   Score.1="rmae"; Score.2="cor"; Score.3="av"

*  3 scores; 12 months; 5 leadtimes
   Nvar=3; Ncol=12; Nrow=6 
   y00=8.0; y01=0.3
   dy=(y00-y01)/Nrow+0.1
   x00=0.5; x01=10.5
   dx=(x01-x00)/Ncol+0.1

   Model="clWRF"
*   Model="RegCM"
   Nmem=6
   MEM=0
While (MEM<Nmem)
   MEM=MEM+1
   "sdfopen /home/tanpv/Progs/562_VNU/test/Verification_scores_" Model "_mem_" MEM ".nc"

 ivar=0
 While (ivar<Nvar)
   ivar=ivar+1; irow=0
   While (irow<Nrow)
     irow=irow+1
     if (irow=1); y2=y00; y1=y2-dy; y11=y1; endif
     if (irow>1); y2=y11+0.08; y1=y2-dy; y11=y1; endif

     icol=0; 
     While (icol<Ncol)
       icol=icol+1; ilead=irow-1
       if (icol=1); x1=x00; x2=x1+dx; x11=x2; endif
       if (icol>1); x1=x11-0.15; x2=x1+dx; x11=x2; endif

       "set parea "x1" "x2" "y1" "y2
       "set mpdraw off"; "set frame off"; "set xlab off";"set ylab off"
       "set t "icol
       if (ivar=1); 
          R_rmae(); "d "Score.ivar"*100"; 
       endif
       if (ivar=2);
          Cor(); "d "Score.ivar;
*          "hatch pval 0 0.05 -angle 45,135 -int 0.05 -color 1 -thickness 0.01 -type 1 -density 0.01" 
*         draw_ver_tit(x00-0.1, (y1+y2)/2, L.irow); endif
       endif
       if (ivar=3); 
          Cor(); "d "Score.ivar; 
*          if (icol=1); draw_ver_tit(x00-0.1, (y1+y2)/2, L.irow); endif
       endif
       if (icol=1); draw_ver_tit(x00-0.0, (y1+y2)/2, L.irow); endif
       if (irow=1); draw_tit((x1+x2)/2-0.05, y00+0.0, Mon.icol); endif
       Sub_Reg2

*       HSTS
*       "Title_lonlat -lon " 112.8 " -lat " 14.500 " -siz " 0.07 " -ang "0 " -tit "HSa""  
*       "Title_lonlat -lon " 113.5 " -lat "  8.250 " -siz " 0.07 " -ang "0 " -tit "TSa""         

*       "hatch pv 0 0.05 -angle 135 -int 0.05 -color 1 -thickness 0.01 -type 1 -density 0.01"
*       draw_tit(x2-0.55, y2-0.2, Varn.ivar)
     Endwhile
   Endwhile
   "xcbar "x01-0.45" "x01-0.30" "y01+0.5" "y00-0.5" -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on "
*   "xcbar "x00" "x01-0.5" 0.3 0.45 -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on -edge triangle "
*   "xcbar "x01-0.45" "x01-0.25" "y01-0.2" "y00-dy-0.2" -fs 1 -xo 0.1 -fw 0.08 -fh 0.1 -fx -0.0 -line on "
   Main_tit((x00+x01)/2, y00+0.35, Model ": " Var.ivar " Member_" MEM)
   if (ivar=2); Main_tit((x00+x01)/2, y00+0.15, "(Hatching: p-value<0.5)"); endif
   "gxprint " Dat_out Score.ivar "_test_" Model "_mem_" MEM ".pdf pdf"

*   pull dd
   "c"
 Endwhile
 "close 1"
Endwhile
*******************************************************
function Main_tit(x, y, title)
  "set string 1 c 5.5 0"; "set strsiz 0.11 0.11"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return
*******************************************************
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 draw_ver_tit(x, y, title)
  "set string 1 c 5.5 90"; "set strsiz 0.09 0.10"
  "draw string "x" "y" "title
  "set string 1 c 5.5 0"
return
*******************************************************
function Line_pl(mark,colo,thick,styl)
"set gxout line"
*"set gxout stnmark"
"set cmark "mark
"set ccolor "colo
"set cthick "thick
"set cstyle "styl
return
********************************************
function Bar_pl(clor,gap)
"set gxout bar"
"set ccolor "clor
"set bargap "gap
"set barbase 0"
"set baropts outline"
"set cthick 2"
"set baropts filled"
*"set ccolor 1"
return
********************************************
   function Cor()
   levels = " -0.9 -0.8 -0.6 -0.4 -0.2 -0.1 -0.05 0 0.05 0.1 0.2 0.4 0.6 0.8 0.9 "
    "run color.gs -gxout grfill -kind blue->lightblue->white-(0)->white->orange->red -levs " levels
   "set grads off"
   return
********************************************
   function R_rmae()
   levels = " 10 20 30 50 70 100 150 200 250 300 500 "
   "run color.gs -gxout grfill -kind lightgreen->yellow->orange->red -levs " levels
   "set grads off"
   return
********************************************
   function R_Clor()
   levels = " 10 30 50 70 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 "
   "run color.gs -gxout grfill -kind lightgreen->yellow->orange->red -levs " levels
   "set grads off"
   return