<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Hi Martin,</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thank you so much; that's clarified a lot of things for me. Can I ask a few follow-up questions? 3 small clarifying ones, and 1 further advice one (much appreciated).</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">(1) I misunderstood that data.cal's "startDate1" corresponded to nIter0=0, instead of
just nIter0. Can I just check that that's correct? See the documentation that says: "</span><span style="letter-spacing: normal; font-family: Lato, proxima-nova, "Helvetica Neue", Arial, sans-serif; font-size: 16px; font-weight: 400; color: rgb(64, 64, 64); background-color: rgb(252, 252, 252);">Note
that the first record read and used by the EXF package corresponds to the value ’startDate1’ set in data.cal.</span><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">"
(<a href="https://mitgcm.readthedocs.io/en/latest/phys_pkgs/exf.html#compile-time-options" id="LPlnk837073" class="OWAAutoLink">https://mitgcm.readthedocs.io/en/latest/phys_pkgs/exf.html#compile-time-options</a>)</span></div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">(2) I also missed that MITgcm interpolated between records provided in T_air.bin! Would
that maybe explain why I was getting an error "lib-4016 : UNRECOVERABLE library error. A READ operation tried to read a nonexistent record (1002)." in a 1000yr long run? (However, my T_air.bin file in that case was 1001 records long, so I'm still not sure
what the problem is.</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">(3) Also didn't know about "repeatPeriod" in &EXF_NML_01 of data.exf - thank you!! This
will likely do exactly what I want to do. So if I set this to 31104000000.0 (1000yr in seconds) and set atempperiod = 31104000.0 (360 days in seconds) and have all other exf periods set as 0, does MITgcm know that that repeatPeriod parameter should only be
applied to T_air.bin and won't throw an error? (And in this scenario, the T_air input file doesn't need an additional entry like you mentioned, correct?)</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">(4) Do you have any advice for running a simulation with time-varying T_air whose period
is greater than 1000yr, given that I run simulations that last 1000 years, one after another, continuing on from pickup files. Say I want T_air with a periodicity of 10,000yr. And my simulation will be something like 20 or 40,000 years in total. I could make
10 separate T_air files to represent years 1-1000, 1001-2000, 2001-3000, and so on to 10,000 when the temp reaches its starting point again, and then just be very careful to just use the correct T_air file with each simulation. But there must be an easier
way than this? (However your method of making use of the "repeatPeriod" parameter seems perfect for my experiments where T_air varies with periods of 500yr and 1000yr.)</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">Thank you so much again, can't say how appreciative I am.</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">M. Shankle</span></div>
<div class="elementToProof"><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);"><br>
</span></div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div class="elementToProof" style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="appendonsend"></div>
<hr style="display: inline-block; width: 98%;">
<div id="divRplyFwdMsg" dir="ltr"><span style="font-family: Calibri, sans-serif; font-size: 11pt; color: rgb(0, 0, 0);"><b>From:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of mitgcm-support-request@mitgcm.org <mitgcm-support-request@mitgcm.org><br>
<b>Sent:</b> Friday, January 5, 2024 5:00 PM<br>
<b>To:</b> mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>Subject:</b> MITgcm-support Digest, Vol 247, Issue 1</span>
<div> </div>
</div>
<div><span style="font-family: Aptos, Aptos_EmbeddedFont, Aptos_MSFontService, Calibri, Helvetica, sans-serif; font-size: 11pt; color: rgb(0, 0, 0);">Send MITgcm-support mailing list submissions to<br>
mitgcm-support@mitgcm.org<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" id="OWAa6552110-1e16-956c-3aab-a05311973ecf" class="OWAAutoLink" data-auth="NotApplicable" data-loopstyle="linkonly"><font color="red"><b>MailScanner has detected a possible fraud attempt from "mailman.mitgcm.org" claiming to be</b></font>
https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C02%7Cmgs23%40st-andrews.ac.uk%7C4cd314c936914946bbd808dc0e0fd17c%7Cf85626cb0da849d3aa5864ef678ef01a%7C0%7C0%7C638400709577522679%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C62000%7C%7C%7C&sdata=P75cnnCxfiSNR%2FZouPo5ILUqeFVD4oQH6Q1lQe7labo%3D&reserved=0</a><br>
or, via email, send a message with subject or body 'help' to<br>
mitgcm-support-request@mitgcm.org<br>
<br>
You can reach the person managing the list at<br>
mitgcm-support-owner@mitgcm.org<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of MITgcm-support digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: Time-varying forcing and EXF package (Martin Losch)<br>
2. Re: MITgcm: Mixing in Vermix Experiment (Martin Losch)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 5 Jan 2024 10:05:41 +0100<br>
From: Martin Losch <Martin.Losch@awi.de><br>
To: MITgcm Support <mitgcm-support@mitgcm.org><br>
Subject: Re: [MITgcm-support] Time-varying forcing and EXF package<br>
Message-ID: <53EEB5FE-584D-4936-B328-9391AB856B7D@awi.de><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Madison,<br>
<br>
startDate_1/2 ind data.cal specify the point in time to which nIter0=0 refers (interally this is called baseTime), so with your data.cal, your year has 360days (TheCalendar=?model?), where nIter0=0 refers to Jan01,39500,0:00 (assuming that startdate_2=000000).
Alternatively you could set the startDate to 00010101 and for your pickup in year 39500 use an nIter0 that corresponds to this year, ie. 39500/deltaTclock<br>
<br>
All constant-in-time forcing fields need to have a ?period? of 0. (e.g. ?uwindperiod=0.?). I am not sure if you need a startDate for those, but it wouldn?t hurt.<br>
For your time varying forcing your setting says (again assuming atempstartdate2=000000), that the first record in the forcing file T_air.bin is on Jan01,39500,0:00, the second on Jan01,39501,0:00, etc, and the model interpolates bi-linearily between them. I
don?t know about the details of your experiment, but you may want to shift this start time to the middle of the year (and then you?d need an additional record for year0-1 for the bilinear interpolation, but see below for an easy way out).<br>
<br>
You can also set a repeat period, e.g. for your 1000 y file:<br>
<br>
repeatPeriod = 31104000000.0, ( 86400*360*1000 )<br>
<br>
in &EXF_NML_01 to make the model start again with the first record, if it goes beyond 1000ys. So in the ?scenario? where your nIter0=39500/deltaTclock, you could use atempstartdate1=00010101 with repeatPeriod = 31104000000.0,<br>
<br>
If your cycle is only 500y, then you could get away with a shorter forcing file and a repeatPeriod of 500y.<br>
<br>
Hope that helps,<br>
<br>
Martin<br>
<br>
> On 30. Dec 2023, at 20:15, Madison Shankle <mgs23@st-andrews.ac.uk> wrote:<br>
><br>
> Dear all,<br>
><br>
> Seeking help setting up time-varying forcing. I unfortunately can't submit jobs at the moment, so apologies for not doing more troubleshooting on my own, but at present it would be very helpful to just have my understanding of a few key parameters checked.<br>
><br>
> Background: I have a time-varying T_air.bin forcing file (dimensions nTime, nY, nX = 1000, 140, 72) where the time dimension consists of annual air temp profiles (so the entire file represents a 1000yr-long record, that varies in time). All other forcing
files are kept constant. T_air varies cyclically with a periodicity of 500yr, so two cycles are contained in T_air.bin. (I made the file 1000yrs long as I typically run my simulations for 1000yrs at a time, which takes about 24hr, and I can start the next
simulation the next day.)<br>
><br>
> Is my understanding of the following parameters correct, and are there any other parameters to be changed that I have missed out?<br>
><br>
> data.exf<br>
> --> atempstartdate1 = 395000101<br>
> --> atempperiod = 31104000<br>
> Here "atempstartdate1" specifies that I have read in a T_air.bin file whose record starts on Jan 1, year 39500. (I am starting my simulations off of a previous study's output, which ran ~32000yr, from which it then took an additional ~7000yr to spin-up my
simulations on ARCHER2 with some slight edits to the model configuration. So, with my model configuration properly spun up, I'm now ready to start the experimental (time-varying temp) simulation at year 39500).<br>
><br>
> data.cal<br>
> --> TheCalendar = 'model'<br>
> --> startDate1 = 395000101<br>
> Here "startDate1" tells the model which record to read first from T_air.bin. At this point, I can say that my 1000yr-long T_air record starts at 395000 (in data.exf) and this parameter tells the model to start reading from that first time step. So, I was
planning on just updating these two parameters each time I start a new 1000yr-long run (e.g. to 40500, 41500, 42500, etc. each time I start a new run). However, I suppose I could also simply write a 40,000yr-long T_air file (40000,ny,nx) (as my goal is to
run my simulation this long), in which case I could leave atempstartdate1 as 395000101, and only have to update data.cal's startDate1 parameter with each new model run.<br>
><br>
><br>
> Is that all correct and am I missing anything else that would need to be changed to set up a time-varying run? Apologies again; I nearly got this working before the holidays but as I said can't submit jobs at the moment, so making some progress by getting
feedback from this group while I sort that out would be very appreciated.<br>
><br>
> Thanks and best wishes for the new year,<br>
> M. Shankle<br>
><br>
><br>
><br>
><br>
><br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" id="OWA4224d202-e74c-61be-1731-ab7d68349249" class="OWAAutoLink" data-auth="NotApplicable" data-loopstyle="linkonly"><font color="red"><b>MailScanner has detected a possible fraud attempt from "mailman.mitgcm.org" claiming to be</b></font>
https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C02%7Cmgs23%40st-andrews.ac.uk%7C4cd314c936914946bbd808dc0e0fd17c%7Cf85626cb0da849d3aa5864ef678ef01a%7C0%7C0%7C638400709577522679%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C62000%7C%7C%7C&sdata=P75cnnCxfiSNR%2FZouPo5ILUqeFVD4oQH6Q1lQe7labo%3D&reserved=0</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 5 Jan 2024 11:20:34 +0100<br>
From: Martin Losch <Martin.Losch@awi.de><br>
To: MITgcm Support <mitgcm-support@mitgcm.org><br>
Subject: Re: [MITgcm-support] MITgcm: Mixing in Vermix Experiment<br>
Message-ID: <AA0AB524-7741-4E47-A198-139023A860C8@awi.de><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Martin,<br>
<br>
currently ?RiLimit? is not a namelist of pkg/pp81. You can just add it in pp81_readparms.F<br>
<br>
Similarly you could add a RiMax to the namelist defined in my82_readparms.F<br>
<br>
In either case I do not know if modifying these parameters makes physical sense, as they are hard coded in the original implementations (I think it was SPEM a predecessor of ROMS) from which the MITgcm versions are derived.<br>
<br>
Martin<br>
<br>
> On 4. Dec 2023, at 18:17, Martin COUBARD <martin.coubard@imt-atlantique.net> wrote:<br>
><br>
> Hello MITGCM Users,<br>
> I am implementing the vermix experiment available on Github (<a href="https://github.com/MITgcm/MITgcm/tree/master/verification/vermix" id="OWA19b12c5a-5dc8-644b-c7b3-2bede0f66094" class="OWAAutoLink" data-auth="NotApplicable" data-loopstyle="linkonly"><font color="red"><b>MailScanner has detected a possible fraud attempt from "github.com" claiming to be</b></font> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FMITgcm%2FMITgcm%2Ftree%2Fmaster%2Fverification%2Fvermix&data=05%7C02%7Cmgs23%40st-andrews.ac.uk%7C4cd314c936914946bbd808dc0e0fd17c%7Cf85626cb0da849d3aa5864ef678ef01a%7C0%7C0%7C638400709577522679%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C62000%7C%7C%7C&sdata=m3zNzYYdup3%2FKK4WUTIOa8Ys3O2FeUPWVM1C8QjFW0A%3D&reserved=0</a>)
and <a href="https://mitgcm.readthedocs.io/en/latest/examples/examples.html" id="OWA961fb1f4-ca0b-210b-fe21-e7dea9ac1e25" class="OWAAutoLink" data-auth="NotApplicable" data-loopstyle="linkonly"><font color="red"><b>MailScanner has detected a possible fraud attempt from "mitgcm.readthedocs.io" claiming to be</b></font>
https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmitgcm.readthedocs.io%2Fen%2Flatest%2Fexamples%2Fexamples.html&data=05%7C02%7Cmgs23%40st-andrews.ac.uk%7C4cd314c936914946bbd808dc0e0fd17c%7Cf85626cb0da849d3aa5864ef678ef01a%7C0%7C0%7C638400709577522679%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C62000%7C%7C%7C&sdata=37rk6oM6U7lgpOdllx7N5uPFj4pZT5WERJxxxNWNox4%3D&reserved=0</a>.<br>
> Everything is running smoothly with the KPP option, and we do see a mixed layer forming over timescales of a few days. Furthermore, a parameter may be tuned to increase/decrease mixing (Rincr, which may be find in data.kpp), which helps tune the model to
our expectations.<br>
> <Capture d?e?cran 2023-12-04 a? 14.59.53.png><br>
><br>
> <Capture d?e?cran 2023-12-04 a? 15.13.24.png><br>
> <Capture d?e?cran 2023-12-04 a? 15.08.17.png> However, we have a problem when using the PP81 and MY82 packages, when it comes to manipulating mixing. I changed RiMax in data.my82 and RiLimit in data.pp81, but apparently nothing happens when I adjust these
values.<br>
> Do you have any suggestions on how to effectively change these parameters and hence control the mixing in PP81 and MY82 packages?<br>
> <Capture d?e?cran 2023-12-04 a? 15.09.50.png><br>
> If any of you can assist me with this, I would greatly appreciate it.<br>
><br>
> Thank you very much.<br>
><br>
> Martin Coubard<br>
> <br>
> --<br>
><br>
> Coubard Martin<br>
> FISE 2A<br>
> 06 68 14 10 95 Technop?le Brest-Iroise<br>
> CS 83818<br>
> 29238 Brest Cedex 3<br>
><br>
> Une ?cole de l'IMT<br>
><br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" id="OWA4ad4518b-2bee-0bc6-12fd-2f8ebf2bdc0b" class="OWAAutoLink" data-auth="NotApplicable" data-loopstyle="linkonly"><font color="red"><b>MailScanner has detected a possible fraud attempt from "mailman.mitgcm.org" claiming to be</b></font>
https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C02%7Cmgs23%40st-andrews.ac.uk%7C4cd314c936914946bbd808dc0e0fd17c%7Cf85626cb0da849d3aa5864ef678ef01a%7C0%7C0%7C638400709577522679%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C62000%7C%7C%7C&sdata=P75cnnCxfiSNR%2FZouPo5ILUqeFVD4oQH6Q1lQe7labo%3D&reserved=0</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
MITgcm-support@mitgcm.org<br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" id="OWAd8bd7a3a-7569-42f6-de75-2ec0a6275f51" class="OWAAutoLink" data-auth="NotApplicable" data-loopstyle="linkonly"><font color="red"><b>MailScanner has detected a possible fraud attempt from "mailman.mitgcm.org" claiming to be</b></font> https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C02%7Cmgs23%40st-andrews.ac.uk%7C4cd314c936914946bbd808dc0e0fd17c%7Cf85626cb0da849d3aa5864ef678ef01a%7C0%7C0%7C638400709577522679%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C62000%7C%7C%7C&sdata=P75cnnCxfiSNR%2FZouPo5ILUqeFVD4oQH6Q1lQe7labo%3D&reserved=0</a><br>
<br>
<br>
------------------------------<br>
<br>
End of MITgcm-support Digest, Vol 247, Issue 1<br>
**********************************************</span></div>
</body>
</html>