How to use the fwrite function to get a specific time range and variable?

[Henri Pinheiro]

Dear grads users,

I have a binary data set and would want an output from a specific time range
and variable.
The input has 1459 time steps and 16 variables and I used a fwrite script to
get a smaller period.

The error message was
SET FWrite Error:  fwrite file is open
Use DISABLE FWRITE command to close file

Please, does anybody have any suggestion?
the fwrite.gs and file.ctl files are beneath

Best Regards
Henri

## fwrite.gs

 'reinit'
'open diaro.ctl'
'set gxout fwrite'
'set y 1 96'
'set x 1 192'

t=237
while (t<=604)

'set t 't

'set fwrite u300_MAM1989.dat'
'd u300'
t=t+1
endwhile

'disable fwrite'
'set gxout shaded'
'set gxout contour'

'quit'



## file.ctl

dset /rede/cycledi/diarioAGCM50/17/17/19500117_1989.bin
title file created by grads using lats4d
undef 1e+20
options big_endian template
xdef 192 linear 0.000000 1.875000
ydef 96 levels
-88.572 -86.723
-84.862 -82.999 -81.135 -79.271 -77.406 -75.541 -73.676 -71.811 -69.946
-68.081
-66.216 -64.351 -62.486 -60.620 -58.755 -56.890 -55.025 -53.160 -51.294
-49.429
-47.564 -45.699 -43.833 -41.968 -40.103 -38.238 -36.372 -34.507 -32.642
-30.777
-28.911 -27.046 -25.181 -23.316 -21.450 -19.585 -17.720 -15.855 -13.989
-12.124
-10.259  -8.394  -6.528  -4.663  -2.798  -0.933   0.933   2.798   4.663
6.528
  8.394  10.259  12.124  13.989  15.855  17.720  19.585  21.450  23.316
25.181
 27.046  28.911  30.777  32.642  34.507  36.372  38.238  40.103  41.968
43.833
 45.699  47.564  49.429  51.294  53.160  55.025  56.890  58.755  60.620
62.486
 64.351  66.216  68.081  69.946  71.811  73.676  75.541  77.406  79.271
81.135
 82.999  84.862  86.723
88.572
zdef 1 linear 1 1
tdef 1459 linear 00Z01jan1989  6hr
vars 16
ghmt500      0  99   TIME MEAN GEOPOTENTIAL HEIGHT500hpa (M ) [] []
ghmt200      0  99   TIME MEAN GEOPOTENTIAL HEIGHT200hpa (M ) [] []
spmt         0  99   TIME MEAN SEA LEVEL PRESSURE (Mb ) []
u850         0  99   TIME MEAN ZONAL WIND (U850hpa) (M/Sec ) []
u300         0  99   TIME MEAN ZONAL WIND (U300hpa) (M/Sec ) []
u200         0  99   TIME MEAN ZONAL WIND (U200hpa) (M/Sec ) []
v850         0  99   TIME MEAN MERIDIONAL WIND (U850hpa) (M/Sec ) []
v300         0  99   TIME MEAN MERIDIONAL WIND (U300hpa) (M/Sec ) []
v200         0  99   TIME MEAN MERIDIONAL WIND (U200hpa) (M/Sec ) []
fcmt200      0  99   TIME MEAN STREAM FUNCTION (M**2/Sec ) []
pvmt200      0  99   TIME MEAN VELOCITY POTENTIAL (M**2/Sec ) []
prec         0  99   TOTAL PRECIPITATION (Kg M**-2 Day**-1) []
rhmt850      0  99   TIME MEAN RELATIVE HUMIDITY850hpa (No Dim ) []
rhmt700      0  99   TIME MEAN RELATIVE HUMIDITY700hpa (No Dim ) []
swos         0  99   SOIL WETNESS OF SURFACE (No Dim ) []
role         0  99   OUTGOING LONG WAVE AT TOP (W M**-2 ) []
endvars



-- 
Henri Rossi Pinheiro
Meteorology Master Student
National Institute for Space Research (INPE)
Weather Forecast and Climate Studies Center (CPTEC)
Cachoeira Paulista, SP/Brazil, 12630-000 Tel.:+55-12-3186-8607
e-mails: henri.pinheiro at cptec.inpe.br
              henri.inpe at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100408/a2b4cc22/attachment.html