[gradsusr] VIMT

Sahidul sahidul at gmail.com
Tue Aug 25 10:33:18 EDT 2015


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150825/ea4fab59/attachment.html 


More information about the gradsusr mailing list