Hallo NCL,
Anybody can tell me why did I get the constant values using the output from
fourier filter:
(0) X2 contains all constant values at nx = 14
(0) X2 contains all constant values at nx = 29
(0) X2 contains all constant values at nx = 44
(0) X2 contains all constant values at nx = 59
(0) X2 contains all constant values at nx = 74
(0) X2 contains all constant values at nx = 89
etc
--
Here is my mini script:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
in = addfile("hgt.anom.nc","r") ; climatology: 1989-2005
z = in->Z300 ; (time,lat,lon)
dimz = dimsizes(z)
ntim = dimz(0)
nlat = dimz(1)
mlon = dimz(2)
;****************************************
CF = ezfftf (z)
printVarSummary(CF)
cf = CF
cf(:,:,:,1:) = 0.0
z_wave1 = ezfftb (cf, 0.0)
copy_VarMeta(z, z_wave1)
X2 = onedtond( ndtooned(z_wave1), (/ *nxy*, *ntim* /) )
if(any(ismissing(X2))) then
print("X2 contains one or more missing values, cannot continue!!.")
return
end if
xxx = dimsizes(X2)
nx = xxx(0)
nt = xxx(1)
do i=0,nx-1
pmin = min(X2(i,:))
pmax = max(X2(i,:))
if(pmin.eq.pmax) then
print("X2 contains all constant values at nx = " + i)
end if
end do
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Aug 4 07:48:22 2013
This archive was generated by hypermail 2.1.8 : Wed Aug 07 2013 - 13:59:00 MDT