[gradsusr] Dear users
rajuattada
raju_attadamsc at rediffmail.com
Fri Jul 23 06:19:50 EDT 2010
Dear Grads Users,
I am Raju A.
I have IMD gridded rainfall.
I want to calculate Standard deviation for JJAS.
How can i do.
Please help me
Thanks in advance
RAJU A
Research Fellow
Indian Institute Of Tropical Meteorology
Dr. Homi Bhabha Road,
PUNE-411008
On Fri, 23 Jul 2010 15:10:54 +0530 wrote
>Send gradsusr mailing list submissions to
>gradsusr at gradsusr.org
>
>To subscribe or unsubscribe via the World Wide Web, visit
>http://gradsusr.org/mailman/listinfo/gradsusr
>or, via email, send a message with subject or body 'help' to
>gradsusr-request at gradsusr.org
>
>You can reach the person managing the list at
>gradsusr-owner at gradsusr.org
>
>When replying, please edit your Subject line so it is more specific
>than "Re: Contents of gradsusr digest..."
>
>
>Today's Topics:
>
> 1. Re: how to use vint function in GrADS (sushant puranik)
> 2. Re: how to use vint function in GrADS (MY)
> 3. Re: how to use vint function in GrADS (MY)
> 4. Re: Obtaining total values of a dimension on a netCDFfile
> (Becker, Bernd)
>
>
>----------------------------------------------------------------------
>
>Message: 1
>Date: Fri, 23 Jul 2010 11:04:01 +0530
>From: sushant puranik <sushantpuranik at gmail.com>
>Subject: Re: [gradsusr] how to use vint function in GrADS
>To: GrADS Users Forum <gradsusr at gradsusr.org>
>Message-ID:
><AANLkTingGqD37gATrPBoOtCYi0XXLwcYd8LgPcC+=MjW at mail.gmail.com>
>Content-Type: text/plain; charset="utf-8"
>
>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-0001.html
>
>------------------------------
>
>Message: 2
>Date: Fri, 23 Jul 2010 15:52:37 +0800
>From: MY <funnytyphoon at gmail.com>
>Subject: Re: [gradsusr] how to use vint function in GrADS
>To: GrADS Users Forum <gradsusr at gradsusr.org>
>Message-ID:
><AANLkTi=_Cf-t2P0qUkjDL3Q0rd2ETxXQvmkuFuf-dkax at mail.gmail.com>
>Content-Type: text/plain; charset=GB2312
>
>Hi, Sushant,
>
>I think 'd a' is enough. Please refer to the document for vint().
>
>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>:
>> 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
>>
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
>
>------------------------------
>
>Message: 3
>Date: Fri, 23 Jul 2010 15:52:37 +0800
>From: MY <funnytyphoon at gmail.com>
>Subject: Re: [gradsusr] how to use vint function in GrADS
>To: GrADS Users Forum <gradsusr at gradsusr.org>
>Message-ID:
><AANLkTi=_Cf-t2P0qUkjDL3Q0rd2ETxXQvmkuFuf-dkax at mail.gmail.com>
>Content-Type: text/plain; charset=GB2312
>
>Hi, Sushant,
>
>I think 'd a' is enough. Please refer to the document for vint().
>
>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>:
>> 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
>>
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
>
>------------------------------
>
>Message: 4
>Date: Fri, 23 Jul 2010 09:37:41 +0100
>From: "Becker, Bernd"
>Subject: Re: [gradsusr] Obtaining total values of a dimension on a
>netCDFfile
>To: "GrADS Users Forum" <gradsusr at gradsusr.org>
>Message-ID:
><E209E2066317644FB35E5AA0E8C1FB7F079994 at EXXMAIL02.desktop.frd.metoffice.com>
>
>Content-Type: text/plain; charset="us-ascii"
>
>Ok.
>you are getting there.
>
>try:
>set t 1
>
>d sum(data,t=1,t=32)
>
>any numbers looking familiar?
>Cheers,
>Bernd.
>
>
>________________________________
>
>From: gradsusr-bounces at gradsusr.org
>[mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Thiago Veloso
>Sent: 22 July 2010 18:06
>To: GrADS Users Forum
>Subject: Re: [gradsusr] Obtaining total values of a dimension on a
>netCDFfile
>
>
> Dear Bernd,
>
> Thanks for your reply.
>
> After trying the procedure you suggested I received lots of
>information on display. Each information set corresponding to a specific
>day.
>
> For example, setting t 1 3 generates:
>
>ga-> set t 1 3
>Time values set: 2009:1:1:0 2009:1:3:0
>ga-> d cmorph
>Data Type = grid
>Dimensions = 0 1
>I Dimension = 0 to 4949 Linear -0.0363783 0.0727567
>J Dimension = 1 to 1649 Linear -59.9636 0.0727714
>Sizes = 4950 1649 8162550
>Undef value = -9.99e+08
>Undef count = 114789 Valid count = 8047761
>Min, Max = 0 241
>Cmin, cmax, cint = 0 240 20
>Stats[sum,sumsqr,root(sumsqr),n]: 4.20946e+06 8.7492e+07 9353.71
>8047761
>Stats[(sum,sumsqr,root(sumsqr))/n]: 0.523059 10.8716 3.29721
>Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 0.523059 10.8716 3.29721
>Stats[(sigma,var)(n)]: 3.25546 10.598
>Stats[(sigma,var)(n-1)]: 3.25546 10.598
>Data Type = grid
>Dimensions = 0 1
>I Dimension = 0 to 4949 Linear -0.0363783 0.0727567
>J Dimension = 1 to 1649 Linear -59.9636 0.0727714
>Sizes = 4950 1649 8162550
>Undef value = -9.99e+08
>Undef count = 89326 Valid count = 8073224
>Min, Max = 0 217
>Cmin, cmax, cint = 0 200 20
>Stats[sum,sumsqr,root(sumsqr),n]: 3.78218e+06 7.38863e+07 8595.71
>8073224
>Stats[(sum,sumsqr,root(sumsqr))/n]: 0.468485 9.15202 3.02523
>Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 0.468485 9.15202 3.02523
>Stats[(sigma,var)(n)]: 2.98874 8.93254
>Stats[(sigma,var)(n-1)]: 2.98874 8.93254
>Data Type = grid
>Dimensions = 0 1
>I Dimension = 0 to 4949 Linear -0.0363783 0.0727567
>J Dimension = 1 to 1649 Linear -59.9636 0.0727714
>Sizes = 4950 1649 8162550
>Undef value = -9.99e+08
>Undef count = 177477 Valid count = 7985073
>Min, Max = 0 200
>Cmin, cmax, cint = 0 200 20
>Stats[sum,sumsqr,root(sumsqr),n]: 3.8412e+06 9.12816e+07 9554.14
>7985073
>Stats[(sum,sumsqr,root(sumsqr))/n]: 0.481048 11.4315 3.38105
>Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 0.481048 11.4315 3.38106
>Stats[(sigma,var)(n)]: 3.34666 11.2001
>Stats[(sigma,var)(n-1)]: 3.34666 11.2001
>ga->
>
> Is this information enough?
>
> Best wishes,
>
> Thiago.
>
>
>--- On Thu, 22/7/10, Becker, Bernd
>wrote:
>
>
>
>From: Becker, Bernd
>Subject: Re: [gradsusr] Obtaining total values of a dimension on
>a netCDF file
>To: "GrADS Users Forum" <gradsusr at gradsusr.org>
>Date: Thursday, 22 July, 2010, 12:20
>
>
>sdfopen file.nc
>set gxout stats
>
>set t 1 32
>d var
>
>what does that produce?
>
>Bernd.
>
>________________________________
>
>From: gradsusr-bounces at gradsusr.org
>[mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Thiago Veloso
>Sent: 22 July 2010 15:18
>To: gradsusr at gradsusr.org
>Subject: [gradsusr] Obtaining total values of a dimension on a
>netCDF file
>
>
> Dear GrADS colleagues,
>
> I am working with a netCDF file which contains precipitation values
>for 31 days. What I am trying to do is to use NCO to obtain the sum of
>the "time" dimension (thus, the 31 days), but with no success until the
>present moment.
>
> Could anybody give some directions?
>
> Best regards,
>
> Thiago.
>
>
>-----Inline Attachment Follows-----
>
>
>_______________________________________________
>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/ed4915ee/attachment.html
>
>------------------------------
>
>_______________________________________________
>gradsusr mailing list
>gradsusr at gradsusr.org
>http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>End of gradsusr Digest, Vol 5, Issue 50
>***************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100723/48dfc5d7/attachment-0003.html
More information about the gradsusr
mailing list