<P>
&nbsp; <BR>
Dear Grads-user,<BR>
&nbsp; &nbsp;  I am trying to calculate vertical integrated moisture convergence from our MM5 model output.I have written a script for it.But it is not working.I am getting the picture but it is not the right one.<BR>
&nbsp;  <BR>
'open domain3.ctl'<BR>
'set map 1 1 9'<BR>
'set gxout shaded'<BR>
'set t 25'<BR>
'define qu=q*u'<BR>
'define qv=q*v'<BR>
'define qq = -hdivg(qu,qv)'<BR>
'd vint(prs(z=3)/100.,qq,200)'<BR>
&nbsp; &nbsp;  Please suggest me to calculate this correctly. I have pressure level MM5 model output with following variables.<BR>
<BR>
&nbsp; &nbsp; ps 0 99 pstar (Pa)<BR>
&nbsp; &nbsp; pslv 0 99 sea level prs (mb)<BR>
&nbsp; &nbsp; u 20 99 u wind (m/s)<BR>
&nbsp; &nbsp; v 20 99 v wind (m/s)<BR>
&nbsp; &nbsp; q 20 99 mixing ratio (kg/kg)<BR>
&nbsp; &nbsp; z 20 99 height AGL (m)<BR>
&nbsp; &nbsp; td 20 99 dewpoint temp (C)<BR>
&nbsp; &nbsp; rh 20 99 rel humidity (%)<BR>
&nbsp; &nbsp; prs 20 99 pressure (Pa)<BR>
&nbsp; &nbsp; <BR>
I am also in doubt which of the above pressure&nbsp; variable should use to calculate moisture convergence.<BR>
&nbsp; &nbsp; Thanking you.<BR>
<BR>
Regards<BR>
Satyaban<BR>

</P>
<br><br>
<a href="http://adworks.rediff.com/cgi-bin/AdWorks/sigclick.cgi/www.rediff.com/signature-home.htm/1507191490@Middle5?PARTNER=3"><IMG SRC="http://adworks.rediff.com/cgi-bin/AdWorks/sigimpress.cgi/www.rediff.com/signature-home.htm/1963059423@Middle5?OAS_query=null&PARTNER=3" BORDER=0 VSPACE=0 HSPACE=0></a>