[gradsusr] vertically integrated moisture flux

Fri Jun 2 03:27:05 EDT 2017

thanks
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 ?)

On Fri, Jun 2, 2017 at 6:14 AM, Sridhara Nayak <
sridharanayakiitkgp at gmail.com> wrote:

> How about the following?
>
> sdfopen shum.nc
> sdfopen uwnd.nc
> sdfopen vwnd.nc
>
> p1000 = const(shum.1,1000,-a)
> define qu=vint(p1000,shum.1*uwnd.2,100)
> define qv=vint(p1000,shum.1*vwnd.3,100)
> d -hdivg(qu,qv)
>
> Best Regards
> Sridhara
>
>
> On Fri, Jun 2, 2017 at 10:01 AM, afwande juliet <afwandej965 at gmail.com>
> wrote:
>
>> Dear users of grads
>> did someone manage to plot vertically integrated moisture flux.
>> I have tried the function vint but gets the below error
>> http://cola.gmu.edu/grads/gadoc/gradfuncvint.html
>> ga-> q file
>> File 1 : mean daily NMC reanalysis (1997)
>>   Descriptor: shum.1997.nc
>>   Binary: shum.1997.nc
>>   Type = Gridded
>>   Xsize = 144  Ysize = 73  Zsize = 8  Tsize = 365  Esize = 1
>>   Number of Variables = 1
>>      shum  8  t,z,y,x  mean Daily specific humidity
>> ga-> set t 1
>> Time values set: 1997:1:1:0 1997:1:1:0
>> ga-> set z 1
>> LEV set to 1000 1000
>> ga-> set z 2
>> LEV set to 925 925
>> ga-> set z 3
>> LEV set to 850 850
>> ga-> set z 8
>> LEV set to 300 300
>> ga-> d vint(850,shum.1*lterp(uwnd.2,shum.1),600)
>> Error from VINT:  Incompatible grids.
>> Contouring: -0.25 to 0.2 interval 0.05
>> *** Error in `grads': free(): invalid pointer: 0x000000000145a430 ***
>> ======= Backtrace: =========
>> /lib/x86_64-linux-gnu/libc.so.6(+0x777e5)[0x7f41ce33c7e5]
>> /lib/x86_64-linux-gnu/libc.so.6(+0x7fe0a)[0x7f41ce344e0a]
>> /lib/x86_64-linux-gnu/libc.so.6(cfree+0x4c)[0x7f41ce34898c]
>> grads(gagfre+0x5a)[0x449e5a]
>> grads(gagrel+0x58)[0x493be8]
>> grads(gadspl+0x668)[0x4ab318]
>> grads(gacmd+0x17ae)[0x4aa33e]
>> grads(main+0x62e)[0x41277e]
>> /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f41ce2e5830]
>> grads(_start+0x29)[0x412e79]
>> ======= Memory map: ========
>> 00400000-004e8000 r-xp 00000000 08:05 524842
>> /usr/bin/grads
>> 006e7000-006e8000 r--p 000e7000 08:05 524842
>> /usr/bin/grads
>> 006e8000-006ec000 rw-p 000e8000 08:05 524842
>> /usr/bin/grads
>> 006ec000-00992000 rw-p 00000000 00:00 0
>> 012a4000-014b8000 rw-p 00000000 00:00 0
>>
>> On Wed, May 31, 2017 at 1:00 AM, Jennifer M Adams <jadams21 at gmu.edu>
>> wrote:
>>
>>> Usage Note #5: The vint function operates only in an X-Y varying
>>> dimension environment.
>>>
>>> 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:
>>>
>>> set z 1
>>> set t 182
>>> d vint(etc. etc.)
>>>
>>> If that works and you are sure your expression is what you want, then
>>> you can
>>>
>>> set t 182 274
>>> define yourvname = vint(etc. etc.)
>>>
>>> … and GrADS will loop over T under the hood to create your defined
>>> variable that varies X, Y, and T.
>>> —Jennifer
>>>
>>>
>>> On May 30, 2017, at 5:43 PM, Pablo Camarena <jcamarena006 at gmail.com>
>>> wrote:
>>>
>>> >
>>> > Thanks Jennifer,
>>> >
>>> > Now I'm trying the following:
>>> >
>>> > sdfopen shum.mon.mean.nc
>>> > sdfopen uwind.mon.mean.nc
>>> > sdfopen vwind.mon.mean.nc
>>> > set lon 225 227
>>> > set lat 23 24
>>> > set t 182 274
>>> > set z 1 3
>>> > d vint(1000,shum.1*lterp(uwnd.2,shum.1),850)
>>> > And I get the following error:
>>> > "Operation error:  Invalid dimension environment
>>> >   Too many varying dimensions"
>>> >
>>> > I think the error is with vertical levels.
>>> >
>>> > 2017-05-30 14:08 GMT-06:00 Jennifer M Adams <jadams21 at gmu.edu>:
>>> > Hi, Pablo —
>>> > 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:
>>> >
>>> > ga-> d vint(1000,shum.1*lterp(uwnd.2,shum.1),600)
>>> >
>>> > In the script lines you posted, I noticed that the syntax here is
>>> incorrect:
>>> > > qu=vint(1000,shu,.1*uwnd.2,300)
>>> > > qv=vint(1000,shu,.1*vwnd.3,300)
>>> >
>>> > change "shu,.1" to "shum.1"
>>> >
>>> > 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 —ftp://cola.gmu.edu/grads/scripts/zinterp.gs.  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 (http://cola.gmu.edu/grads/gadoc/gradfuncvint.html) 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.
>>> > —Jennifer
>>> >
>>> >
>>> > On May 30, 2017, at 2:18 PM, Pablo Camarena <jcamarena006 at gmail.com>
>>> wrote:
>>> >
>>> > > Hi,
>>> > >
>>> > > I want to calculate Vertically integarted moisture flux convergence
>>> (VIMFC)
>>> > >
>>> > > But when I give the following instruction it marks me errors:
>>> > >
>>> > > sdfopen shum.mon.mean.nc
>>> > > sdfopen uwind.mon.mean.nc
>>> > > sdfopen vwind.mon.mean.nc
>>> > >
>>> > > set lon 180 320
>>> > > set lat 0 60
>>> > > set t 763 766
>>> > > set lev 1000 300
>>> > > qu=vint(1000,shu,.1*uwnd.2,300)
>>> > > qv=vint(1000,shu,.1*vwnd.3,300)
>>> > >
>>> > > vimf=-hdivg(qu,qv)
>>> > > d vimf
>>> > >
>>> > > ga-> qu=vint(1000, shum.1*uwnd.2,600)
>>> > >
>>> > > Error from VINT:  Incompatible grids.
>>> > > *** glibc detected *** grads: double free or corruption (!prev):
>>> 0x0000000003707c40 ***
>>> > > ======= Backtrace: =========
>>> > > /lib/x86_64-linux-gnu/libc.so.
>>> > > 6(+0x7e846)[0x7f898ac4c846]
>>> > > grads(gree+0x2bb)[0x4fab92]
>>> > > grads(gadef+0x126b)[0x518144]
>>> > > grads(gacmd+0x16c)[0x52af24]
>>> > > grads(main+0xa7a)[0x4975ab]
>>> > > /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed)[0x7f
>>> 898abef76d]
>>> > > grads[0x4960a9]
>>> > > ======= Memory map: ========
>>> > > 00400000-00bcb000 r-xp 00000000 08:05 1573111
>>> > > 7f898bdc2000-7f898bdc3000 r--p 00017000 08:05 1572886
>>>     /lib/x86_64-linux-gnu/libpthread-2.15.soAborted (core dumped)
>>> > >
>>> > > Does anyone know how to fix it?
>>> > > Best,
>>> > > Pablo Camarena
>>> > >
>>> > > _______________________________________________
>>> > > gradsusr mailing list
>>> > > gradsusr at gradsusr.org
>>> > > http://gradsusr.org/mailman/listinfo/gradsusr
>>> >
>>> > --
>>> > Jennifer Miletta Adams
>>> > Center for Ocean-Land-Atmosphere Studies (COLA)
>>> > George Mason University
>>> >
>>> >
>>> >
>>> >
>>> > _______________________________________________
>>> > gradsusr mailing list
>>> > gradsusr at gradsusr.org
>>> > http://gradsusr.org/mailman/listinfo/gradsusr
>>> >
>>> > _______________________________________________
>>> > gradsusr mailing list
>>> > gradsusr at gradsusr.org
>>> > http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>> --
>>> Jennifer Miletta Adams
>>> Center for Ocean-Land-Atmosphere Studies (COLA)
>>> George Mason University
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170602/4ddf89ff/attachment-0001.html 