[gradsusr] Diurnal Climatology

[Mehwish]

Dear Joshua,

Thank you detailed explanation. 
I found an error in control file so I think after fixing that your given procedure will work well :)

Thank you again

With Best Regards,

Mehwish 


On 20/09/2014, at 2:34 PM, Ngaina Joshua <jngaina at gmail.com> wrote:

> 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
> 
> _______________________________________________
> 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/baaf0d96/attachment.html