[gradsusr] stream function in pressure level

[resmi ea]

Dear users
 i have calculated stream function using open grads from 1000-10 hPa level
using NCEP/NCAR reanalysis data ex.1dec2009
The grads script give the grd file containing three variable like psi,upsi
and vpsi.
Here i have attached the figure. i find mistake in the figure...the
programme doesnot undergone integration part, at each level it calculate
three variable and finally dispalys.
 i want to know that whether it is correct or not. how can i modify the
programme
'reinit'
'sdfopen G:\windows\ncep\uwnd\uwnd.2009.nc'
'sdfopen G:\windows\ncep\uwnd\vwnd.2009.nc'
'set lon 0 357.5'
'set lat -90 90'
prompt 'key-in file name for psi:'
pull fname
'set fwrite 'fname
tt=1
while (tt<=1)
'set t 'tt
zz=1
while (zz<=17)
'set z 'zz
'define psi= fish_psi(uwnd,vwnd.2)'
'define one=1+0*lat'
'define upsi=mvadv(one,psi)'
'define vpsi=-muadv(one,psi)'
'set gxout fwrite'
'd psi/1e7'
'd upsi'
'd vpsi'
zz=zz+1
endwhile
tt=tt+1
endwhile
'disable fwrite'



-- 
Resmi E A
Research Scholar
Cochin university of science and technology
Cochin-16
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120522/83e6b86f/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1dec2008psi.jpeg
Type: image/jpeg
Size: 15795 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120522/83e6b86f/attachment-0001.jpeg