<div dir="ltr"><div><div><div>Dear users <br></div>One thing being I have been wondering about vertically Intergarted moisture flux is that when I use this  formular suggested by Jenniffer here, the moisture will be constant even if I vary the time but the winds can be seen as changing.<br>Can someone help me. I want  to plot moisure flux at different times i,e <br><ol><li>  start of the season</li><li>peak of the season</li><li>cessation of the season etc . That means time must vary<br></li></ol>&#39;p925 = const(shum.1,925,-a)&#39;<br>&#39;define qu=vint(p925,shum.1*uwnd.2,850)&#39;<br>&#39;define qv=vint(p925,shum.1*vwnd.3,850)&#39;<br>&#39;define quave=ave(qu,t=1,t=2)&#39;<br>&#39;define qvave=ave(qv,t=1,t=2)&#39;<br><br><br></div>see my script and figures plotted which shows that even if I vary time, the moisture is the same <br><br></div>thanks<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jun 2, 2017 at 10:27 AM, afwande juliet <span dir="ltr">&lt;<a href="mailto:afwandej965@gmail.com" target="_blank">afwandej965@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>thanks<br></div>This seems to generate no error but am not sure what it does since I wanted to plot moisture flux vertically intergrated at 850mb and 500mb and 200mb. My shum data varies from 1000-300mb , I am not sure why here its defined as  constant, -a(negative acceleration ?)<br></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jun 2, 2017 at 6:14 AM, Sridhara Nayak <span dir="ltr">&lt;<a href="mailto:sridharanayakiitkgp@gmail.com" target="_blank">sridharanayakiitkgp@gmail.com</a><wbr>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>How about the following?<br><br>sdfopen <a href="http://shum.nc" target="_blank">shum.nc</a><br>sdfopen <a href="http://uwnd.nc" target="_blank">uwnd.nc</a><br>sdfopen <a href="http://vwnd.nc" target="_blank">vwnd.nc</a><br><br>p1000 = const(shum.1,1000,-a)<br>define qu=vint(p1000,shum.1*uwnd.2,10<wbr>0)<br>define qv=vint(p1000,shum.1*vwnd.3,10<wbr>0)<br>d -hdivg(qu,qv)<br><br></div>Best Regards<span class="m_-5435658607204180593HOEnZb"><font color="#888888"><br></font></span></div><span class="m_-5435658607204180593HOEnZb"><font color="#888888">Sridhara<br><div><div><br></div></div></font></span></div><div class="m_-5435658607204180593HOEnZb"><div class="m_-5435658607204180593h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jun 2, 2017 at 10:01 AM, afwande juliet <span dir="ltr">&lt;<a href="mailto:afwandej965@gmail.com" target="_blank">afwandej965@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Dear users of grads<br></div>did someone manage to plot vertically integrated moisture flux.<br></div>I have tried the function vint but gets the below error<br><a href="http://cola.gmu.edu/grads/gadoc/gradfuncvint.html" target="_blank">http://cola.gmu.edu/grads/gado<wbr>c/gradfuncvint.html</a> <br>ga-&gt; q file<br>File 1 : mean daily NMC reanalysis (1997)<br>  Descriptor: <a href="http://shum.1997.nc" target="_blank">shum.1997.nc</a><br>  Binary: <a href="http://shum.1997.nc" target="_blank">shum.1997.nc</a><br>  Type = Gridded<br>  Xsize = 144  Ysize = 73  Zsize = 8  Tsize = 365  Esize = 1<br>  Number of Variables = 1<br>     shum  8  t,z,y,x  mean Daily specific humidity<br>ga-&gt; set t 1 <br>Time values set: 1997:1:1:0 1997:1:1:0 <br>ga-&gt; set z 1 <br>LEV set to 1000 1000 <br>ga-&gt; set z 2 <br>LEV set to 925 925 <br>ga-&gt; set z 3 <br>LEV set to 850 850 <br>ga-&gt; set z 8 <br>LEV set to 300 300 <br>ga-&gt; d vint(850,shum.1*lterp(uwnd.2,s<wbr>hum.1),600)<span><br>Error from VINT:  Incompatible grids. <br></span>Contouring: -0.25 to 0.2 interval 0.05 <br>*** Error in `grads&#39;: free(): invalid pointer: 0x000000000145a430 ***<br>======= Backtrace: =========<br>/lib/x86_64-linux-gnu/libc.so.<wbr>6(+0x777e5)[0x7f41ce33c7e5]<br>/lib/x86_64-linux-gnu/libc.so.<wbr>6(+0x7fe0a)[0x7f41ce344e0a]<br>/lib/x86_64-linux-gnu/libc.so.<wbr>6(cfree+0x4c)[0x7f41ce34898c]<br>grads(gagfre+0x5a)[0x449e5a]<br>grads(gagrel+0x58)[0x493be8]<br>grads(gadspl+0x668)[0x4ab318]<br>grads(gacmd+0x17ae)[0x4aa33e]<br>grads(main+0x62e)[0x41277e]<br>/lib/x86_64-linux-gnu/libc.so.<wbr>6(__libc_start_main+0xf0)[0x7f<wbr>41ce2e5830]<br>grads(_start+0x29)[0x412e79]<br>======= Memory map: ========<br>00400000-004e8000 r-xp 00000000 08:05 524842                        <wbr>     /usr/bin/grads<br>006e7000-006e8000 r--p 000e7000 08:05 524842                        <wbr>     /usr/bin/grads<br>006e8000-006ec000 rw-p 000e8000 08:05 524842                        <wbr>     /usr/bin/grads<br>006ec000-00992000 rw-p 00000000 00:00 0 <br>012a4000-014b8000 rw-p 00000000 00:00 0   <br></div><div class="m_-5435658607204180593m_-530491365356031727HOEnZb"><div class="m_-5435658607204180593m_-530491365356031727h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 31, 2017 at 1:00 AM, Jennifer M Adams <span dir="ltr">&lt;<a href="mailto:jadams21@gmu.edu" target="_blank">jadams21@gmu.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Usage Note #5: The vint function operates only in an X-Y varying dimension environment.<br>
<br>
So, you will have to fix A, and if you must have T varying, it will probably work with ‘define’ instead of ‘display’, but I would try fixing T first just to make sure you’re getting the result you intend:<br>
<br>
set z 1<br>
set t 182<br>
d vint(etc. etc.)<br>
<br>
If that works and you are sure your expression is what you want, then you can<br>
<br>
set t 182 274<br>
define yourvname = vint(etc. etc.)<br>
<br>
… and GrADS will loop over T under the hood to create your defined variable that varies X, Y, and T.<br>
<span class="m_-5435658607204180593m_-530491365356031727m_-3965267870390237335HOEnZb"><font color="#888888">—Jennifer<br>
</font></span><div class="m_-5435658607204180593m_-530491365356031727m_-3965267870390237335HOEnZb"><div class="m_-5435658607204180593m_-530491365356031727m_-3965267870390237335h5"><br>
<br>
On May 30, 2017, at 5:43 PM, Pablo Camarena &lt;<a href="mailto:jcamarena006@gmail.com" target="_blank">jcamarena006@gmail.com</a>&gt; wrote:<br>
<br>
&gt;<br>
&gt; Thanks Jennifer,<br>
&gt;<br>
&gt; Now I&#39;m trying the following:<br>
&gt;<br>
&gt; sdfopen <a href="http://shum.mon.mean.nc" rel="noreferrer" target="_blank">shum.mon.mean.nc</a><br>
&gt; sdfopen <a href="http://uwind.mon.mean.nc" rel="noreferrer" target="_blank">uwind.mon.mean.nc</a><br>
&gt; sdfopen <a href="http://vwind.mon.mean.nc" rel="noreferrer" target="_blank">vwind.mon.mean.nc</a><br>
&gt; set lon 225 227<br>
&gt; set lat 23 24<br>
&gt; set t 182 274<br>
&gt; set z 1 3<br>
&gt; d vint(1000,shum.1*lterp(uwnd.2,<wbr>shum.1),850)<br>
&gt; And I get the following error:<br>
&gt; &quot;Operation error:  Invalid dimension environment<br>
&gt;   Too many varying dimensions&quot;<br>
&gt;<br>
&gt; I think the error is with vertical levels.<br>
&gt;<br>
&gt; 2017-05-30 14:08 GMT-06:00 Jennifer M Adams &lt;<a href="mailto:jadams21@gmu.edu" target="_blank">jadams21@gmu.edu</a>&gt;:<br>
&gt; Hi, Pablo —<br>
&gt; Grads shouldn’t be dumping core on you, but at least you got a clue before it crashed. The grids of your shum.1 and uwnd.2 files apparently don’t match. Try using the lterp() function to interpolate in the lat/lon dimenions, something like this:<br>
&gt;<br>
&gt; ga-&gt; d vint(1000,shum.1*lterp(uwnd.2,<wbr>shum.1),600)<br>
&gt;<br>
&gt; In the script lines you posted, I noticed that the syntax here is incorrect:<br>
&gt; &gt; qu=vint(1000,shu,.1*uwnd.2,300<wbr>)<br>
&gt; &gt; qv=vint(1000,shu,.1*vwnd.3,300<wbr>)<br>
&gt;<br>
&gt; change &quot;shu,.1&quot; to &quot;shum.1&quot;<br>
&gt;<br>
&gt; If the problem is that the vertical levels are different in your data files, then you need to fix that. There’s a script in the library that may be helpful —<a href="ftp://cola.gmu.edu/grads/scripts/zinterp.gs" rel="noreferrer" target="_blank">ftp://cola.gmu.edu/grads/scri<wbr>pts/zinterp.gs</a>.  I believe the vint() function is a bit fussy about the Z axis. The lterp() function does not do vertical interpolation. Please read the vint() function Usage Notes carefully (<a href="http://cola.gmu.edu/grads/gadoc/gradfuncvint.html" rel="noreferrer" target="_blank">http://cola.gmu.edu/grads/gad<wbr>oc/gradfuncvint.html</a>) to make sure you’ve got everything right. If you’re still having trouble, please also include the output from &quot;ncdump -c” on your data files, the version of GrADS you are using, and the shortest possible script that illustrates the error.<br>
&gt; —Jennifer<br>
&gt;<br>
&gt;<br>
&gt; On May 30, 2017, at 2:18 PM, Pablo Camarena &lt;<a href="mailto:jcamarena006@gmail.com" target="_blank">jcamarena006@gmail.com</a>&gt; wrote:<br>
&gt;<br>
&gt; &gt; Hi,<br>
&gt; &gt;<br>
&gt; &gt; I want to calculate Vertically integarted moisture flux convergence (VIMFC)<br>
&gt; &gt;<br>
&gt; &gt; But when I give the following instruction it marks me errors:<br>
&gt; &gt;<br>
&gt; &gt; sdfopen <a href="http://shum.mon.mean.nc" rel="noreferrer" target="_blank">shum.mon.mean.nc</a><br>
&gt; &gt; sdfopen <a href="http://uwind.mon.mean.nc" rel="noreferrer" target="_blank">uwind.mon.mean.nc</a><br>
&gt; &gt; sdfopen <a href="http://vwind.mon.mean.nc" rel="noreferrer" target="_blank">vwind.mon.mean.nc</a><br>
&gt; &gt;<br>
&gt; &gt; set lon 180 320<br>
&gt; &gt; set lat 0 60<br>
&gt; &gt; set t 763 766<br>
&gt; &gt; set lev 1000 300<br>
&gt; &gt; qu=vint(1000,shu,.1*uwnd.2,300<wbr>)<br>
&gt; &gt; qv=vint(1000,shu,.1*vwnd.3,300<wbr>)<br>
&gt; &gt;<br>
&gt; &gt; vimf=-hdivg(qu,qv)<br>
&gt; &gt; d vimf<br>
&gt; &gt;<br>
&gt; &gt; ga-&gt; qu=vint(1000, shum.1*uwnd.2,600)<br>
&gt; &gt;<br>
&gt; &gt; Error from VINT:  Incompatible grids.<br>
&gt; &gt; *** glibc detected *** grads: double free or corruption (!prev): 0x0000000003707c40 ***<br>
&gt; &gt; ======= Backtrace: =========<br>
&gt; &gt; /lib/x86_64-linux-gnu/libc.so.<br>
&gt; &gt; 6(+0x7e846)[0x7f898ac4c846]<br>
&gt; &gt; grads(gree+0x2bb)[0x4fab92]<br>
&gt; &gt; grads(gadef+0x126b)[0x518144]<br>
&gt; &gt; grads(gacmd+0x16c)[0x52af24]<br>
&gt; &gt; grads(main+0xa7a)[0x4975ab]<br>
&gt; &gt; /lib/x86_64-linux-gnu/libc.so.<wbr>6(__libc_start_main+0xed)[0x7f<wbr>898abef76d]<br>
&gt; &gt; grads[0x4960a9]<br>
&gt; &gt; ======= Memory map: ========<br>
&gt; &gt; 00400000-00bcb000 r-xp 00000000 08:05 1573111<br>
&gt; &gt; 7f898bdc2000-7f898bdc3000 r--p 00017000 08:05 1572886                    /lib/x86_64-linux-gnu/libpthre<wbr>ad-2.15.soAborted (core dumped)<br>
&gt; &gt;<br>
&gt; &gt; Does anyone know how to fix it?<br>
&gt; &gt; Best,<br>
&gt; &gt; Pablo Camarena<br>
&gt; &gt;<br>
&gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; gradsusr mailing list<br>
&gt; &gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/li<wbr>stinfo/gradsusr</a><br>
&gt;<br>
&gt; --<br>
&gt; Jennifer Miletta Adams<br>
&gt; Center for Ocean-Land-Atmosphere Studies (COLA)<br>
&gt; George Mason University<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/li<wbr>stinfo/gradsusr</a><br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/li<wbr>stinfo/gradsusr</a><br>
<br>
--<br>
Jennifer Miletta Adams<br>
Center for Ocean-Land-Atmosphere Studies (COLA)<br>
George Mason University<br>
<br>
<br>
<br>
<br>
______________________________<wbr>_________________<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/li<wbr>stinfo/gradsusr</a><br>
</div></div></blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<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/li<wbr>stinfo/gradsusr</a><br>
<br></blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<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/li<wbr>stinfo/gradsusr</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br></div>