[gradsusr] VIMT

[sushant puranik]

Hi,
vint command will solve your issue. Following is an example showing
computation of meridional moisture between 1000 to 600hPa level

set lev 1000
d vint(lev,q*v,600)

Regards
Sushant
On 25 Aug 2015 20:05, "Sahidul" <sahidul at gmail.com> wrote:

Dear Help,


I want to calculate VIMF (Vertical Integrated Moisture flux) using WRF
data. I am not sure how to do this. I am using the following script to do
...Here PSFC is the SFC PRESSURE (Pa) and QVAPOR is Water vapor mixing
ratio (kg kg-1)


Could you please help on this regard.


************************************************

'reinit'

'clear'

'open /path/wrf/inp_grds.ctl'



'set grads off'

'set grid off'

'set xlopts 1 6 0.16'

'set ylopts 1 6 0.16'

'set mpdset hires'

'set lat 15 35'

'set lon 60  85'



'set gxout shaded'

'define q3=ave(vint(psfc*0.01,QVAPOR,100),t=13,t=36,1)'

'd q3'

'run cbarn.gs 1 1 8.9 4.2'

'set gxout vector'

*'set arrscl 0.6 1000'

'define q1=ave(vint(psfc*0.01,QVAPOR*skip(u,4),100),t=13,t=36,1)'

'define q2=ave(vint(psfc*0.01,QVAPOR*skip(v,4),100),t=13,t=36,1)'

'd q1;q2'

'draw title Vertically Integrated Moisture Transport(kg/m/s) '

'print file.eps'

'!convert +antialias -density 150 -rotate +90 -trim file.eps VIMT.jpg'

'quit'

****************************************************************************


With regards,


Sahidul

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150825/d3321315/attachment.html