[gradsusr] Problem in calculating moisture convergence at 200 hpa level.

[Arlindo Meque]

Dear all.


 
I'm trying to compute moisture
convergence  at 200 hpa. I've successfully done the calculations for
850 and 500 hpa.  But for 200 hpa I am getting a very strange output,
please see the attachment. 

 

This is the small script I am using for
that purpose:
 
'reinit' 
'reset' 
'sdfopen era.uq200.nc ' 
'sdfopen era.vq200.nc' 
'set t 1 last' 
'define diveuq= -(fish_div(u.1,v.2))' 
'define mconv=diveuq' 
'set sdfwrite ERA.mconv.200.nc' 
'sdfwrite mconv'
"era.uq200.nc" is a  U wind component  at
200 hpa level multiplied by  the specific humidiry athe the same
level.


Any idea?


Arlindo 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130201/07a92a1f/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.mconv.200.png
Type: image/png
Size: 108934 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20130201/07a92a1f/attachment-0003.png