<div dir="ltr">Dear Sir Dennis and NCL -experts,<div><br></div><div>Thank you for this! </div><div><br></div><div>I have one more question regarding fourier_info and
ezfftf in NCL.</div><div><br></div><div>Is it correct that If I want to get the smoothed time series for the first, second, and third harmonics, I have to use the <b>ezfftf</b>?</div><div>(I would like to get the smoothed time series but I can only think of the ezfftf and save to a netcdf file)</div><div><br></div><div>Does this mean that the fourier_info for the first three harmonics is applied on the following smoothened time series?</div><div> cf(:,:,:, 0:1 ) = 0.0 ; First harmonic<br></div><div> cf(:,:,:, 0:2 ) = 0.0 ; Second harmonic<br></div><div> cf(:,:,:, 0:3 ) = 0.0 ; Third harmonic<br></div><div><br></div><div>Code is below</div><div>-----------------------------------------------------</div><div><br></div><div> Y = y(lat|:,lon|:,time|:) ; convenience<br> printVarSummary(Y)<br> cf = ezfftf (Y) ; ==> cf(2,nlat,mlon,ntim/2)<br> printVarSummary(cf) ; [2] x [73] x [144] x [36]<br> ;print(cf@xbar)<br> do n=0,nlat*mlon-1<br> cf@xbar(n) = 0.0 ; reconstruct with harmonics only<br> end do<br> cf(:,:,:, 0:3 ) = 0.0 ; set non-desired harmonics to 0.0<br></div><div>------------------------------------------------------------------------------------</div><div><br></div><div><br></div><div>Thank you so much for the help.</div><div><br></div><div>Sincerely,</div><div>Lyndz</div><div><br></div><div><br></div><div><br><div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><table border="0" cellspacing="0" cellpadding="0" style="border-collapse:collapse"><tbody></tbody></table></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Feb 19, 2022 at 6:36 AM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><a href="https://www.ncl.ucar.edu/Applications/fouranal.shtml" target="_blank"><b>https://www.ncl.ucar.edu/Applications/fouranal.shtml</b></a></div><div><br></div><div>See the example<b>. fanal_4.ncl ...</b> in particular, the figure.<br><b></b></div><div><b><br></b></div><div>The phase component gives the index (subscript) of the maxima.<b><br></b></div><div><pre> (<b>1</b>,0) 72 <=== Phases: <b>location of first maximum</b>
(<b>1</b>,1) 36
(<b>1</b>,2) 12<br>for 'latitude', then the first maxima occur at <br>finfo(1,lat(72)), finfo(1,lat(36)), finfo(1,lat(12))
</pre><pre>If this were 'time', then the first maxima occur at <br>finfo(1,time(72)), finfo(1,time(36)), finfo(1,time(12))
</pre></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Feb 17, 2022 at 7:42 AM Lyndz via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear NCL-experts,</div><div><br></div><div>I am using the <b>fourier_info </b>to extract the amplitude of the first and second harmonics of a 24 time series (<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/fourier_info.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/fourier_info.shtml</a>)</div><div><br></div><div><b>;[time | 24] x [lon | 42] x [lat | 68]</b><br></div><div><b>finfo_pa = fourier_info(pa(lon|:,lat|:,time|:),2,1.)<br></b></div><div><br></div><div>But, I am only getting the amplitude, phase, and variance explained from the output. </div><div><br></div><div>How can I extract the specific time/s (for first harmonics, it is one time, for the 2nd harmonic it is two) when the wave reaches its peak (or the time when I got the amplitude) for the first and second harmonics?</div><div><br></div><div>I'll appreciate any help on this.</div><div><br></div><div>Sincerely,</div><div>Lyndz</div><div><br></div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div>