<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 style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Thank you very much Martin. A constant period of 2629800 seems to do the job!</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
All the best,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
Kaitlin</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<div id="Signature">
<div></div>
<div></div>
<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top: 0px; margin-bottom: 0px;margin-top:0px; margin-bottom:0px; margin-top:0; margin-bottom:0">
</p>
<p style="margin-top: 0px; margin-bottom: 0px;margin-top:0px; margin-bottom:0px; font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols; font-size:16px">
</p>
<div><span style="font-size:11pt; color:rgb(117,123,128)">
<p style="margin-top: 0px; margin-bottom: 0px;margin:0cm 0cm 0.0001pt; font-size:11pt; font-family:Calibri,sans-serif; color:rgb(32,31,30); text-align:start; background-color:rgb(255,255,255)">
<b><span style="margin:0px; color:rgb(31,73,125)">Dr Kaitlin Naughten</span></b><b><span style="margin:0px; font-size:9pt; color:rgb(31,73,125)"><span> </span></span></b><span style="margin:0px; font-size:9pt; color:rgb(31,73,125)">|<b><span> </span></b></span><span style="margin:0px; color:rgb(31,73,125)">Ocean-Ice
Modeller | British Antarctic Survey</span></p>
<p style="margin-top: 0px; margin-bottom: 0px;margin:0cm 0cm 0.0001pt; font-size:11pt; font-family:Calibri,sans-serif; color:rgb(32,31,30); text-align:start; background-color:rgb(255,255,255)">
<span style="margin:0px; color:rgb(31,73,125)">High Cross, Madingley Road, Cambridge CB3 0ET</span></p>
<p style="margin-top: 0px; margin-bottom: 0px;margin:0cm 0cm 0.0001pt; font-size:11pt; font-family:Calibri,sans-serif; color:rgb(32,31,30); text-align:start; background-color:rgb(255,255,255)">
<span style="margin:0px; color:rgb(31,73,125)">Email:<span> </span></span><u><span style="margin:0px; color:rgb(68,114,196)"><a href="mailto:kaight@bas.ac.uk" target="_blank" rel="noopener noreferrer" title="mailto:kaight@bas.ac.uk" style="margin:0px; color:rgb(149,79,114); text-decoration:underline"><span style="color:rgb(12,100,192)">kaight@bas.ac.uk</span></a></span></u></p>
</span></div>
<p style="margin-top: 0px; margin-bottom: 0px;margin-top:0px; margin-bottom:0px">
</p>
</div>
</div>
<div id="appendonsend"></div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Martin Losch <Martin.Losch@awi.de><br>
<b>Sent:</b> 25 October 2019 12:28<br>
<b>To:</b> MITgcm Support <mitgcm-support@mitgcm.org><br>
<b>Subject:</b> Re: [MITgcm-support] Monthly averaged OBCS with useOBCSYearlyFields</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="PlainText">Forgot about your option 2: This does not work with yearly fields, i.e. the period = -12. always assumes a climatological year with 12 records that does not change with year. Hence the model expects the correct filename and does not add
the year, i.e. if you data.obcs says “UVEL_piControl.OBCS_N” together with *period=-12, then the model looks for a file with this name (“UVEL_piControl.OBCS_N”), but not “UVEL_piControl.OBCS_N_1979”<br>
<br>
Martin<br>
<br>
<br>
> On 25. Oct 2019, at 13:23, Martin Losch <Martin.Losch@awi.de> wrote:<br>
> <br>
> Hi Kaitlin,<br>
> <br>
> as far as I know the calendarDumps option in data.cal on refers to the “dumps”, so the output of the model. For the input, this functionality is not implemented, not even for the exf pack. (I am assuming here, that you do use pkg/exf)<br>
> <br>
> That means to the *Period will always be interpreted as such, i.e. your 30days, will lead to reading new obcs fields (or exf fields) every 30 days, which will lead to a phase error and lack of data at the end (the record determination happens in S/R EXF_GetFFieldRec,
which is the same for obcs and exf.<br>
> With a Gregorian calendar, I normally compute the perod as<br>
> <br>
> 0.25*(366 + 3*365)*86400/12 = 2629800.<br>
> <br>
> Obviously, this will only work with a regular cycle of leap years. If you want to be more accurate, you’ll have to divide the entire length of your run (in seconds) by the total number of months of your simulation. This will lead to a small inaccuracy in
loading forcing data for a specific month, but there will be no phase error per year.<br>
> <br>
> Martin<br>
> <br>
> <br>
>> On 25. Oct 2019, at 13:04, Christoph Voelker <christoph.voelker@awi.de> wrote:<br>
>> <br>
>> Dear Kaitlin,<br>
>> I think your first method should work, BUT you need to provide also the monthly fields for the adjacent years, so the model has something to interpolate. In your case, was there a file UVEL_piControl.OBCS_N_1980 present?<br>
>> Best regards, Christoph<br>
>> <br>
>> Am 25.10.19 um 11:37 schrieb Naughten, Kaitlin A.:<br>
>>> Hi everyone,<br>
>>> <br>
>>> I am trying to run with monthly-averaged OBCS which change every year (useOBCSYearlyFields=true) and a Gregorian calendar meaning the months are different lengths. I am not sure that the code allows me to do this. I have tried two different things:<br>
>>> • Set obcs*period to 2592000 (30 days) and hope that the code turns it into a "real" month, as it does for the diagnostics. This doesn't seem to work as the model dies near the end of the first year when wants a 13th record to interpolate to, and can't
find it.<br>
>>> • Set obcs*period to -12 as you do for a repeating monthly climatology. Now the code dies because it's constructing a weird file name which doesn't exist (seems to be either UVEL_piControl.OBCS_N or UVEL_piControl.OBCS_N.001.001.data when it should
just be UVEL_piControl.OBCS_N_1979). Looking through the code in obcs_exf_load, it seems that obcs*period=-12 isn't set up to work with useOBCSYearlyFields, because the variables year0 and year1 are never initialised.<br>
>>> What do you suggest I do? Is there a way to set a non-constant OBCS period so the months can be different lengths? Or, an easy way to add this functionality to the code?<br>
>>> <br>
>>> Many thanks,<br>
>>> Kaitlin Naughten<br>
>>> <br>
>>> Dr Kaitlin Naughten | Ocean-Ice Modeller | British Antarctic Survey<br>
>>> High Cross, Madingley Road, Cambridge CB3 0ET<br>
>>> Email: kaight@bas.ac.uk<br>
>>> <br>
>>> <br>
>>> This email and any attachments are intended solely for the use of the named recipients. If you are not the intended recipient you must not use, disclose, copy or distribute this email or any of its attachments and should notify the sender immediately and
delete this email from your system. <br>
>>> UK Research and Innovation has taken every reasonable precaution to minimise risk of this email or any attachments containing viruses or malware but the recipient should carry out its own virus and malware checks before opening the attachments. UK Research
and Innovation does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses.
<br>
>>> Opinions, conclusions or other information in this message and attachments that are not related directly to UK Research and Innovation business are solely those of the author and do not represent the views of UK Research and Innovation.<br>
>>> <br>
>>> <br>
>>> <br>
>>> _______________________________________________<br>
>>> MITgcm-support mailing list<br>
>>> <br>
>>> MITgcm-support@mitgcm.org<br>
>>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
>> -- <br>
>> Christoph Voelker<br>
>> Alfred Wegener Institute <br>
>> Helmholtz Centre for Polar and Marine Research<br>
>> Am Handelshafen 12<br>
>> 27570 Bremerhaven, Germany<br>
>> e: <br>
>> Christoph.Voelker@awi.de<br>
>> <br>
>> t: +49 471 4831 1848<br>
>> <br>
>> _______________________________________________<br>
>> MITgcm-support mailing list<br>
>> MITgcm-support@mitgcm.org<br>
>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
MITgcm-support@mitgcm.org<br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</div>
</span></font></div>
<br>
<br>
This email and any attachments are intended solely for the use of the named recipients. If you are not the intended recipient you must not use, disclose, copy or distribute this email or any of its attachments and should notify the sender immediately and delete
this email from your system. <br>
UK Research and Innovation has taken every reasonable precaution to minimise risk of this email or any attachments containing viruses or malware but the recipient should carry out its own virus and malware checks before opening the attachments. UK Research
and Innovation does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses.
<br>
Opinions, conclusions or other information in this message and attachments that are not related directly to UK Research and Innovation business are solely those of the author and do not represent the views of UK Research and Innovation.<br>
<br>
</body>
</html>