Sujata,<br><br>Assume the file names for non-leap year & leap year climatology as:<br>filename.leap.year.climo<br>filename.non-leap.year.climo<br><br>set up soft-links for non-leap years in a do loop as (in cshell) as:<br>
<br>foreach year (1979 1981 1982 1983 1985 1986 1987 1989 1990 1991)<br>ln -s filename.non-leap.year.climo climatology.$year<br>end<br><br>similarly set up soft-links for leap years in a do loop as (in cshell) as:<br><br>
foreach year (1980 1984 1988 1992)<br>
ln -s filename.leap.year.climo climatology.$year<br>
end<br><br>add as many non-leap/leap years as necessary within the parentheses for the data <br>period in your analysis. Note: the newly created common (for both leap & non-leap year)<br>climatology file name for each YEAR of analysis is:<br>
climatology.YEAR<br><br>now for the grads control file use for dset as:<br><br>dset climatology.%y4<br>options template<br>....<br>....<br>...<br><br>I hope this helps. If it is still not year, please talk to someone<br>who knows unix and grads. I am myself not an expert in these <br>
things either :-)<br><br>Muthu<br><br><div class="gmail_quote">On Mon, Sep 10, 2012 at 2:38 AM, Sujata Mandke <span dir="ltr"><<a href="mailto:amin@tropmet.res.in" target="_blank">amin@tropmet.res.in</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Respecetd Muthuvel sir,<br>
<br>
Thanks a lot for help.<br>
<br>
I need further assistance in this problem, as I understood<br>
your suggestions partly. Follows the details of the part of<br>
your suggestions that I understood and the part that<br>
is not clear.<br>
<br>
> Create two different climatology data sets, one for<br>
> and 365 days another for leap year with 366 days.<br>
<br>
This is clear.<br>
<br>
> Set up repeated soft links for each year of data set respectively<br>
>to non-leap years and >leap years in analysis period (very easy<br>
> to do, couple of lines of unix scripts). Create a >single grads<br>
> ctl file or all these climo years.<br>
<br>
I do not know how to create soft links (in unix scripts)<br>
and how to create a single grads ctl file for two separate datasets<br>
(one for non-leap years climo and other for leap years climo).<br>
<br>
>Open up your total field control field and this climo control file,<br>
> then you are in business. We did like this at CPC more than<br>
> two decades ago and we continue to do like this now as we work<br>
> with daily cdas fields. I hope this helps. Please email me<br>
> separately if you need further assistance. Good day!<br>
<br>
Opening two 'ctl' files and calculating difference between<br>
control fields and climo field looks easier to do in grads script.<br>
<br>
<br>
Regards and bets wishes<br>
Dr. Sujata Mandke<br>
Scientist,IITM,Pune-8<br>
<br>
----- Original Message -----<br>
From: Muthuvel Chelliah <<a href="mailto:muthuvel.chelliah@noaa.gov">muthuvel.chelliah@noaa.gov</a>><br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Sent: Fri, 07 Sep 2012 20:22:55 +0530 (IST)<br>
Subject: Re: [gradsusr] daily anomaly<br>
<br>
Sujata,<br>
Create two different climatology data sets, one for 365 days and another<br>
for leap year with 366 days.<br>
Set up repeated soft links for each year of data set respectively to<br>
non-leap years and leap years<br>
in analysis period (very easy to do, couple of lines of unix scripts).<br>
Create a single grads ctl file<br>
for all these climo years. Open up your total field control field and this<br>
climo control file, then you are in business. We did like this at CPC more<br>
than two decades ago and we continue to do like this now as we work with<br>
daily cdas fields. I hope this helps. Please email me separately if you<br>
need further assistance. Good day!<br>
Muthu<br>
<br>
On Fri, Sep 7, 2012 at 2:52 AM, Sujata Mandke <<a href="mailto:amin@tropmet.res.in">amin@tropmet.res.in</a>> wrote:<br>
<br>
> Respected Narapusetty Sir,<br>
><br>
> Thanks for prompt reply and willingness to provide<br>
> further assistance in my problem. I will go thorugh your<br>
> article and then try to solve the problem.<br>
><br>
> However, you had solved the problem by<br>
> calculating climatology. I am having longterm daily<br>
> climatology dataset, so there is no need to calculate climatology.<br>
> Regards<br>
> Sujata Mandke<br>
> IITM,PUNE-8,INDIA.<br>
><br>
> ----- Original Message -----<br>
> From: Balachandrudu Narapusetty <<a href="mailto:bnarapus@gmu.edu">bnarapus@gmu.edu</a>><br>
> To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
> Cc: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> Sent: Thu, 06 Sep 2012 20:10:31 +0530 (IST)<br>
> Subject: Re: [gradsusr] daily anomaly<br>
><br>
> Hi Dr. Sujata,<br>
><br>
> I am happy to see somebody has same problem that I had some time ago. I<br>
> solved it by estimating the climatology (and thereafter anomalies) using<br>
> few harmonics.<br>
><br>
> I am sure this article is useful in dealing with your problem:<br>
><br>
> Narapusetty, Balachandrudu, Timothy DelSole, Michael K. Tippett, 2009:<br>
> Optimal Estimation of the Climatological Mean. J. Climate, 22, 4845–4859.<br>
> doi: <a href="http://dx.doi.org/10.1175/2009JCLI2944.1" target="_blank">http://dx.doi.org/10.1175/2009JCLI2944.1</a><br>
><br>
> I am happy to provide further assistance if you have questions.<br>
><br>
> Best regards,<br>
> Bala<br>
><br>
><br>
> ----- Original Message -----<br>
> From: Sujata Mandke <<a href="mailto:amin@tropmet.res.in">amin@tropmet.res.in</a>><br>
> Date: Thursday, September 6, 2012 7:24 am<br>
> Subject: [gradsusr] daily anomaly<br>
><br>
> > Dear Grads Users<br>
> > I wish to calculate daily OLR anomaly from<br>
> > the following datasets.<br>
> > I have daily OLR data (netcdf) for the period<br>
> > 1Jun1974-31May2012. Also long term daily OLR<br>
> > climatology (netcdf) data for 365 days.<br>
> ><br>
> > As such writting script for calculating daily<br>
> > anomaly is very simple. But how to tackle following<br>
> > two problems while calculating anomaly<br>
> ><br>
> > (i) Actual data is on Gregorian(real) calendar<br>
> > and thus has data for 366 days in leap years.<br>
> > (ii) What to do with missing values?<br>
> ><br>
> > Does anybody has grads script to calculate daily<br>
> > anomaly by taking care of above two problems.<br>
> > Regards<br>
> > Dr. Sujata Mandke,Scientist,<br>
> > IITM,PUNE-8,INDIA.<br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> ><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
><br>
<br>
<br>
<br>
--<br>
<br>
------------------------------------------------------------------<br>
Dr.Muthuvel Chelliah, <a href="mailto:eMail%3AMuthuvel.Chelliah@noaa.gov">eMail:Muthuvel.Chelliah@noaa.gov</a><br>
Rm.3008, Climate Prediction Center, Ph:<a href="tel:%28301%29%20683-3393" value="+13016833393">(301) 683-3393</a><br>
NCEP/NWS/NOAA/U.S.Dept.of Commerce,<br>
5830, University Research Court, NCWCP<br>
College Park, MD 20740.<br>
<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div><br><br clear="all"><br>-- <br><pre cols="72"><font style="font-family:tahoma,sans-serif" size="2">------------------------------------------------------------------<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>