[ncl-talk] Fwd: Extracting Fourier-transformed data from wkSpaceTime function

Dennis Shea shea at ucar.edu
Wed Oct 7 12:48:47 MDT 2020


The following was sent offline to me and Karin Meier-Fleischer.
The 2nd email states that the issue was solved using an (unfortunately)
undocumented option

   opt at netCDF = True

---------- Forwarded message ---------
From: Mckinnon-Gray, Elliot <elliot.mckinnon-gray at mpimet.mpg.de>
Date: Tue, Oct 6, 2020 at 11:01 AM
Subject: Extracting Fourier-transformed data from wkSpaceTime function



Hi Karin (Cc Dennis because of the time of day I tend to work),

What I am trying to do is to calculate correlation and difference between
the spectra of modelled and observed precipitation. For this I need to
write to disk and save Fourier-transformed data. So as far as my
understanding goes, in the function script, data are subjected to a 2D FFT
and then this transformed field is power spectra in wavenumber frequency
space, and the plots are created using this data and according axes.

The function script I refer to I *think* is called wkspacetime_cam.. maybe
you will be able to tell from the attached script. I no longer remember
since I work with a modified version which I will attach - this originated
from Karin actually since it was modified to change psumanti 0.0 values to
_Fillvalue [stops it getting stuck in an infinite loop for certain data]

I’m not particularly familiar with FFT in NCL, or temp file creation, so
would I be able to ask for some guidance as to where in the script I might
find the transformed data and how I can write that to a file? I hope I
communicated my question clearly!

Thanks,

Elliot
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201007/ea94462c/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: function_wkSpacetime.ncl
Type: application/octet-stream
Size: 35321 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201007/ea94462c/attachment.obj>


More information about the ncl-talk mailing list