*subsetter
*Jean Pierre Arabonis
*19 August 2005
'open ./nww3.all.grb.ctl'
'q file'
var1 = htsgwsfc
var2 = windsfc
t1 = 1
t2 = 61

'disable fwrite'
'set gxout fwrite'
'set fwrite autotest'

'set lat -35'
'set lon 18'


i = 1
while i <= t2
'set t 'i
say 'printing to grib file at time 'i
'd 'var1
'd 'var2
i = i + 1
endwhile

'disable fwrite'
'close 1'