<html><head></head><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:13px"><div id="yui_3_16_0_ym19_1_1496188863224_31224" dir="ltr"><span id="yui_3_16_0_ym19_1_1496188863224_31406">Hi Pablo, here I attach my vertical intergarted moisture flux convergence scripts for your reference.</span></div><div id="yui_3_16_0_ym19_1_1496188863224_31224" dir="ltr">Hope it help.</div><div id="yui_3_16_0_ym19_1_1496188863224_31224" dir="ltr"><br></div><div id="yui_3_16_0_ym19_1_1496188863224_31224" dir="ltr">Thank you.</div><div id="yui_3_16_0_ym19_1_1496188863224_31224" dir="ltr"><br></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: verdana, helvetica, sans-serif; font-size: 13px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> On Wednesday, May 31, 2017 5:46 AM, Pablo Camarena <jcamarena006@gmail.com> wrote:<br></font></div> <br><br> <div class="y_msg_container"><div id="yiv3258757972"><div><div dir="ltr"><div><div><div><br clear="none"></div>Thanks Jennifer,<br clear="none"><br clear="none"></div><span class="yiv3258757972gmail-short_text" id="yiv3258757972gmail-result_box" lang="en"><span class="yiv3258757972gmail-">Now I'm trying the following:<br clear="none"><br clear="none"></span></span><div style="margin-bottom:0in;">sdfopen <a rel="nofollow" shape="rect" target="_blank" href="http://shum.mon.mean.nc/">shum.mon.mean.nc</a></div>sdfopen <a rel="nofollow" shape="rect" target="_blank" href="http://uwind.mon.mean.nc/">uwind.mon.mean.nc</a><br clear="none"></div>sdfopen <a rel="nofollow" shape="rect" target="_blank" href="http://vwind.mon.mean.nc/">vwind.mon.mean.nc</a><br clear="none"><div><a rel="nofollow" shape="rect" target="_blank" href="http://uwind.mon.mean.nc/"></a>
<div style="margin-bottom:0in;">set lon 225 227<br clear="none">set lat 23 24<br clear="none"></div>set t 182 274<br clear="none">set z 1 3<br clear="none">
<div style="margin-bottom:0in;">d vint(1000,shum.1*lterp(uwnd.2,shum.1),850)</div><div style="margin-bottom:0in;"><span class="yiv3258757972gmail-short_text" id="yiv3258757972gmail-result_box" lang="en"><span class="yiv3258757972gmail-">And I get the following error:</span></span></div>"Operation error: Invalid dimension environment<br clear="none"> Too many varying dimensions"<br clear="none"><br clear="none"></div><div>I think the error is with vertical levels.<br clear="none"></div><div><div class="yiv3258757972gmail-a3s yiv3258757972gmail-aXjCH yiv3258757972gmail-m15c5afa820bbcfab" id="yiv3258757972gmail-:kb"></div></div></div><div class="yiv3258757972gmail_extra"><br clear="none"><div class="yiv3258757972gmail_quote">2017-05-30 14:08 GMT-06:00 Jennifer M Adams <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:jadams21@gmu.edu" target="_blank" href="mailto:jadams21@gmu.edu">jadams21@gmu.edu</a>></span>:<br clear="none"><div class="yiv3258757972yqt9728563319" id="yiv3258757972yqt61599"><blockquote class="yiv3258757972gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi, Pablo —<br clear="none">
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 clear="none">
<br clear="none">
ga-> d vint(1000,shum.1*lterp(uwnd.2, shum.1),600)<br clear="none">
<br clear="none">
In the script lines you posted, I noticed that the syntax here is incorrect:<br clear="none">
> qu=vint(1000,shu,.1*uwnd.2, 300)<br clear="none">
> qv=vint(1000,shu,.1*vwnd.3, 300)<br clear="none">
<br clear="none">
change "shu,.1" to "shum.1"<br clear="none">
<br clear="none">
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 rel="nofollow" shape="rect" target="_blank" href="ftp://cola.gmu.edu/grads/scripts/zinterp.gs">ftp://cola.gmu.edu/grads/ 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 rel="nofollow" shape="rect" target="_blank" href="http://cola.gmu.edu/grads/gadoc/gradfuncvint.html">http://cola.gmu.edu/grads/ 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 clear="none">
—Jennifer<br clear="none">
<div><div class="yiv3258757972h5"><br clear="none">
<br clear="none">
On May 30, 2017, at 2:18 PM, Pablo Camarena <<a rel="nofollow" shape="rect" ymailto="mailto:jcamarena006@gmail.com" target="_blank" href="mailto:jcamarena006@gmail.com">jcamarena006@gmail.com</a>> wrote:<br clear="none">
<br clear="none">
> Hi,<br clear="none">
><br clear="none">
> I want to calculate Vertically integarted moisture flux convergence (VIMFC)<br clear="none">
><br clear="none">
> But when I give the following instruction it marks me errors:<br clear="none">
><br clear="none">
> sdfopen <a rel="nofollow" shape="rect" target="_blank" href="http://shum.mon.mean.nc/">shum.mon.mean.nc</a><br clear="none">
> sdfopen <a rel="nofollow" shape="rect" target="_blank" href="http://uwind.mon.mean.nc/">uwind.mon.mean.nc</a><br clear="none">
> sdfopen <a rel="nofollow" shape="rect" target="_blank" href="http://vwind.mon.mean.nc/">vwind.mon.mean.nc</a><br clear="none">
><br clear="none">
> set lon 180 320<br clear="none">
> set lat 0 60<br clear="none">
> set t 763 766<br clear="none">
> set lev 1000 300<br clear="none">
> qu=vint(1000,shu,.1*uwnd.2, 300)<br clear="none">
> qv=vint(1000,shu,.1*vwnd.3, 300)<br clear="none">
><br clear="none">
> vimf=-hdivg(qu,qv)<br clear="none">
> d vimf<br clear="none">
><br clear="none">
> ga-> qu=vint(1000, shum.1*uwnd.2,600)<br clear="none">
><br clear="none">
> Error from VINT: Incompatible grids.<br clear="none">
> *** glibc detected *** grads: double free or corruption (!prev): 0x0000000003707c40 ***<br clear="none">
> ======= Backtrace: =========<br clear="none">
> /lib/x86_64-linux-gnu/libc.so.<br clear="none">
> 6(+0x7e846)[0x7f898ac4c846]<br clear="none">
> grads(gree+0x2bb)[0x4fab92]<br clear="none">
> grads(gadef+0x126b)[0x518144]<br clear="none">
> grads(gacmd+0x16c)[0x52af24]<br clear="none">
> grads(main+0xa7a)[0x4975ab]<br clear="none">
> /lib/x86_64-linux-gnu/libc.so. 6(__libc_start_main+0xed)[ 0x7f898abef76d]<br clear="none">
> grads[0x4960a9]<br clear="none">
> ======= Memory map: ========<br clear="none">
> 00400000-00bcb000 r-xp 00000000 08:05 1573111<br clear="none">
> 7f898bdc2000-7f898bdc3000 r--p 00017000 08:05 1572886 /lib/x86_64-linux-gnu/ libpthread-2.15.soAborted (core dumped)<br clear="none">
><br clear="none">
> Does anyone know how to fix it?<br clear="none">
> Best,<br clear="none">
> Pablo Camarena<br clear="none">
><br clear="none">
</div></div>> ______________________________ _________________<br clear="none">
> gradsusr mailing list<br clear="none">
> <a rel="nofollow" shape="rect" ymailto="mailto:gradsusr@gradsusr.org" target="_blank" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br clear="none">
> <a rel="nofollow" shape="rect" target="_blank" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/ listinfo/gradsusr</a><br clear="none">
<br clear="none">
--<br clear="none">
Jennifer Miletta Adams<br clear="none">
Center for Ocean-Land-Atmosphere Studies (COLA)<br clear="none">
George Mason University<br clear="none">
<br clear="none">
<br clear="none">
<br clear="none">
<br clear="none">
______________________________ _________________<br clear="none">
gradsusr mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:gradsusr@gradsusr.org" target="_blank" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/ listinfo/gradsusr</a><br clear="none">
</blockquote></div></div><br clear="none"></div></div></div><div class="yqt9728563319" id="yqt38015">_______________________________________________<br clear="none">gradsusr mailing list<br clear="none"><a shape="rect" ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br clear="none"><a shape="rect" href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br clear="none"></div><br><br></div> </div> </div> </div></div></body></html>