<div dir="ltr"><div><div><div><br></div>Thanks Jennifer,<br><br></div><span id="gmail-result_box" class="gmail-short_text" lang="en"><span class="gmail-">Now I'm trying the following:<br><br></span></span><p style="margin-bottom:0in">sdfopen <a target="_blank" href="http://shum.mon.mean.nc/">shum.mon.mean.nc</a></p>sdfopen <a target="_blank" href="http://uwind.mon.mean.nc/">uwind.mon.mean.nc</a><br></div>sdfopen <a href="http://vwind.mon.mean.nc">vwind.mon.mean.nc</a><br><div><a target="_blank" href="http://uwind.mon.mean.nc/"></a>
<p style="margin-bottom:0in">set lon 225 227<br>set lat 23 24<br></p>set t 182 274<br>set z 1 3<br>
<p style="margin-bottom:0in">d vint(1000,shum.1*lterp(uwnd.2,shum.1),850)</p><p style="margin-bottom:0in"><span class="gmail-short_text" id="gmail-result_box" lang="en"><span class="gmail-">And I get the following error:</span></span></p>"Operation error: Invalid dimension environment<br> Too many varying dimensions"<br><br></div><div>I think the error is with vertical levels.<br></div><div><div class="gmail-a3s gmail-aXjCH gmail-m15c5afa820bbcfab" id="gmail-:kb"><wbr></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-05-30 14:08 GMT-06:00 Jennifer M Adams <span dir="ltr"><<a href="mailto:jadams21@gmu.edu" target="_blank">jadams21@gmu.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi, Pablo —<br>
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>
<br>
ga-> d vint(1000,shum.1*lterp(uwnd.2,<wbr>shum.1),600)<br>
<br>
In the script lines you posted, I noticed that the syntax here is incorrect:<br>
> qu=vint(1000,shu,.1*uwnd.2,<wbr>300)<br>
> qv=vint(1000,shu,.1*vwnd.3,<wbr>300)<br>
<br>
change "shu,.1" to "shum.1"<br>
<br>
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/<wbr>scripts/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/<wbr>gadoc/gradfuncvint.html</a>) to make sure you’ve got everything right. If you’re still having trouble, please also include the output from "ncdump -c” on your data files, the version of GrADS you are using, and the shortest possible script that illustrates the error.<br>
—Jennifer<br>
<div><div class="h5"><br>
<br>
On May 30, 2017, at 2:18 PM, Pablo Camarena <<a href="mailto:jcamarena006@gmail.com">jcamarena006@gmail.com</a>> wrote:<br>
<br>
> Hi,<br>
><br>
> I want to calculate Vertically integarted moisture flux convergence (VIMFC)<br>
><br>
> But when I give the following instruction it marks me errors:<br>
><br>
> sdfopen <a href="http://shum.mon.mean.nc" rel="noreferrer" target="_blank">shum.mon.mean.nc</a><br>
> sdfopen <a href="http://uwind.mon.mean.nc" rel="noreferrer" target="_blank">uwind.mon.mean.nc</a><br>
> sdfopen <a href="http://vwind.mon.mean.nc" rel="noreferrer" target="_blank">vwind.mon.mean.nc</a><br>
><br>
> set lon 180 320<br>
> set lat 0 60<br>
> set t 763 766<br>
> set lev 1000 300<br>
> qu=vint(1000,shu,.1*uwnd.2,<wbr>300)<br>
> qv=vint(1000,shu,.1*vwnd.3,<wbr>300)<br>
><br>
> vimf=-hdivg(qu,qv)<br>
> d vimf<br>
><br>
> ga-> qu=vint(1000, shum.1*uwnd.2,600)<br>
><br>
> Error from VINT: Incompatible grids.<br>
> *** glibc detected *** grads: double free or corruption (!prev): 0x0000000003707c40 ***<br>
> ======= Backtrace: =========<br>
> /lib/x86_64-linux-gnu/libc.so.<br>
> 6(+0x7e846)[0x7f898ac4c846]<br>
> grads(gree+0x2bb)[0x4fab92]<br>
> grads(gadef+0x126b)[0x518144]<br>
> grads(gacmd+0x16c)[0x52af24]<br>
> grads(main+0xa7a)[0x4975ab]<br>
> /lib/x86_64-linux-gnu/libc.so.<wbr>6(__libc_start_main+0xed)[<wbr>0x7f898abef76d]<br>
> grads[0x4960a9]<br>
> ======= Memory map: ========<br>
> 00400000-00bcb000 r-xp 00000000 08:05 1573111<br>
> 7f898bdc2000-7f898bdc3000 r--p 00017000 08:05 1572886 /lib/x86_64-linux-gnu/<wbr>libpthread-2.15.soAborted (core dumped)<br>
><br>
> Does anyone know how to fix it?<br>
> Best,<br>
> Pablo Camarena<br>
><br>
</div></div>> ______________________________<wbr>_________________<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/<wbr>listinfo/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">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
</blockquote></div><br></div>