[gradsusr] Areal Average over Depth

[Andrew Friedman]

To clarify: Does ‘mean' account for non-linear spacing in the vertical?

Otherwise, this would only return the vertical mean if your depth levels are regularly spaced. You can use ‘vint' for vertical integration.

-Andrew

> On Aug 21, 2018, at 4:08 PM, Joseph Basconcillo <jbasconcillo at pagasa.dost.gov.ph> wrote:
> 
> I would like to share the reply of Davide as reference to similar issues. Thank you.
> 
> set x 1
> set y 1
> set z 1
> set t 1 last
> d tloop( mean(aave(var, global),z=10,z=20) 
> 
> Joseph
> 
> 
> ----- Original Message -----
> From: "Joseph Basconcillo" <jbasconcillo at pagasa.dost.gov.ph>
> To: "GrADS Users Forum" <gradsusr at gradsusr.org>, "davide sacchetti" <davide.sacchetti at arpal.gov.it>
> Sent: Tuesday, August 21, 2018 10:58:27 PM
> Subject: Re: [gradsusr] Areal Average over Depth
> 
> Hi, David:
> 
> It did not result to averaged values from 100-300m depth. Instead, it resulted hovmoller diagram that varies with depth and time. I want integration of subsurface ocean temperature from 100-300m depth levels. Thanks.
> 
> Yours,
> Joseph
> 
> ----- Original Message -----
> From: "Davide Sacchetti" <davide.sacchetti at arpal.gov.it>
> To: gradsusr at gradsusr.org
> Sent: Tuesday, August 21, 2018 5:32:29 PM
> Subject: Re: [gradsusr] Areal Average over Depth
> 
> try this:
> 
> define test=tloop(aave(var, global))
> d test
> 
> it should work
> bye
> Davide
> 
> On Tue, 2018-08-21 at 15:20 +0800, Joseph Basconcillo wrote:
>> Good day! I want to create a timeseries data of subsurface ocean
>> temperature from 100-300m using global domain. I suppose the function
>> aave can produce timeseries in one depth level. Using global domain
>> does not make sense but I want to test it first then proceed to
>> regional domain.
>> 
>> Say, if I do..
>> 
>> set x 1
>> set y 1
>> set z 10 20
>> set t 1 last
>> d tloop(aave(var, global))
>> 
>> Error says the Z cannot vary. So I want to take the average value of
>> Z from 10th-20th depth level. Is there a way to do this in GrADS? I
>> looked at the function vint. How can I used that in my problem? Thank
>> you.
>> 
>> Yours,
>> Joseph Basconcillo
>> 
>> _______________________________________________
>> 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
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr