<div dir="ltr">Check example #2 in the documentation for the vint function:<div><br></div><div><p style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium">2. This is an artificial example that demonstrates a vertical integration from a fixed lower bound of 1000mb to the top of the atmosphere, and integrating a field of all 1's. This gives an answer of 10204.1 (or 100000/9.8) which is the mass of air (in kg) of a 1 meter squared column when the surface pressure is 1 bar and the accelleration due to gravity is assumed to be exactly 9.8m/sec**2 over the entire column.</p>
<ul style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium"><code>vint(const(ps,1000),const(t,1),0)</code><code><div><br></div></code></ul><font color="#000000" face="arial, helvetica, sans-serif">Grads is probably having a hard time figuring out what you mean by 'lev' in the first argument. My guess is the same numerical values of lev are being used in each of the three commands. The value is probably 1000, so in effect you are integrating over progressively deeper layers, which is why you are getting fqu1 < fqu2 < fqu3 for example. I would try changing the first argument to something like const(p,1000), const(p,700), const(p,500) etc.</font></div>
<div><font color="#000000" face="arial, helvetica, sans-serif"><br></font></div><div><font color="#000000" face="arial, helvetica, sans-serif">Jeff Duda<br></font><ul style="color:rgb(0,0,0);font-family:'Times New Roman';font-size:medium">
<code><br></code></ul></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Jan 23, 2014 at 1:17 PM, °° Maytę °° . <span dir="ltr"><<a href="mailto:mayteleal@hotmail.com" target="_blank">mayteleal@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div dir="ltr"><div><div dir="ltr"><br><div>
<div dir="ltr">Hi,<br><br><span lang="en"><div class="im"><span><span lang="en"><span>I'm trying to</span> <span>calculate the</span> <span>transport</span> <span>and convergence</span> <span>of</span> <span>moisture</span> <span>flux</span> <span>integrated</span> <span>in three</span> <span>layers</span> <span>of the troposphere</span><span>. </span><span>Define</span> <span>the layers</span> <span>as follows:</span> </span></span><span></span><br>
<br> <span>* Lower troposphere</span></div><span>:<br>'set lev 1000'<br></span> <span>'defines</span> <span>fqu1</span> <span>vint</span> <span>=</span> <span>(lev, q</span><span>*u</span><span></span><span>,</span> <span>850)</span><span>'</span><span> </span><span></span><br>
<div class="im"><br> <span>* Middle troposphere</span><span>:</span><br> <span>'set</span> <span>lev</span> <span>700</span><span>'</span><br> <span>'</span><span>defines</span> <span>fqu2</span> <span>vint</span> <span>=</span> <span>(</span><span>lev,</span> <span>u*</span><span>q</span><span></span><span>,</span> <span>500</span><span>)</span><span>'</span><br>
<br> <span>* High</span> <span>Troposphere</span><span>:</span><br> <span>'set</span> <span>lev</span> <span>500</span><span>'</span><br> <span>'</span><span>defines</span> <span>fqu3</span> <span>vint</span> <span>=</span> <span>(</span><span>lev,</span> <span>u</span>*<span>q</span><span></span><span>,</span> <span>300</span><span>)</span><span>'</span><br>
</div></span><div class="im"><br><span lang="en"><span lang="en"><span>The same for</span> <span>fqv1, fqv2 and fqv3...<br><br></span></span><span>For then</span><span>, the values at</span> <span>low levels</span> <span>are</span> <span>lower than</span> <span>in the</span> <span>middle and upper</span> <span>troposphere.</span><br>
<span></span></span><span lang="en"><span><span lang="en"><span>Is there another</span> <span>way to calculate</span><span>?</span> </span></span><span></span><br><br> <span>Something is</span> <span>wrong, since</span> <span>I reviewed</span> <span>all the steps</span> <span>and the only</span> <span>change</span> <span>I make</span> <span>is</span> <span>the</span> <span>'</span><span>lev</span><span>'</span><span>.</span><br>
<span>Before</span> <span>had done</span> <span>the integration</span> <span>from 1000 to</span> <span>300 hPa</span> <span>and not</span> <span>had any problem</span><span>.</span><br><br> <span>thank you</span><span>,</span><br>
<span>May</span></span><br> </div></div></div> </div></div> </div></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" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div>