[ncl-talk] calculate monthly anomaly

[Dennis Shea]

Attached is an *un*documented and *un*supported function that (hopefully)
will 'pad' a series.

If you have x(lat,lon,time) .... it will have to be reordered because the
function currently pad only the leftmost dimension.

load "dim_pad.ncl"

      .....
   xr = x(time|:,lat|:,lon|:)

   X = *dim_pad*(xr, *12, 0*)  ; xr[*164*,nlat,mlon] ===> dim_pad(xr,12, 0)
===> X[*168* ][nlat][mlon]]

    x_clm = *clmMonTLL*(X)    ; *(12,lat,lon)*
    printVarSummary(x_clm)
    delete([/X, xr/])             ; delete temporary/work arrays

On Wed, Dec 13, 2017 at 1:55 PM, Debasish Hazra <debasish.hazra5 at gmail.com>
wrote:

> Hi,
> I have 14 years of monthly satelitte data, but there are some missing
> months in between (total 164 files, inplace of 168). I am trying to
> calculate monthly climatology and anomaly for each of the months. But as I
> am using either clmMonLLT or calcMonAnomLLT it needs "time dimension must
> be multiple of 12". Is there any other function available to calculate
> anomalies when time dimension is not exactly divisible by 12.
> Thanks.
> Debasish
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171213/45bf6a5b/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dim_pad.ncl
Type: application/octet-stream
Size: 2532 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171213/45bf6a5b/attachment.obj>