[gradsusr] level dimension didn't vary in gs file
resmi ea
reshmiea at gmail.com
Tue Jan 31 03:46:49 EST 2012
Dear all users
this is my gs file; the ctl file is given below. i want to file u-wind
anomaly fro the grid 0-90N and 0-357.5
at 1000 t0 10 hPa level for a single day(time= 1feb2008)
The grd file has shows the output at a single level
The problem is output result didnot vary at different levels
'reinit'
'sdfopen /media/resmi-1/windows/ncep/uwnd/200h_uwnd/uwnd.2000_10.nc'
'set lat 0 90'
'set lon 0 357.5'
*'set lev 850'
'set time 1feb2008'
'set gxout fwrite'
'set fwrite stj850_26j.grd'
n=1
while(n<=17)
'set z 'n
'define tt=uwnd(time=1feb2000)'
'define tt1=uwnd(time=1feb2001)'
'define tt2=uwnd(time=1feb2002)'
'define tt3=uwnd(time=1feb2003)'
'define tt4=uwnd(time=1feb2004)'
'define tt5=uwnd(time=1feb2005)'
'define tt6=uwnd(time=1feb2006)'
'define tt7=uwnd(time=1feb2007)'
'define tt8=uwnd(time=1feb2008)'
'define tt9=uwnd(time=1feb2009)'
'define tt10=uwnd(time=1feb2010)'
'define daymean=(tt+tt1+tt2+tt3+tt4+tt5+tt6+tt7+tt8+tt9+tt10)/11'
*'d daymean'
'define dayanom=tt9-daymean'
*'d dayanom'
n=n+1
endwhile
'disable fwrite'
Control file
DSET /media/resmi-1/windows/ncep/uwnd/200h_uwnd/stj850_26j.grd
undef -999.99
xdef 144 linear 0 2.5
ydef 37 linear 0 2.5
zdef 17 levels 1000 925 850 700 600 500 400 300 200 250 150 100 70 50 30
20 10
tdef 1 linear 01feb2008 1dy
TITLE stj850_26j.grd 144x37 grd
vars 1
uwnd 0 -999 uwnd temperature
endvars
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120131/8950c7ae/attachment-0003.html
More information about the gradsusr
mailing list