[MITgcm-support] High-frequency time series output at a single domain point for computation of a power spectrum

[Burns, Paul]

Hi Martin,

Apologies for the delay, I've been distracted by other tasks.  I've just got around to sorting this out again.  I have a couple of remaining questions:

  *   If I set nSetRegMskFile=1, and so compute my statistics at only 1 level, how do I choose the level?  I can't see an obvious way to choose the level.
  *
  *   To try and understand what is happening, I downloaded the binary file regMask_lat24.bin and read it into Python (using dat = np.fromfile('regMask_lat24.bin', dtype='>f4') ).  The array 'dat' has length Nx*Ny*2, where Nx and Ny are the number of points along x and y coordinates, respectively.  First of all, following your advice I was expecting dat to have length Nx*Ny..?  The values in dat are strange, with non-integers appearing, whereas I was expecting only 0s and 1s.  So either I have read in the binary incorrectly, or I am missing some information about how this works.


Many thanks,
Paul





Dr Paul Burns MInstP FRMetS

Research Fellow

Geophysical and Astrophysical Fluid Dynamics Group

Department of Mathematics

Harrison Building

University of Exeter

________________________________
From: MITgcm-support <mitgcm-support-bounces at mitgcm.org> on behalf of Martin Losch <Martin.Losch at awi.de>
Sent: 27 January 2022 07:00
To: mitgcm-support at mitgcm.org <mitgcm-support at mitgcm.org>
Subject: Re: [MITgcm-support] High-frequency time series output at a single domain point for computation of a power spectrum

CAUTION: This email originated from outside of the organisation. Do not click links or open attachments unless you recognise the sender and know the content is safe.


The mask array “regMask_lat24.bin” (it’s in input) has the same horizontal dimensions as the model domain, plus a 3rd dimension of the length nSetRegMskFile. In the examples this is just one and there is only 1 level (you can see that because in input the 2D bathyme​try files have the same size as regMask_lat24.bin). In theory you could have multiple levels in case you want to do averages of overlapping regions.

Martin

________________________________
From: MITgcm-support <mitgcm-support-bounces at mitgcm.org> on behalf of Burns, Paul <P.Burns2 at exeter.ac.uk>
Sent: Wednesday, January 26, 2022 11:07
To: MITgcm-support at mitgcm.org
Subject: Re: [MITgcm-support] High-frequency time series output at a single domain point for computation of a power spectrum

Hi Martin,

That's great.  It looks like the stats part of the Diagnostics package is worth a try.  It might be worth me looking at time series at different isolated domain points in the future.

I had a look at the global_ocean.cs32x15 example.  Do you know if the mask array should have the same dimensions as the model domain?

Many thanks for your help,
Paul




Dr Paul Burns MInstP FRMetS

Research Fellow

Geophysical and Astrophysical Fluid Dynamics Group

Department of Mathematics

Harrison Building

University of Exeter

________________________________
From: MITgcm-support <mitgcm-support-bounces at mitgcm.org> on behalf of Martin Losch <Martin.Losch at awi.de>
Sent: 25 January 2022 18:30
To: MITgcm Support <mitgcm-support at mitgcm.org>
Subject: Re: [MITgcm-support] High-frequency time series output at a single domain point for computation of a power spectrum

CAUTION: This email originated from outside of the organisation. Do not click links or open attachments unless you recognise the sender and know the content is safe.


Hi Paul,

you can use the statistics part of the diagnostics package and specify a mask for your spefict grid point. The masking is used in verification/global_ocean.cs32x15/input.seaice/data.diagnostics (DIAG_STATIS_PARMS) , but I agree the diagnostics package is not always straightforward.

If it’s just this one-time thing, I would probably just add a print statement at the right place (o:

Martin

> On 25. Jan 2022, at 13:17, Burns, Paul <P.Burns2 at exeter.ac.uk> wrote:
>
> Hi,
>
> I am using MITgcm and subscribed to the mailing list.  I need to write out timeseries at high time frequency and at a single point in the domain, in order to compute a power spectrum.  I am currently using the Diagnostics package to write out state fields, but I am unclear how to modify my setup to write out timeseries data at a single point.
>
> I would be grateful if someone could advise me about how to go about this?  The Diagnostics package is complex.
>
>
> Many thanks,
> Paul
>
>
>
>
> Dr Paul Burns MInstP FRMetS
> Research Fellow
> Geophysical and Astrophysical Fluid Dynamics Group
> Department of Mathematics
> Harrison Building
> University of Exeter
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=04%7C01%7CP.Burns2%40exeter.ac.uk%7Ce0d16f86b9f446d29cdf08d9e030d7ad%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C637787322844758474%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=UlMaeZBBwEBPNNxJXe31QVk9%2FrGERTToeMs5NsLDsGQ%3D&reserved=0<https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=04%7C01%7CP.Burns2%40exeter.ac.uk%7Cb45157bd9e88439653b208d9e162f18a%7C912a5d77fb984eeeaf321334d8f04a53%7C0%7C0%7C637788637597968271%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=68j81eFjo41MzNL3VPM9SeB882hhAqA0uYCiLDhpbpQ%3D&reserved=0>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20220328/7446b2f7/attachment-0001.html>