[gradsusr] Vertical gradient of the saturation mixing ratio

[Morwakoma Matabane]

Good day,

I have managed to compute the saturation vapor pressure of two levels, say z-1 and z+1. And after computing the mixing ratios for the respective levels, I cannot seem to be able to calculate the dQs/dz (the vertical gradient of the mixing ratios). All I get whenever I attempt to compute the dQs/dz is an ERROR.
...

'set z 3'

*At temperature t1=t(z+1)

'define qs1=0.622*es1/lev'

*At temperature t2=t(z-1)

'define qs2=0.622*es2/lev'

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

'd dq'

OUTPUTS AN ERROR

How can I avoid the erroneous situation? Your help is highly appreciated.


Regards
Morwakoma
South African Weather Service (SAWS) Scientist


Morwakoma M Matabane BSc Witwatersrand, H. Dipl. Comp. Sc Witwatersrand, BSc Hon Pretoria, MSc Reading,  Research Scientist
South African Weather Service 442 Rigel Avenue South Erasmusrand Pretoria 0181
Tel: +27(0)12 367 6234 Fax:+27(0)12 367 6189
E-mail: morwakoma.matabane at weathersa.co.za Website:www.weathersa.co.za
[cid:image001.png at 01CE044D.69CE1730]
"Power our future with weather, climate and water"
"Big whorls have little whorls that feed on their velocity, and little whorls have lesser whorls and so on to viscosity", Lewis F Richardson.







-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130206/963bc8e0/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 4963 bytes
Desc: image001.png
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130206/963bc8e0/attachment-0003.png