**** NOTE
**** use portrait ******


var=prec000120
input_title='HAFS Ensemble\mean TOTAL PRECIPITATION (in) DURING DAY 1\Storm: DELT(26L) init 2020100700'

var='prec024048'
input_title='HAFS Ensemble\MEAN TOTAL PRECIPITATION (in) DURING 24-48 h\Storm: NINE(09L) init 2024092400'
outfile='prec024048.png'

say 'draw 'var

'open ens_prob.ctl'
'c'
'set parea 1.0 7.5 1.4 7.9'
*'set lon 'lon1' 'lon2
*'set lat 'lat1' 'lat2 
'set grads off'
'set mpdset mres'
'run rgbset.gs'
'set mpt 0 79 3 3'
'set mpt 1 79 3 3'
'set mpt 2 79 3 3'
'set stat on'
'set gxout shaded'
'set xlopts 1 5 0.18'
'set ylopts 1 5 0.18'
'run newcolor.gs'
'set clevs  1    2   4  8  16  24  32 '
'set ccols 0  31  32  33 34  35  36  37 '

* draw a faked data, in case all are badval
'd rainbow'


'draw xlab LON(degree)'
'draw ylab LAT(degree)'
**'draw title 'input_title1'\'input_title2 ' Max='smax
'draw title 'input_title

** find x, y range on current displayed plot
'query gxinfo'
rec3 = sublin(result,3)
rec4 = sublin(result,4)
xlim1 = subwrd(rec3,4)
xlim2 = subwrd(rec3,6)
ylim1 = subwrd(rec4,4)
ylim2 = subwrd(rec4,6)
say 'xlim1 = 'xlim1'  ylim1 = 'ylim1
say 'xlim2 = 'xlim2'  ylim2 = 'ylim2

** draw color bar
xbar=(xlim1+xlim2)/2.0
ybar=ylim1-1.0
'run cbarn 0.9 0 'xbar' 'ybar

* need to set clevs ccols again
'set clevs  1    2   4  8  16  24  32 '
'set ccols 0  31  32  33 34  35  36  37 '
*'d 'var
'd smth9('var')'
stat=result
line8=sublin(stat,8)
smin=subwrd(line8,4)
smax=subwrd(line8,5)

** draw max value
xs=xlim1+0.5
ys=ylim1-0.6
smaxformat=math_format('%-7.3f',smax)
'draw string 'xs' 'ys'  Max='smaxformat'(in)'



* draw mean track if available 
** member id#
*'set line 15 1 1'
*'run oplot_xy data/track_mean' 0'


* draw tracks of members

** member id#
start_mem=0
end_mem=30
m=start_mem
while(m <= end_mem)
say 'm='m
smem=m
 if(m<10) ; smem='0'm ; endif
** if(t >=10 & t<100) ; ts='0't ; endif

*'set line 15 1 1'
'set line 1 1 1'
*'run oplot_xy data/track_mem'smem' 3'
'run oplot_xy plot_track_mem'smem' 0'

m=m+1
endwhile

'printim 'outfile' x900y900 white'


***************************************