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

[Burns, Paul]

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/20220126/19ae4c97/attachment-0001.html>