[gradsusr] vertically integrated moisture flux

Thu Jun 8 10:05:25 EDT 2017

Dear users
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.
Can someone help me. I want  to plot moisure flux at different times i,e

   1.   start of the season
   2. peak of the season
   3. cessation of the season etc . That means time must vary

'p925 = const(shum.1,925,-a)'
'define qu=vint(p925,shum.1*uwnd.2,850)'
'define qv=vint(p925,shum.1*vwnd.3,850)'
'define quave=ave(qu,t=1,t=2)'
'define qvave=ave(qv,t=1,t=2)'

see my script and figures plotted which shows that even if I vary time, the
moisture is the same

thanks

On Fri, Jun 2, 2017 at 10:27 AM, afwande juliet <afwandej965 at gmail.com>
wrote:

> 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/20170608/7111783e/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mflux_sond.eps
Type: application/postscript
Size: 697365 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20170608/7111783e/attachment-0001.eps 