[ncl-talk] Perform Fourier Analysis to Pentad Series
Lyndon Mark Olaguera
olagueralyndonmark429 at gmail.com
Tue Jan 3 22:28:16 MST 2017
Dear Sir Dennis,
Many many thanks for the reply.
I got it working already. :-)
Sincerely,
*Lyndz*
On Wed, Jan 4, 2017 at 6:52 AM, Dennis Shea <shea at ucar.edu> wrote:
> ntim = 73
> nlat = 73
> mlon = 144
>
> y = random_uniform(0,100,(/ntim,nlat,mlon/))
> y!0 = "time"
> y!1 = "lat"
> y!2 = "lon"
>
> printVarSummary(y) ; [time | 73] x [lat | 73] x [lon |
> 144]
> print("==========")
>
> Y = y(lat|:,lon|:,time|:) ; convenience
> printVarSummary(Y)
> cf = ezfftf (Y) ; ==> cf(2,nlat,mlon,ntim/2)
> printVarSummary(cf) ; [2] x [73] x [144] x [36]
>
> ;print(cf at xbar)
>
> do n=0,nlat*mlon-1
> cf at xbar(n) = 0.0 ; reconstruct with harmonics only
> end do
>
> cf(:,:,:, 0:4 ) = 0.0 ; set non-desired harmonics to 0.0
> cf(:,:,:, 6:16) = 0.0
> cf(:,:,:,18: ) = 0.0
>
> YY = ezfftb(cf,cf at xbar)
> copy_VarMeta(Y, YY)
> YY at info = "Only harmonics 5 and 18 retained: means removed"
> printVarSummary(YY) ; (lat,lon,time)
>
>
> You can apply to your file.
> Also, *you* cn determine how to write the file out.
>
> Method 1 at
> http://www.ncl.ucar.edu/Applications/o-netcdf.shtml
>
>
> On Mon, Jan 2, 2017 at 7:41 PM, Lyndon Mark Olaguera <
> olagueralyndonmark429 at gmail.com> wrote:
>
>> Dear All,
>>
>> I would like to perform fourier analysis to a pentad (nonoverlapping
>> five day average) time series(t=1 to t=73). I would like to get the sum of
>> the 5th and 18th harmonics and output as a separate netcdf file. Is this
>> possible to do in ncl?I am new in ncl. Ill appreciate if anyone can point
>> me to the right functions that I can use here.
>>
>> Here is the link to the file that I'm working with.
>>
>> https://drive.google.com/file/d/0B9faET7Bc2o8cWpmRE1yNldNVE9
>> SQ1JNNF9XdXJPVUdqbUc0/view?usp=sharing
>>
>> Many thanks in advance,
>>
>> Lyndz
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170104/b939521d/attachment.html
More information about the ncl-talk
mailing list