| Since it looks like you have more than 1 year of data, you can get ferret to do the monthly averaging for you, using the original data directly: set dat climatological_axes; can dat climatological_axes let mydata_cl = mydata[gt=month_irreg@mod] shade mydata_cl This is especially useful if your data time axis is irregular, or extends over partial years. In those cases, the average of mydata_cl: mydata_cl[l=@ave] is a better estimate of the mean than the simple mean: mydata[l=@ave] ! this avoids weighting parts of the year with more observation. @mod makes the axis extend indefinitely, thus allowing plotting the climatology over the original data: plot mydata ! more than 1 year, possibly irregular plot/over mydata_cl Or finding anomalies: let mydata_cl2 = mydata_cl[gt=mydata] ! put the climatology onto the original time axis let mydata_anom = mydata - mydata_cl2 Billy K Sent from my phone .... Please excuse my brevity .... and my typing.
|