<p dir="ltr">Hi,<br>
vint command will solve your issue. Following is an example showing computation of meridional moisture between 1000 to 600hPa level</p>
<p dir="ltr">set lev 1000<br>
d vint(lev,q*v,600)</p>
<p dir="ltr">Regards<br>
Sushant</p>
<div class="gmail_quote">On 25 Aug 2015 20:05, "Sahidul" <<a href="mailto:sahidul@gmail.com">sahidul@gmail.com</a>> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">Dear Help,</p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">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)<br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">Could you please help on this regard.<br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">************************************************<br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'reinit'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'clear'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'open /path/wrf/inp_grds.ctl'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"> </p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set grads off'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set grid off'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set xlopts 1 6 0.16'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set ylopts 1 6 0.16'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set mpdset hires'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set lat 15 35'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set lon 60<span> </span>85'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"> </p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set gxout shaded'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'define q3=ave(vint(psfc*0.01,QVAPOR,100),t=13,t=36,1)'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'d q3'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'run <a href="http://cbarn.gs" target="_blank">cbarn.gs</a> 1 1 8.9 4.2'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'set gxout vector'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">*'set arrscl 0.6 1000'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'define q1=ave(vint(psfc*0.01,QVAPOR*skip(u,4),100),t=13,t=36,1)'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'define q2=ave(vint(psfc*0.01,QVAPOR*skip(v,4),100),t=13,t=36,1)'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'d q1;q2'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'draw title Vertically Integrated Moisture Transport(kg/m/s) '</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'print file.eps'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'!convert +antialias -density 150 -rotate +90 -trim file.eps VIMT.jpg'</p>
<p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">'quit'</p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">****************************************************************************</p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">With regards,</p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal"><br></p><p class="MsoNormal" style="margin-bottom:0.0001pt;line-height:normal">Sahidul<br></p>
</div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div>