'reinit'
'set display color white'
'c'
*'open CST/JACKSON/SAIDAMODELIBIS.S001.ctl'
*'open CST/MILENA/SAIDAMODELIBIS.S001.ctl'
'open Nwp/JACKSON/SAIDAMODELIBIS.S001.ctl'
'open Nwp/MILENA/SAIDAMODELIBIS.S001.ctl'


'set z 1'
'set time 03z01feb2015'
'q dims'
lin=sublin(result,5)
tti=subwrd(lin,9)

'set time 03z28feb2015'
'q dims'
lin=sublin(result,5)
ttf=subwrd(lin,9)
xmin=1.0 
xmax=10.5 
ymin=1 
ymax=7.5
'set parea 'xmin' 'xmax' 'ymin' 'ymax''
'set font 5'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set grads off'
'set grid off'
'set xlabs 0|3|6|9|12|15|18|21|'
'set t 'tti' 'tti+23''
'set vrange -20 900'
'd tloop(ave(swrad.1 ,t+1,t='ttf',24)) '
'd tloop(ave(evap.1  ,t+1,t='ttf',24)) '
'd tloop(ave(evap.2  ,t+1,t='ttf',24)) '
'd tloop(ave(slhf.2  ,t+1,t='ttf',24)) '
'set strsiz 0.25 0.25'
'set string 1 c 0.25 0'
'draw string  'xmin+4.5' 'ymin-0.7' Local Time '
'draw string  'xmin+5.0' 'ymax+0.7' [Drie Season] Latent Heat Fluxo (W/m^2)  '

'q pos'
'c'
'set grads off'
'set grid off'
'set xlabs 0|3|6|9|12|15|18|21|'
'set t 'tti' 'tti+23''
'set vrange -20 900'
'd tloop(ave(swrad.1 ,t+1,t='ttf',24)) '
'd tloop(ave(sens.1  ,t+1,t='ttf',24)) '
'd tloop(ave(sens.2  ,t+1,t='ttf',24)) '
'd tloop(ave(sshf.2  ,t+1,t='ttf',24)) '
'set strsiz 0.25 0.25'
'set string 1 c 0.25 0'
'draw string  'xmin+4.5' 'ymin-0.7' Local Time '
'draw string  'xmin+5.0' 'ymax+0.7' [Drie Season] Sensible Heat Fluxo (W/m^2)  '
'q pos'
'c'
'set grads off'
'set grid off'
'set xlabs 0|3|6|9|12|15|18|21|'
'set t 'tti' 'tti+23''
'set vrange -20 500'
'd tloop(ave(swrad.2 ,t+1,t='ttf',24)) '
'd tloop(ave((swrad.1-swup.1) +(lwrad.1 -firb.1)-evap.1-sens.1 ,t+1,t='ttf',24)) '
*'d tloop(ave((swrad.2-swup.2) +(lwrad.2 -firb.2)-evap.2-sens.2 ,t+1,t='ttf',24)) '
'set strsiz 0.25 0.25'
'set string 1 c 0.25 0'
'draw string  'xmin+4.5' 'ymin-0.7' Local Time '
'draw string  'xmin+5.0' 'ymax+0.7' [Drie Season] Net Radiation (W/m^2)  '