[gradsusr] daily climo akin to monthly climo?

[Muthuvel Chelliah]

Oho! I meant to say,
>>I*T* will be nice if this issue is addressed.........  :-)

I am always nice to developers!!

On Wed, Oct 17, 2012 at 11:55 AM, Muthuvel Chelliah <
muthuvel.chelliah at noaa.gov> wrote:

> 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.
>
>
>


-- 

------------------------------------------------------------------
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/f808d224/attachment-0003.html