[gradsusr] Mixing ratio gradient

[Morwakoma Matabane]

Attached is the entire script and the output. The output "CONSTANT FIELD. VALUE -0".


'reinit'
'sdfopen http://monsoondata.org:9090/dods/gfsens/gfsens.2013020200'
'set lat -60 0'
'set lon -30 60'

'set z 3'

'set gxout shaded'
'define t=t+273'

'define es1=611*exp((2.5e6/461)*((1/t(z+1))-(1/273)))'
'define qs1=0.622*es1/(lev-es1)'

'define es2=611*exp((2.5e6/461)*((1/t(z-1))-(1/273)))'
'define qs2=0.622*es2/(lev-es2)'


'define dq=(qs1-qs2)/(lev(z+1)-lev(z-1))'

'd dq'

'cbarn 1.0 0 5.5 0.3'
'printim AA.gif gif -x1400 -y1200 white'
'quit'



Help highly appreciated.



Regards

Morwakoma

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130206/254f3a9e/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: AA.gif
Type: image/gif
Size: 8991 bytes
Desc: AA.gif
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130206/254f3a9e/attachment-0003.gif