>I am using calcMonAnomLLLT to calculate monthly anomalies for a variable
>from its monthly climatology. There are missing values in both the data
>files for the same grid points. So for these grid points the anomaly is
>calculated to be 0. I would like these points to be missing values (9e20 in
>this case) rather than 0. Is it possible to do this using the calcMonAnom
>function?
Hello
NCL automatically sets all expressions where
one component has a value _FillValue attribute
of the array being processed.
Does your data (say, x) have the attribute _FillValue
set to 9e20 ? If not, you can readily set it via
x@_FillValue = 9e20
--------------------------------------
Most commonly, let x(time,lev,lat,lon)
x@_FillValue = 9e20 ; if not already set
; ntim bust be divisible by 12
; reorder using named dimensions
xClm = clmMonLLLT( x(lev|:,lat|:,lon|:,time|:) ) ; xClm(lev,lat,lon)
xAnom= calcMonAnomLLLT (x, xClm)
printVarSummary(xClm)
-------------------------------------
good luck
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Tue Jan 25 2005 - 16:27:55 MST