[ncl-talk] calculate monthly anomaly

Dennis Shea shea at ucar.edu
Tue Dec 26 11:59:57 MST 2017


Function uses 'elseif' introduced in 6.4.0. 

Sent from my iPhone

> On Dec 26, 2017, at 12:12 PM, Debasish Hazra <debasish.hazra5 at gmail.com> wrote:
> 
> 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 undocumented and unsupported 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/ab181002/attachment.html>


More information about the ncl-talk mailing list