[gradsusr] Diurnal Climatology

[Ngaina Joshua]

Mehwish


For diurnal climatology, you have to define your hour by hour climatology
for the 26 years before plotting

try this out


#### define your hour climatology from hour 1 to hour 24
define hour1cli = ave(pratesfc*86400,t=1,t=227785,24)
define hour2cli = ave(pratesfc*86400,t=1,t=227786,24)
define hour3cli = ave(pratesfc*86400,t=1,t=227787,24)
define hour4cli = ave(pratesfc*86400,t=1,t=227788,24)
define hour5cli = ave(pratesfc*86400,t=1,t=227789,24)
define hour6cli = ave(pratesfc*86400,t=1,t=227790,24)
..
### define upto hour 24

define hour21cli = ave(pratesfc*86400,t=21,t=227805,24)
define hour22cli = ave(pratesfc*86400,t=22,t=227806,24)
define hour23cli = ave(pratesfc*86400,t=23,t=227807,24)
define hour24cli = ave(pratesfc*86400,t=24,t=227808,24)

### Your can then plot your 24 hour climatology
define hour24cli = (hour1cli+hour2cli+...............+hour24cli)

set gxout shaded
##### set clevs 0 5 15 30 60 100 250
#####set ccols 7 12 10 3 5 11 4

d hour24cli
draw title Diurnal Climatology
####   cbar

I hope it works

Regards
Joshua

On Sat, Sep 20, 2014 at 7:58 AM, mehwish ramzan <mehwish.ramzan at gmail.com>
wrote:

> Dear Jennifer,
>
> I have tried your kind suggestions but error is still the same.
>
>
>
> With Best Regards,
>
> Mehwish
>
> On Sat, Sep 20, 2014 at 4:56 AM, Jennifer Adams <jma at cola.iges.org> wrote:
>
>> I think you need to change the third argument in your ave() declaration
>> from “t+227808” to “t=227808”.
>> —Jennifer
>>
>>
>> On Sep 19, 2014, at 3:13 AM, mehwish ramzan <mehwish.ramzan at gmail.com>
>> wrote:
>>
>>
>> Dear GrADS Users,
>>
>> I am trying to calculate the Diurnal Climatology for 26 years "hourly
>> data" as follows;
>>
>> open r_pgb_hourly.ctl
>>
>> set t 1 24
>>
>> set x 1
>>
>> set y 1
>>
>> define p1= ave(pratesfc*60*60*24, t-0, t+227808, 24)
>>
>>
>>
>> Averaging.  dim = 3, start = 1, end = 227808
>>
>> but after averaging it suddenly gives *Segmentation fault error*
>>
>>
>>  Also if i want to calculate just the JJA diurnal Climatology for 26
>> years then what changes i need to make in above script?
>>
>>
>> Please guide me.
>>
>> Thank You
>>
>> With Best Regards,
>>
>> Mehwish
>>
>>
>>  _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>> --
>> Jennifer M. Adams
>> Center for Ocean-Land-Atmosphere Studies (COLA)
>> 111 Research Hall, Mail Stop 2B3
>> George Mason University
>> 4400 University Drive
>> Fairfax, VA 22030
>>
>>
>>
>>
>>
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>>
>>
>
>
> --
> Mehwish RAMZAN
> Doctoral Candidate,
> ====================================================================
>  Oki Laboratory / Yoshimura Laboratory
> Available at:
> Atmosphere and Ocean Research Institute, Kashiwa Campus, University of
> Tokyo,Japan
> *Also at:*Institute of Industrial Science, Komaba Campus, University of
> Tokyo,Japan
> phone : +81-80-4806-3962
> email: mehwish at rainbow.iis.u-tokyo.ac.jp
>
>
>
>
> _______________________________________________
> 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/20140920/08907057/attachment-0001.html