Re: [ferret_users] Climatology for each month at hourly resolution

[Allison Hogikyan - NOAA Affiliate <allison.hogikyan@xxxxxxxx>]

Ansley did answer my initial question, but to clarify:

The end goal is to calculate the monthly diurnal cycle anomaly, and add it to the smoothed daily climatology.  The data has a diurnal cycle that is comparable to the seasonal cycle, and I don't want to smooth that out of the climatology.

My sticking point was picking out just the Januaries (Februaries, etc.) from each year efficiently, especially since I'm using the gregorian calendar.  

However, the monthly-average amplitude you mention may be a quicker means to the same end.  I would be interested to try it, if you have a sample handy!

On Fri, Jul 8, 2016 at 12:32 PM, Billy Kessler <william.s.kessler@xxxxxxxx> wrote:

Hi Allison -

Aside from the syntax error Ansley pointed out, I don't see how this computes a diurnal cycle.

It's unclear what result you want.

Is it:
- the monthly-average amplitude of the diurnal cycle, for each month? (a monthly time series)
- the RMS variability within each day? (a daily time series, perhaps grouped by month)

Please clarify, then we'd be better able to find a solution. I have a good way to do the first of these, no looping.

Billy

> On Jul 8, 2016, at 8:05 AM, Allison Hogikyan - NOAA Affiliate <allison.hogikyan@xxxxxxxx> wrote:
>
> Hi,
> I am trying to include the diurnal cycle (which changes seasonally) in a climatology.  I have 11 years of data, and am trying to calculate the average diurnal cycle in each month.  I'm having trouble figuring out the best way to construct this loop in Ferret, though.  I think there should be a way to call VAR[t="1-jan-'YEAR' 00:00:00":"31-jan-'YEAR' 23:00:00], where you loop through the 11 years in variable YEAR.
>
> This is the error I am getting:
> repeat/range=2004:2014/name=m (let januaries=alb[t="1-jan-'m':00:00:00":"31-jan-'m':23:00:00":1];save/append/file=jan.cdf januaries)
> !-> REPEAT: M:2004
>  **ERROR: command syntax: 1-jan-'m':00:00:00
> LIST/FORMAT=CDF/append/file=jan.cdf januaries
> Command file, command group, or REPEAT execution aborted
>
>
> I appreciate any help or other ideas!
>
> Allison

--

Hollings Scholar Class of 2015

University of Michigan

College of Engineering

BSE candidate, Climate Science