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

[Martin Losch]

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20220127/d283aea7/attachment.html>