Re: 6-hour climatologies

From: Maria Gehne <mgehne_at_nyahnyahspammersnyahnyah>
Date: Mon Sep 16 2013 - 10:46:44 MDT

Hi Mike,

I don't think there is a script. But it's straightforward to compute the
climatology yourself. For example, if your data is X(time, lat, lon),
then you could try something like this:

dims = dimsizes(X)
nlat = dims(1)
nlon = dims(2)
ndays = 366 (or 365, depending on how many days your year has)

Xclim = new( (/4*ndays,nlat,nlon/), typeof(X) )

do i = 0,4*ndays-1
     Xclim(i,:,:) = dim_avg_n( X(i::4*ndays,:,:) , 0 )
end do

And then copy meta data if you need to.

Hope this helps,
Maria

On 09/16/2013 10:11 AM, Michael Ventrice wrote:
> Hello,
>
> I am trying to calculate 6-hour climatologies of surface temperature
> and am wondering if there is a script to help with this? I have
> realized that using the clmDayTLL function is just for daily-averaged
> fields.
>
> Thanks,
> Mike
>
> --
> * Michael J. **Ventrice, Ph.D. *
> Medium Range Operations Scientist
> MVentrice@wsi.com <mailto:MVentrice@wsi.com>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
Maria Gehne
Postdoctoral Fellow
NCAR/CGD/CAS
P.O. Box 3000
Boulder, CO 80307

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Sep 16 10:46:53 2013

This archive was generated by hypermail 2.1.8 : Mon Sep 16 2013 - 13:43:52 MDT