<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><div id="yiv1841082659">Hi group memeber,<br><br>I write a gs file to create a data set which is time series changing with height.<br>My original data is NARR netCDF data sets (8 times per day)<br>The ddf file as below:<br>==========<br>DSET d:\gradsfile\ncdata\uwnd.200812.nc<br>DTYPE netcdf<br>TITLE NARR 4D field Lambert Conformal Conic projection<br>UNDEF 32766 missing_value<br>UNPACK scale_factor=0.01f add_offset=47.66f<br>PDEF 349 277 lcc 1 -145.5 1 1 50 50 -107 32463.41 32463.41<br>XDEF 615 linear 150 0.3333<br>YDEF 255 linear 2 0.3333<br>ZDEF 29 levels 1000, 975, 950, 925, 900, 875, 850, 825, 800, 775, 750, 725, 700,<br> 650, 600, 550, 500, 450, 400, 350, 300, 275, 250, 225, 200, 175, 150, 125, 100<br>tdef 248 linear 00:00z01dec2008 3hr<br>vars 3<br>uwnd 29 t,z,y,x grid_eastward_wind<br>lat=>glat 0 y,x Latitude<br>lon=>glon 0 y,x
Longitude<br>endvars<br>===============<br>vwnd200812.ddf is simillar with above.<br><br>My gs file as below:<br>=============<br>'reinit'<br>'open d:\gradsfile\ncdata\uwnd200812.ddf'<br>'open d:\gradsfile\ncdata\vwnd200812.ddf'<br>*'set grads off'<br>*'set grid off'<br>'set gxout fwrite'<br>'set fwrite d:\gradsfile\ncdata\wind-height09z28dec2008.grd'<br>'set dfile 1'<br>'set lon 283.3'<br>'set lat 43.66'<br>*'set time 09z28dec2008'<br>'set dfile 2'<br>'set lon 283.3'<br>'set lat 43.66'<br>*'set time 09z28dec2008'<br>h=1000<br>t0=218<br>while(t0<222)<br>'set t 't0''<br>while(h>850)<br>'set lev 'h''<br>'u=0.01*uwnd.1+47.66'<br>'v=0.01*vwnd.2+47.66'<br>'wind=mag(u,v)'<br>'d wind'<br>h=h-25<br>endwhile<br>t0=t0+1<br>endwhile<br>say 'ok'<br>===============<br>My ctl file is as below:<br>==========<br>DSET d:\gradsfile\ncdata\wind-height09z28dec2008.grd'<br>DTYPE grib<br>TITLE Wind-Height 09z12dec2008<br>UNDEF 32766 missing_value<br>XDEF 1 linear 1
1<br>YDEF 1 linear 1 1<br>ZDEF 7 levels 1000, 975, 950, 925, 900, 875, 850<br>tdef 5 linear 03z28dec2008 3hr<br>vars 1<br>wind 7 99 wind_height<br>endvars<br>=========<br>When I open the ctl file from Grads, it shows:<br> can not open station/index map file (null)<br>the data file was not open.<br><br>I don't know what's wrong with it. Is there any one can tell me the reason? Thanks.<br><br>Regards<br><br><br>Emily<br>
York University, Toronto<br></div></td></tr></table>