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