*
* credit for technique: Vaughan Phillips
*
'reinit'
'open ifamean.ctlt'
'set x 1'
'set y 1'
'set z 1'
'set t 393 440'
*
* note: prec in units of (mm/h), dt = 3 hours, paccum in units of (mm)
*
dt=3
'paccum = 'dt'*tloop(sum(prec,t=393,t+0))'
'enable print prec_accum.gx'
'set vrange 0. 200.'
'd paccum'
'draw title TOGA Accumulated Precipitation'
'draw ylab (mm)'
'print'
'disable print'
'!gxps -c -d -i prec_accum.gx -o prec_accum.ps'
'!gzip -f prec_accum.ps'
'!rm prec_accum.gx'
