Oho! I meant to say,<br>>>I<b style="background-color:rgb(255,0,0)">T</b> will be nice if this issue is addressed......... :-)<br><br>I am always nice to developers!!<br><br><div class="gmail_quote">On Wed, Oct 17, 2012 at 11:55 AM, Muthuvel Chelliah <span dir="ltr"><<a href="mailto:muthuvel.chelliah@noaa.gov" target="_blank">muthuvel.chelliah@noaa.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thanks Krishnamohan.<br><br>oops.. Yes, I do see now that you had raised the same issue before :-)<br><br>I will be nice if this issue is addressed, since especially after all I made<br>
sure that I am looking at only the 365 days in ALL years with the same <br>
365_day_calendar option (ignoring the etc 1 day in all intervening leap years).<br>after all, the climo is already calculated in the first step. <br>>>set t 1 365<br>define dailyclimo=ave(dailydata,t+0,t=M,365) <<<br>
and this above step works. The problem is in the modify seasonal part,<br>which just repeats the same first climo values in the first 365 days (already<br>there!!) to the successive 365 days chunks. Right?<br><br>So, it is just like the 12 months case. Hence, technically, this also should be feasible to do,<br>
if the developers are kind enough to consider, and make it part of the standard package, since so <br>many may benefit from this. But then, easy for me to say though ..:-)<br><br>Thanks<br>Muthu<br><br><div class="gmail_quote">
On Tue, Oct 16, 2012 at 11:41 PM, krishnamohan <span dir="ltr"><<a href="mailto:krishmet@gmail.com" target="_blank">krishmet@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I have asked a similar question some time ago in grads user forum. The link is given below<br><br><a href="http://gradsusr.org/pipermail/gradsusr/2010-December/012072.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/2010-December/012072.html</a><br>
<br>and the answer from Jennifer Adams, IGES/COLA is quoted below<br><br>"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."<br>
<br><br>-- <br>Krishnamohan.K.S<br>CSIR Research Fellow<br>Department Of Atmospheric Sciences<br>Cochin University of Science and Technology<br>Cochin,India<br><br><br><br><br><div class="gmail_quote">On Tue, Oct 16, 2012 at 10:48 PM, Muthuvel Chelliah <span dir="ltr"><<a href="mailto:muthuvel.chelliah@noaa.gov" target="_blank">muthuvel.chelliah@noaa.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Grads Gurus,<br><br>Here is something to ponder.<br>In grads if we have long monthly data, we all know<br>we can compute monthly climo by, for example,<br><br>set t 1 12<br>define monclimo=ave(monthlydata,t+0,t=N,12)<br>
set t 1 N<br>modify monclimo seasonal.<br><br>In the same way, if I want to compute a quick and dirty daily climo,<br>and if had daily data for many years, I did,<br><br>set t 1 365<br>define dailyclimo=ave(dailydata,t+0,t=M,365)<br>
<br>the above step works, and the values look correct.<br>But the following step does not work!!<br><br>modify dailyclimo seasonal<br><br>It says: <br>MODIFY Error: Invalid time scaling in defined variable<br><br>Any ideas/suggestions? or is this not possible to do :-(<br>
<br>Thanks in advance!!<br><br>Muthu<span><font color="#888888"><br>-- <br><pre cols="72"><font style="font-family:tahoma,sans-serif">------------------------------------------------------------------<br>Dr.Muthuvel Chelliah, <a href="mailto:eMail%3AMuthuvel.Chelliah@noaa.gov" target="_blank">eMail:Muthuvel.Chelliah@noaa.gov</a><br>
Rm.3008, Climate Prediction Center, Ph:<a href="tel:%28301%29%20683-3393" value="+13016833393" target="_blank">(301) 683-3393</a><br>NCEP/NWS/NOAA/U.S.Dept.of Commerce,<br>5830, University Research Court, NCWCP<br>College Park, MD 20740.</font><br>
</pre><br>
</font></span><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><br>-- <br><pre cols="72"><font style="font-family:tahoma,sans-serif">------------------------------------------------------------------<br>
Dr.Muthuvel Chelliah, <a href="mailto:eMail%3AMuthuvel.Chelliah@noaa.gov" target="_blank">eMail:Muthuvel.Chelliah@noaa.gov</a><br>
Rm.3008, Climate Prediction Center, Ph:<a href="tel:%28301%29%20683-3393" value="+13016833393" target="_blank">(301) 683-3393</a><br>NCEP/NWS/NOAA/U.S.Dept.of Commerce,<br>5830, University Research Court, NCWCP<br>College Park, MD 20740.</font><br>
</pre><br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><pre cols="72"><font style="font-family:tahoma,sans-serif">------------------------------------------------------------------<br>Dr.Muthuvel Chelliah, <a href="mailto:eMail%3AMuthuvel.Chelliah@noaa.gov" target="_blank">eMail:Muthuvel.Chelliah@noaa.gov</a><br>
Rm.3008, Climate Prediction Center, Ph:(301) 683-3393<br>NCEP/NWS/NOAA/U.S.Dept.of Commerce,<br>5830, University Research Court, NCWCP<br>College Park, MD 20740.</font><br></pre><br>