[gradsusr] Script for Streamfunction

[sushant puranik]

Hi
I have made a script for plotting stream function. The script is as

'reinit'
'set display color white'
'clear'
'sdfopen D:/ECMWF_Monthly_Mean_data/SP/SPjun.nc'
'sdfopen D:/ECMWF_Monthly_Mean_data/wind/windjun.nc'
count = 1
while (count <=5)
  'clear'
  'set grads off'
  'set t 'count
  'q time'
date=subwrd(result,3)
hour=substr(date,1,3)
day=substr(date,4,2)
month=substr(date,6,3)
year=substr(date,9,10)
date_str=hour' 'day'/'month'/'year
'clear'
'set display color white'
'set map 1 1 6'
'set line  1 1 1'
'set grads off'
'enable print temp'
'define one=1+0*lat'
'define um=vint(lev,q*u.2,300)'
'define vm=vint(lev,q*v.2,300)'
'define psi=fish_psi(um,vm)'
'define upsi=-muadv(one,psi)'
'define vpsi=-mvadv(one,psi)'
'set cint 2'
'd psi/1e7'
'set gxout vector'
'd skip(upsi,4,4);skip(vpsi,4,4)'
*'d upsi;vpsi'
**'d skip(um,4,4);skip(vm,4,4)'
'draw title 'year'.'month''
'draw string 4.5 0.3 Streamfunction for vertically integrated water vapour
flux Unit:(kgs`a-1`n) x 10`a7`n'
'print'
'disable print'
'gxyat -x 1024 -y 840 D:/Transport/Streamplot/new'year'.png'
count = count + 1
endwhile
'clear'

But the results are opposite. Means where it should show negative(positive)
values it shows positive (negative).


any suggestion.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130314/dffccf69/attachment-0003.html