[ncl-talk] NCL fourier Script

[Omondi Augustine]

Dear NCL users,

I am new to NCL and currently running fourier analysis of daily(3hrly)
rainfall estimate on NCL version 6.2.1 on a cygwin. I have modified the NCL
script given in example section as fanal_2 to read;

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"
   f  = addfile("/climatology.nc","r")
   t  = f->cmorph_precip

    N = dimsizes(t)
   nlat = N(0)
   mlon = N(1)
   time = N(2)
    nhar     = 2
  finfo_t = fourier_info (t(lat|:,lon|:,time|:), nhar, 1.)
   copy_VarCoords(t(0,:,:), finfo_t(0,:,:,0))
   finfo_t at long_name = "CMORPH precipitation estimate"
   finfo_t at units     = "mm/hr"
   printVarSummary(finfo_t)

   finfo_t at _FillValue = t at _FillValue
   finfo_t(0,:,:,:) = where(finfo_t(0,:,:,:).eq.0, t at _FillValue,
finfo_t(0,:,:,:))
   print(finfo_t)

The question is that i am not sure whether i am doing the right thing and
therefore kindly asking for experts opinion on the result(amplitude, phase
and percentage variance) it is giving. Attached is my data(long term mean
of 3hrly rainfall estimate on a 20lat by 20lon grid).

Thanks in advance for your advice

Omondi Augustine

Graduate student
College of Atmospheric Science, Nanjing University of Information Science
and Technology, Nanjing 210044, Jiangsu, China
Email: dimosh05 at gmail.com / omondi at kms.or.ke
phone: +86 132 60771131
Skype: omondi_augustine
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161104/05b086d3/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: climatology.nc
Type: application/x-netcdf
Size: 29780 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161104/05b086d3/attachment.nc