[gradsusr] how to use vint function in GrADS

[sushant puranik]

Thank you for suggestion. Actually i want to calculate following integral.

                 1000

Q = (1/g) ∫ (V q)dp.

             400

I used this kind of command in grads

'define ps=1000'
'define a=vint(ps,vwnd*shum,400)'
'd a/9.81'
with this command i am able to make a plot.

whether it is correct or not ?

Thanking you


Sushant

On Fri, Jul 23, 2010 at 10:38 AM, MY <funnytyphoon at gmail.com> wrote:

> Hi,
>
> please try
>
> 'define ps=1000'
> 'd vint(ps,vwnd*shum,400)'
>
> best wishes,
> ++++++++++++++++++++++++++++++++++++++++++++++++++++
> YING, Ming
> Shanghai Typhoon Institute of China Meteorological Administration
> 166 Puxi Road, Xuhui District, Shanghai, 200030
> PR China
> ++++++++++++++++++++++++++++++++++++++++++++++++++++
>
>
>
>
>
> 2010/7/23 sushant puranik <sushantpuranik at gmail.com>:
> > Hi
> > I am trying to find vertical integral in GrADS using "vint" command. For
> > that purpose i opened vwnd and specific humidity file as first and second
> > respectively. The q file command for both the variable are as
> >
> > ga-> q file
> > File 1 : mean daily NMC reanalysis (2008)
> >   Descriptor: /home/Sushant/Desktop/10may08v.nc
> >   Binary: /home/Sushant/Desktop/10may08v.nc
> >   Type = Gridded
> >   Xsize = 29  Ysize = 25  Zsize = 8  Tsize = 1
> >   Number of Variables = 1
> >     vwnd 8 -999 mean Daily V wind
> >
> > ga-> q file 2
> > File 2 : mean daily NMC reanalysis (2008)
> >   Descriptor: /home/Sushant/Desktop/10may08q.nc
> >   Binary: /home/Sushant/Desktop/10may08q.nc
> >   Type = Gridded
> >   Xsize = 29  Ysize = 25  Zsize = 8  Tsize = 1
> >   Number of Variables = 1
> >     shum 8 -999 mean Daily specific humidity
> >
> > after this i used
> >
> > ga->  d vint(1000,vwnd*shum.2,400)
> > gxstart=0, gxstop=29, gystart=0, gystop=25.
> > gzstart=0, gzstop=1, gtstart=0, gtstop=1.
> > gxstart=0, gxstop=29, gystart=0, gystop=25.
> > gzstart=0, gzstop=1, gtstart=0, gtstop=1.
> > Error from VINT:  Incompatible grids.
> > Contouring: -0 to 2.4e+36 interval 3e+35
> > Segmentation fault
> >
> > after showing this message, the grads window gets closed.
> >
> > Can anyone has solution for this perticular problem.
> >
> > Thanks in advance
> >
> >
> > --
> > Sushant
> >
> >
> > _______________________________________________
> > 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/20100723/4be8ea92/attachment-0003.html