[gradsusr] Mixing ratio gradient

[M.Mahakur]

Hi,

I feel t is in degree K and no need to add 273, if you are doing for conversion. secondly, please check your formulation ...

----- Original Message -----
From: Morwakoma Matabane <Morwakoma.Matabane at weathersa.co.za>
To: 'gradsusr at gradsusr.org' <gradsusr at gradsusr.org>
Sent: Wed, 06 Feb 2013 17:33:12 +0530 (IST)
Subject: [gradsusr] Mixing ratio gradient

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 --------------
A non-text attachment was scrubbed...
Name: Matabane.jpeg
Type: image/jpeg
Size: 7806 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130206/91a3c05f/attachment-0003.jpeg