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

Dennis Shea shea at ucar.edu
Wed Oct 7 15:35:18 MDT 2020


When I wrote this library, I tried to add options I speculated users might
find useful. There were lots of options! I guess I neglected to document
the netCDF option. Unfortunately, I no longer have access to the NCL
documentation.

Stay Healthy!

On Wed, Oct 7, 2020 at 9:29 AM Mckinnon-Gray, Elliot <
elliot.mckinnon-gray at mpimet.mpg.de> wrote:

> Replying for information if this is referred to on ncl-talk…
>
>
> I solved this issue and managed to extract an .nc file with all my power
> spectra data simply by writing
> opt at netCDF = True in the wkSpaceTime.ncl function (or whichever version
> you have)
>
> This option I only found after combing through the function from
> diagnostics_cam - might I suggest that is quickly put onto one of the
> information pages for these WK scripts? I couldn’t find information on this
> option at least.
>
> Elliot
>
> On 6. Oct 2020, at 19:01, Mckinnon-Gray, Elliot <
> elliot.mckinnon-gray at mpimet.mpg.de> wrote:
>
> 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
>
> <function_wkSpacetime.ncl>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201007/5fd3de11/attachment.html>


More information about the ncl-talk mailing list