[ncl-talk] calculate monthly anomaly

Debasish Hazra debasish.hazra5 at gmail.com
Tue Dec 26 10:12:35 MST 2017


Thanks Dennis. I tried to use "dim_pad" and it gave the following errors :

fatal:syntax error: line 49 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/dim_pad.ncl before or near then
  elseif (rankx.eq.2) then
-------------------------^
fatal:syntax error: possibly an undefined procedure
fatal:syntax error: line 57 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/dim_pad.ncl before or near then
  elseif (rankx.eq.3) then
-------------------------^
fatal:syntax error: possibly an undefined procedure
fatal:syntax error: line 65 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/dim_pad.ncl before or near then
  elseif (rankx.eq.4) then
-------------------------^
fatal:syntax error: possibly an undefined procedure
fatal:syntax error: line 73 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/dim_pad.ncl before or near then
  elseif (rankx.eq.5) then
-------------------------^
fatal:syntax error: possibly an undefined procedure
fatal:Syntax Error in block, block not executed
fatal:error at line 91 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/dim_pad.ncl

Debasish

On Wed, Dec 13, 2017 at 7:18 PM, Dennis Shea <shea at ucar.edu> wrote:

> 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/20171226/6b1351e9/attachment.html>


More information about the ncl-talk mailing list