[gradsusr] How to Calculate Vertically integrated Moisture Convergence Using Grads

[kalyan chakravarthy]

Dear Users,

Can anyone help me how to calculate vertically integrated moisture
convergence and moisture convergence at a particular level.

Below is my script to calculate Vertically Integrated Moisture Convergence.
sdfopen uwind1.nc (1000-300mb)
sdfopen vwind1.nc  (1000-300mb)
sdfopen sphum.nc (1000-300mb)

set lat 0 38
set lon 60 100
set mpdset rupres
set map 1 1 10
set LEV 1
set t 1 8
define qu=vint(1000,uwnd,300)
define qv=vint(1000,vwnd.2,300)
DEFINE SPHUM=VINT(1000,SHUM.3,500)
define mconv=-hdivg(qu*SPHUM,qv*SPHUM)
set gxout shaded
set display color white
SET CLEVS -3E-07 -2E-07 -1E-07 0 1E-07 2E-07 3E-07 4E-07 5E-07 6E-07
SET CCOLS 9 14 4 11 5 13 3 10 7 12 8 2 6
d mconv
cbar

after running the script it is showing error that incompatible grids.
Can anyone help me.


I also tried to calculate moisture convergence at particular pressure
level. But I am not sure up to what extent it is correct. I have got this
script from grads user forum.I am using NCEP/NCAR reanalysis data.

sdfopen uwind1.nc
sdfopen vwind1.nc
sdfopen sphum.nc
set lat 0 38
set lon 60 100
set mpdset rupres
set map 1 1 10
set LEV 850
set t 1
define mconv=-hdivg(shum.3*uwnd.1,shum.3*vwnd.2)
set gxout shaded
set display color white
SET CLEVS -3E-07 -2E-07 -1E-07 0 1E-07 2E-07 3E-07 4E-07 5E-07 6E-07
SET CCOLS 9 14 4 11 5 13 3 10 7 12 8 2 6
d mconv
cbar
draw title moisture convergence 850MB 13062013
printim moistureconverg13062013.png





-- 

*Regards*
*Y Kalyan Chakravarthy*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140623/1b2485b1/attachment.html