[gradsusr] daily climo akin to monthly climo?

[Muthuvel Chelliah]

Thanks Krishnamohan.

oops.. Yes, I do see now that you had raised the same issue before :-)

I will be nice if this issue is addressed, since especially after all I made
sure that I am looking at only the 365 days in ALL years with the same
365_day_calendar option (ignoring the etc 1 day in all intervening leap
years).
after all, the climo is already calculated in the first step.
>>set t 1 365
define dailyclimo=ave(dailydata,t+0,t=M,365) <<
and this above step works. The problem is in the modify seasonal part,
which just repeats the same first climo values in the first 365 days
(already
there!!) to the successive 365 days chunks. Right?

So, it is just like the 12 months case. Hence, technically, this also
should be feasible to do,
if the developers are kind enough to consider, and make it part of the
standard package, since so
many may benefit from this. But then, easy for me to say though ..:-)

Thanks
Muthu

On Tue, Oct 16, 2012 at 11:41 PM, krishnamohan <krishmet at gmail.com> wrote:

> I have asked a similar question some time ago in grads user forum. The
> link is given below
>
> http://gradsusr.org/pipermail/gradsusr/2010-December/012072.html
>
> and the answer from Jennifer Adams, IGES/COLA is quoted below
>
> "The 'modify' command does not calculate a daily climatology, and there is
> no way to do that directly in GrADS. You must use your own code and make
> decisions about how to handle leap days, filtering, and other
> complications."
>
>
> --
> Krishnamohan.K.S
> CSIR Research Fellow
> Department Of Atmospheric Sciences
> Cochin University of Science and Technology
> Cochin,India
>
>
>
>
> On Tue, Oct 16, 2012 at 10:48 PM, Muthuvel Chelliah <
> muthuvel.chelliah at noaa.gov> wrote:
>
>> Dear Grads Gurus,
>>
>> Here is something to ponder.
>> In grads if we have long monthly data, we all know
>> we can compute monthly climo by, for example,
>>
>> set t 1 12
>> define monclimo=ave(monthlydata,t+0,t=N,12)
>> set t 1 N
>> modify monclimo seasonal.
>>
>> In the same way, if I want to compute a quick and dirty daily climo,
>> and if had daily data for many years, I did,
>>
>> set t 1 365
>> define dailyclimo=ave(dailydata,t+0,t=M,365)
>>
>> the above step works, and the values look correct.
>> But the following step does not work!!
>>
>> modify dailyclimo seasonal
>>
>> It says:
>> MODIFY Error:  Invalid time scaling in defined variable
>>
>> Any ideas/suggestions? or is this not possible to do :-(
>>
>> Thanks in advance!!
>>
>> Muthu
>> --
>>
>> ------------------------------------------------------------------
>> Dr.Muthuvel Chelliah, eMail:Muthuvel.Chelliah at noaa.gov
>>
>>
>>
>> Rm.3008, Climate Prediction Center, Ph:(301) 683-3393
>> NCEP/NWS/NOAA/U.S.Dept.of Commerce,
>> 5830, University Research Court, NCWCP
>> College Park, MD 20740.
>>
>>
>>
>> _______________________________________________
>> 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
>
>


-- 

------------------------------------------------------------------
Dr.Muthuvel Chelliah, eMail:Muthuvel.Chelliah at noaa.gov
Rm.3008, Climate Prediction Center, Ph:(301) 683-3393
NCEP/NWS/NOAA/U.S.Dept.of Commerce,
5830, University Research Court, NCWCP
College Park, MD 20740.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20121017/26841b1e/attachment-0003.html