*********************************************
****** note: adjust vrange as needed
******
*********************************************
*
'reset'
*
'open vrmse1_hdvar.des'
'open vrmse2_hdvar.des'
'open vrmse3_hdvar.des'
'open vbias1_hdvar.des'
'open vbias2_hdvar.des'
'open vbias3_hdvar.des'
*
'set display color white'
'clear'
'set mproj off'
********************************
********************************
********************************
********************************
*
lcount=1
while(lcount<2)
   if (lcount=1)
      x1=1.5
      x2=10.5
      y1=5.0
      y2=8.0
   endif
   if (lcount=2)
      x1=1.5
      x2=10.5 
      y1=1.0
      y2=4.0
   endif

   'set vpage 0.0 11.0 0.0 8.5'
   'set parea 'x1' 'x2' 'y1' 'y2 

****************specific humidity******************
  'set grid on 3 1'
  'set xlopts 1 6 0.16'
  'set ylopts 1 6 0.16'

  'set x 1 193'
  'set xlabs 00Z|12Z|00Z|12Z|00Z|12Z|00Z|12Z|00Z|12Z|00Z|12Z|00Z|12Z|00Z|12Z|'


  if (lcount=1) 
     'set y 3'  
     'set vrange  0.0, 1.5'
     'set ylab %.1f'
     'set yaxis  0.0  1.5 0.5'

     'set grads off'
     'set grid on'

     'set ccolor 9'
     'set cmark 0'
     'set cthick 6'
     'set cstyle 1'
     'set digsiz 0.06'
     'display rmse.1'
     'set string 9 tl 6'
     'set strsiz 0.15 0.15'
     'draw string 'x1+0.3' 'y2-0.15' ' q_bckgrmseavg ' g kg`a-1`n'

     'set ccolor 1'
     'set cmark 0'
     'set cthick 6'
     'set cstyle 1'
     'set digsiz 0.06'
     'display rmse.2'
     'set string 1 tl 6'
     'set strsiz 0.15 0.15'
     'draw string 'x1+2.7' 'y2-0.15' ' q_anlrmseavg ' g kg`a-1`n'

     'set ccolor 2'
     'set cmark 0'
     'set cthick 6'
     'set cstyle 1'
     'set digsiz 0.06'
     'display rmse.3'
     'set string 2 tl 6'
     'set strsiz 0.15 0.15'
     'draw string 'x1+5.1' 'y2-0.15' ' q_cressrmseavg ' g kg`a-1`n'

     'set strsiz 0.12 0.12'
     'set string 4 tl 6'
     'draw string 'x1+7.0' 'y1+0.20' 8-DAY TIME SERIES'

     if (filetype = pos1)
        'draw title USING OBS ASSIMILATED / RMSE FOR SPFH'
     endif
     if (filetype = neg3)
        'draw title CROSS-VALIDATION RMSE FOR SPFH'
     endif

     'draw xlab Time of Day (from 'starttime')'
     'set strsiz 0.15 0.15'
     'draw ylab RMSE (g kg`a-1`n)' 
  endif

  if (lcount=2) 
     'set y 3'  
     'set vrange -1.5 1.5'
     'set ylab %.1f'
     'set yaxis -1.5 1.5 0.5'

     'set grads off'
*    'set grid horizontal'
     'set grid on'

     'set ccolor 9'
     'set cmark 0'
     'set cthick 6'
     'set cstyle 1'
     'set digsiz 0.06'
     'display bias.4'
     'set string 9 tl 6'
     'set strsiz 0.15 0.15'
     'draw string 'x1+0.3' 'y2-0.15' ' q_bckgbiasavg ' g kg`a-1`n'   

     'set ccolor 1'
     'set cmark 0'
     'set cthick 6'
     'set cstyle 1'
     'set digsiz 0.06'
     'display bias.5'
     'set string 1 tl 6'
     'set strsiz 0.15 0.15'
     'draw string 'x1+2.7' 'y2-0.15' ' q_anlbiasavg ' g kg`a-1`n'   

     'set ccolor 2'
     'set cmark 0'
     'set cthick 6'
     'set cstyle 1'
     'set digsiz 0.06'
     'display bias.6'
     'set string 2 tl 6'
     'set strsiz 0.15 0.15'
     'draw string 'x1+5.1' 'y2-0.15' ' q_cressbiasavg ' kg`a-1`n'


     'set strsiz 0.12 0.12'
     'set string 4 tl 6'
     'draw string 'x1+7.0' 'y1+0.20' 8-DAY TIME SERIES'

     'draw xlab Time of Day (from 'starttime')'
     'set strsiz 0.15 0.15'
     'draw ylab BIAS (g kg`a-1`n)' 
  endif

  lcount=lcount+1
 endwhile
*****************************************************

 x1=1.5
 x2=11.0
*y1=0.5
 y1=4.5
 y2=8.5
 'set parea 'x1' 'x2' 'y1' 'y2
 'set string 9 tl 4'
 'set strsiz 0.15 0.15'
 'draw string 'x1+0.3' 'y1-0.3' BCKG IN PURPLE'
 'set string 1 tl 4'
 'draw string 'x1+3.0' 'y1-0.3' GSI ANL IN BLACK'
 'set string 2 tl 4'
 'draw string 'x1+5.8' 'y1-0.3' CRESS ANL IN RED'

 'enable print gg'
 'print'
 'disable print'
 '!gxps  -i -c gg -o  bias_rmse_for_2mQ.ps'
 '!/bin/rm -rf gg'
 'printim gg.gif gif white'
 '!mv gg.gif bias_rmse_for_2mQ.gif'
 'printim gg.png png white'
 '!mv gg.png bias_rmse_for_2mQ.png'
 'quit'
************************************************************