<div dir="ltr"><div class="gmail_default" style="font-family:arial narrow,sans-serif">Hello Jeff,<br clear="all"></div><div class="gmail_default" style="font-family:arial narrow,sans-serif"><br></div><div class="gmail_default" style="font-family:arial narrow,sans-serif">I have used the vint() function from Grads. The pressure is in millibars, q (sp. humidity) is kg/kg, u is in m/s, v is in m/s.</div><div class="gmail_default" style="font-family:arial narrow,sans-serif"><br></div><div class="gmail_default" style="font-family:arial narrow,sans-serif">Please find the attached gs files IVT.gs and <a href="http://vintmoisture.gs">vintmoisture.gs</a></div><div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr"><font color="#0000ff"><b><i><span class="gmail_default" style="font-family:"arial narrow",sans-serif"><br></span></i></b></font></div><div><font color="#0000ff"><b><i><span class="gmail_default" style="font-family:"arial narrow",sans-serif">I have used these functions </span></i></b></font></div><div><font color="#0000ff"><b><i><span class="gmail_default" style="font-family:"arial narrow",sans-serif"><br></span></i></b></font></div><div>'define um=q*u'<br>'define vm=q*v'<br>'define ium=vint(1000,um,<span class="gmail_default" style="font-family:"arial narrow",sans-serif">30</span>0)'<br>'define ivm=vint(1000,vm,<span class="gmail_default" style="font-family:"arial narrow",sans-serif">300</span>)'<br>'set gxout shaded'<br>'d mag(1.e-2*ium,1.e-2*ivm)'<span class="gmail_default" style="font-family:"arial narrow",sans-serif"> (from <a href="http://moisture.gs">moisture.gs</a> file)</span><br>'vimf=-hdivg(ium,ivm)<br></div><div><br></div><div><div class="gmail_default" style="font-family:"arial narrow",sans-serif">Thanks for your help.</div><div class="gmail_default" style="font-family:"arial narrow",sans-serif"><br></div><div class="gmail_default" style="font-family:"arial narrow",sans-serif">Regards</div><div class="gmail_default" style="font-family:"arial narrow",sans-serif"><br></div></div><div dir="ltr"><font color="#0000ff"><b><i>Kalyan Chakravarthy Yesoda</i></b></font><div><br></div><div>Graduate Student, PhD</div><div>University of Nevada, Reno</div><div>Mobile no: 775-420-9434</div><div><br></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Nov 13, 2019 at 8:16 PM Jeff Duda <<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Are you sure your input fields have the correct (SI) units? Also, i don't know what grads formula you're using to perform the integral, but vint() assumes pressure coordinates in <u>millibars</u> rather than Pascals, so that may be a source of your problem.<div><br></div><div>I don't work with IVT, and the formulas I found online don't match what you indicated, so I can't say for sure if the problem is with the formulation or not.</div><div><br></div><div>Jeff Duda</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Nov 13, 2019 at 8:39 PM kalyan chakravarthy <<a href="mailto:kalyan4space@gmail.com" target="_blank">kalyan4space@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif">Hello Everyone,<br clear="all"></div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif">I need some help to calculate Integrated Water vapor transport (IVT) from surface to 300hPa.  I am using ERA-Interim data (uwind (u), vwind(v), specific humidity(q)). I have used the scripts that are published in the grads user forum, but the results seems wrong (ex. the values are from 40000 kg/ms). </div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif">The formula I am using to calculate the IVT is</div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif">IVT=Sqrt((integral(q*u)*dp)^2+((integral(q*v)*dp)^2)</div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif"><br></div><div class="gmail_default" style="color:rgb(34,34,34);font-size:14px;font-family:"arial narrow",sans-serif">Integration is from Surface to 300 hPa.</div><div><br></div><div><br></div><div>Thanks</div><div><br></div><div>Kalyan</div><br><br>-- <br><div dir="ltr"><font color="#0000ff"><b><i>Kalyan Chakravarthy Yesoda</i></b></font><div><br></div><div>Graduate Student, PhD</div><div>University of Nevada, Reno</div><div>Mobile no: 775-420-9434</div><div><br></div></div><br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font size="2"><span>Jeff Duda, Research Scientist</span></font></div><div dir="ltr"><font size="2"><span></span></font>University of Colorado Boulder</div><div dir="ltr"><font size="2"><span></span></font>Cooperative Institute for Research in Environmental Sciences</div><div dir="ltr">NOAA/OAR/ESRL/Global Systems Division<br><font size="2"><span>

<span>Boulder, CO<br></span></span></font>



</div></div></div></div></div></div></div></div>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">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>
</blockquote></div>