<div dir="ltr">Alright, so I was thinking, if I could (seems easy enough with CDO) get each of the years 1981, 1982... into single files, I can then do this:<div><br></div><div><div><font size="1">Assume the file names for non-leap year &amp; leap year climatology as:</font></div><div><font size="1">filename.leap.year.climo</font></div><div><font size="1">filename.non-leap.year.climo</font></div><div><font size="1"><br></font></div><div><font size="1">set up soft-links for non-leap years in a do loop as (in cshell) as:</font></div><div><font size="1"><br></font></div><div><font size="1">foreach year  (1979 1981 1982 1983 1985 1986 1987 1989 1990 1991)</font></div><div><font size="1">ln -s filename.non-leap.year.climo  climatology.$year</font></div><div><font size="1">end</font></div><div><font size="1"><br></font></div><div><font size="1">similarly set up soft-links for leap years in a do loop as (in cshell) as:</font></div><div><font size="1"><br></font></div><div><font size="1">foreach year  (1980 1984 1988 1992)</font></div><div><font size="1">ln -s filename.leap.year.climo  climatology.$year</font></div><div><font size="1">end</font></div><div><font size="1"><br></font></div><div><font size="1">add as many non-leap/leap years as necessary within the parentheses for</font></div><div><font size="1">the data</font></div><div><font size="1">period in your analysis. Note: the newly created common (for both leap &amp;</font></div><div><font size="1">non-leap year)</font></div><div><font size="1">climatology file name for each YEAR of analysis is:</font></div><div><font size="1">climatology.YEAR</font></div><div><font size="1"><br></font></div><div><font size="1">now for the grads control file use for dset as:</font></div><div><font size="1"><br></font></div><div><font size="1">dset climatology.%y4</font></div><div><font size="1">options template</font></div></div><div><br></div><div><br></div><div>----</div><div><br></div><div>open currentrun.ctl</div><div>open climodaa.ctl</div><div><br></div><div>define tmp2mcur = tmp2m.1</div><div>define tmp2mcli = <span style="font-size:12.8px">ave(tmp2m.2,t=1,t=43228,</span><span style="font-size:12.8px">1460)</span></div><div><br></div><div>Which would at least get me SOMEWHERE but I&#39;m still confused on how then I would do: <span style="font-size:12.8px">ave(tmp2m.2,t=1,t=43228,</span><span style="font-size:12.8px">1460)&#39; -or- </span><span style="font-size:12.8px">ave(hgt500mb,t=1,t=43228,</span><span style="font-size:12.8px">1464)&#39; during the leap years to get tmp2m from climo files for Nov 22 00z? Because my original point of the t=1 in that example needs to be the current dates zhour value no? Nov 22 00z is 1301 on a non-leapyear or 1305 during a leapyear... even if you used that in a loop t=i (current loop iteration) wouldn&#39;t do anything magical either.</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Now, unless there&#39;s some sort of crazy magical unicorn stuff happening where-by opening the current model runs .ctl file first, t=1 suddenly becomes Nov 22 00z, t=2 would be Nov 22 06z... I can&#39;t seem to find any information on that.</span></div></div>