<div dir="ltr"><div><div><br></div><div>[1] <br>As noted in the four_info documentation, the dimension being used "must be periodic". 'time' is not periodic.<br></div>This function would work on a *global* (time,lat,lon) array because the longitudes are periodic.<br></div>So, specifically, what you are doing is *not* correct.<br><div><br>[2] <br></div><div>The spatial domain is small: 3.4S-1.4N and 30.6E-35.4E<br></div><div>You have only 8 time steps of 3-hrly data. Hence, one day of data. <br><br></div><div>I suggest simple averaging of all grid points at each time step.<br><br></div><div> cm = f->cmorph_precip <br></div><div> printVarSummary(cm) ; (time,lat,lon) ... (8,20,20); dim numbers (0,1,2)<br><br></div><div> cmAvg = dim_avg_n_Wrap(cm,(/1,2/)) ; cmAvg(time), cmAvg(8) <br></div><div> print(cmAvg)<br><br></div><div>There are only 8 values .... do you see any discernible patterm. EG: diurnal cycle ?<br></div><div>plot the eight value<br><br><br>[3]<br>ncl-talk would like to help BUT ...<br>ncl-talk can not take the place of student-advisor interaction. <br></div><div>I think this is the correct approach for you.<br><br></div><div>Good Luck<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 4, 2016 at 12:34 AM, Omondi Augustine <span dir="ltr"><<a href="mailto:dimosh05@gmail.com" target="_blank">dimosh05@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear NCL users,<br><br></div>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;<br><br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl" <br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"<br> f = addfile("/<a href="http://climatology.nc" target="_blank">climatology.nc</a>","r")<wbr> <br> t = f->cmorph_precip<br><br> N = dimsizes(t)<br> nlat = N(0) <wbr> <br> mlon = N(1)<br> time = N(2) <wbr> <br> nhar = 2 <wbr> <br> finfo_t = fourier_info (t(lat|:,lon|:,time|:), nhar, 1.) <br> copy_VarCoords(t(0,:,:), finfo_t(0,:,:,0)) <br> finfo_t@long_name = "CMORPH precipitation estimate"<br> finfo_t@units = "mm/hr"<br> printVarSummary(finfo_t)<br> <br> finfo_t@_FillValue = t@_FillValue<br> finfo_t(0,:,:,:) = where(finfo_t(0,:,:,:).eq.0, t@_FillValue, finfo_t(0,:,:,:)) <br> print(finfo_t)<br><div><br></div><div>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).<br><br></div><div>Thanks in advance for your advice<br></div><div><br></div><div><div><div><div class="m_-6649930979404777292gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><font face="arial, helvetica, sans-serif" size="2">Omondi Augustine</font></div><div dir="ltr"><font face="arial, helvetica, sans-serif" size="2"><br>Graduate student</font></div><div dir="ltr"><font face="arial, helvetica, sans-serif" size="2"><span style="text-align:center;line-height:115%">College of Atmospheric Science, </span><span style="text-align:center;line-height:115%">Nanjing
University of Information Science and Technology, Nanjing 210044, Jiangsu,
China</span></font></div><div dir="ltr"><font face="arial, helvetica, sans-serif" size="2">Email: <a href="mailto:dimosh05@gmail.com" target="_blank">dimosh05@gmail.com</a> / <a href="mailto:omondi@kms.or.ke" target="_blank">omondi@kms.or.ke</a><br>phone: <a href="tel:%2B86%20132%2060771131" value="+8613260771131" target="_blank">+86 132 60771131</a></font></div><div><font face="arial, helvetica, sans-serif" size="2">Skype: omondi_augustine</font></div></div></div></div></div></div>
</div></div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>