<div dir="ltr">I am running a case exactly the same as the "tutorial_global_oce_latlon", but used the data that made by myself. The input data is from the ECCO "<a href="ftp://ftp-icdc.cen.uni-hamburg.de/EASYInit/ECCO/">ftp://ftp-icdc.cen.uni-hamburg.de/EASYInit/ECCO/</a>", the bathymetry is ETOPO5, and the SST and SSS is from WOA94. Land-sea mask and topography adjustments have been applied by comparing the grid depth data and input data. <div><br></div><div>The model runs into a serious cold bias (the maximum SST drops rapidly to 25C soon and stays stable), though the temperature pattern seems Okay. Another weird thing is the the Eta has values on land! Did anybody face this issue before? Thank you for your help in advance.</div><div><br></div><div>Best,</div><div><br></div><div>James Young<br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 27, 2018 at 7:16 AM,  <span dir="ltr"><<a href="mailto:mitgcm-support-request@mitgcm.org" target="_blank">mitgcm-support-request@mitgcm.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Send MITgcm-support mailing list submissions to<br>
        <a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:mitgcm-support-request@mitgcm.org">mitgcm-support-request@mitgcm.<wbr>org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:mitgcm-support-owner@mitgcm.org">mitgcm-support-owner@mitgcm.<wbr>org</a><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: Controlling MNC output (Jody Klymak)<br>
   2. Re: SST relaxation with temporal variability (Uchida Takaya)<br>
   3. Re: Controlling MNC output (Gus Correa)<br>
   4. Re: SST relaxation with temporal variability (Martin Losch)<br>
   5. Blocking of light by sea ice and plankton (Andrew Twelves)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Mon, 26 Mar 2018 09:07:26 -0700<br>
From: Jody Klymak <<a href="mailto:jklymak@uvic.ca">jklymak@uvic.ca</a>><br>
To: <a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a><br>
Subject: Re: [MITgcm-support] Controlling MNC output<br>
Message-ID: <<a href="mailto:43B28540-9896-4E55-8E7E-DA2249E79A87@uvic.ca">43B28540-9896-4E55-8E7E-<wbr>DA2249E79A87@uvic.ca</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Taylor,<br>
<br>
If you just use the state-dump files, I don?t know of a way to only starts saving after N years, except to stop the run after saving a pickup, and then restart with the pickup.<br>
<br>
If you use the diagnostics package, then you can set the ?phase? of the save to be many years in the future and the saves won?t start until that point.<br>
<br>
Cheers,  Jody<br>
<br>
> On Mar 26, 2018, at  8:58 AM, Taylor Shropshire <<a href="mailto:tas14j@my.fsu.edu">tas14j@my.fsu.edu</a>> wrote:<br>
><br>
> Hello,<br>
><br>
> I am attempting to save computational resources by having the model write output files only after a predetermined amount of spin up. I see variables like deltaTClock and dumpFreq but can't seem to figure out how to make the model write output files only after a few years of a simulation have been ran.<br>
><br>
> Thanks,<br>
> Taylor<br>
><br>
> Taylor Shropshire<br>
> PhD Candidate - Oceanography<br>
> Center for Ocean - Atmospheric Prediction Studies<br>
> Florida State University<br>
> ______________________________<wbr>_________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a> <mailto:<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.<wbr>org</a>><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a> <<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a>><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20180326/1f3179ef/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>pipermail/mitgcm-support/<wbr>attachments/20180326/1f3179ef/<wbr>attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 26 Mar 2018 12:53:28 -0400<br>
From: Uchida Takaya <<a href="mailto:tu2140@columbia.edu">tu2140@columbia.edu</a>><br>
To: <a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a><br>
Subject: Re: [MITgcm-support] SST relaxation with temporal variability<br>
Message-ID: <<a href="mailto:6169FD67-220D-4B0E-8247-4D89C92A5239@columbia.edu">6169FD67-220D-4B0E-8247-<wbr>4D89C92A5239@columbia.edu</a>><br>
Content-Type: text/plain; charset=utf-8<br>
<br>
Hi Martin,<br>
<br>
Thank you very much for your reply and information.<br>
<br>
If I understand correctly what you have mentioned, as long as the forcing data (e.g. `zonalWindFile`, `thetaClimeFile` etc.) all have the same relaxation time scale, cycle and period, the time varying surface forcing should work fine with just prescribing the forcing files along with the `tauThetaClimRelax`, `externForcingCycle` and `externForcingPeriod` parameters in the namelist file `data` without having to resort to the `rbcs` nor `exf` package. Is this understanding correct?<br>
<br>
Best,<br>
Takaya<br>
<br>
> On Mar 26, 2018, at 4:37 AM, Martin Losch <<a href="mailto:Martin.Losch@awi.de">Martin.Losch@awi.de</a>> wrote:<br>
><br>
> Hi Takaya,<br>
><br>
> `externForcingCycle` and `externForcingPeriod` are used for all forcing data that is specified in the namelist file ?data?. That makes this way of reading and applying forcing data very unflexible.<br>
><br>
> You could use the rbcs package for surface restoring, although I would not recommend it as it?s a bit of an overkill. rbcs can be used for 3D restoring and involves a lot of extra storage.<br>
><br>
> The recommended surface forcing package is the exf-package. There you have different forcing periods etc for each forcing field, including surface restoring. The exf-package is more involved in terms of configuration, but in the long run it pays off to use it.<br>
><br>
> Martin<br>
><br>
>> On 20. Mar 2018, at 20:28, Uchida Takaya <<a href="mailto:tu2140@columbia.edu">tu2140@columbia.edu</a>> wrote:<br>
>><br>
>> Dear MITgcm support,<br>
>><br>
>><br>
>> I am trying to include sea-surface temperature (SST) relaxation with temporal variability.<br>
>><br>
>> I found the documentation for which parameters to change for temporally varying wind stress (<a href="http://mitgcm.org/public/r2_manual/latest/online_documents/node104.html" rel="noreferrer" target="_blank">http://mitgcm.org/public/r2_<wbr>manual/latest/online_<wbr>documents/node104.html</a>) but the `externForcingCycle` and `externForcingPeriod` parameters seem only to be incorporated into the momentum equation.<br>
>> Do the parameters also apply to the SST forcing prescribed by `thetaClimFile` in the `data` name file or would I need to use the `rbcs` package by prescribing `rbcsForcingPeriod`, `rbcsForcingCycle` and `tauRelaxT`?<br>
>> In the latter case, should I simply comment out `tauThetaClimRelax` and `thetaClimFile` in the `data` name file?<br>
>><br>
>> Any clarification on this would be helpful. Thank you very much for your support in advance.<br>
>><br>
>><br>
>> Best,<br>
>> Takaya<br>
>> ???????<br>
>> PhD Candidate<br>
>> Physical Oceanography<br>
>> Columbia University in the City of New York<br>
>> <a href="https://roxyboy.github.io/" rel="noreferrer" target="_blank">https://roxyboy.github.io/</a><br>
>><br>
>> ______________________________<wbr>_________________<br>
>> MITgcm-support mailing list<br>
>> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
><br>
> ______________________________<wbr>_________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 26 Mar 2018 13:18:15 -0400<br>
From: Gus Correa <<a href="mailto:gus@ldeo.columbia.edu">gus@ldeo.columbia.edu</a>><br>
To: MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>><br>
Subject: Re: [MITgcm-support] Controlling MNC output<br>
Message-ID: <<a href="mailto:82c9f30f-d1b4-1944-3f93-c24b19324ef0@ldeo.columbia.edu">82c9f30f-d1b4-1944-3f93-<wbr>c24b19324ef0@ldeo.columbia.edu</a><wbr>><br>
Content-Type: text/plain; charset=utf-8; format=flowed<br>
<br>
Hi Taylor<br>
<br>
Piggybacking on Jody's suggestions:<br>
<br>
1) My understanding is that deltaTClock (in seconds)<br>
controls the model time step<br>
(other than the momentum solver time step, controlled by deltaTMom),<br>
is the model fundamental heartbeat,<br>
and in most setups it boils down to the the same as deltaT,<br>
deltaTtracer, etc.<br>
(MITgcm insiders correct me, please correct me if I am wrong.)<br>
deltaTClock doesn't control the output frequency.<br>
<br>
2) You could run the spinup phase first,<br>
with a large value of dumpFreq (in seconds),<br>
say, greater or equal to the model endTime, to avoid output<br>
during spinup (if that is really what you want).<br>
[Having some output during spinup would let you QC the spinup.]<br>
The endTime or nTimeSteps in the data namelist file<br>
will tell the model when to stop.<br>
<br>
Then run again, now picking up at the end of the spinup phase,<br>
now with a smaller dumpFreq,<br>
to output what you need for post-spinup.<br>
<br>
You need to update nIter0 (in data)<br>
for this second phase of the run, to match the final<br>
iteration of the spinup phase.<br>
<br>
The same technique can be applied if you're using the<br>
diagnostics package, by changing the values in data.diagnostics<br>
between the two phases.<br>
BTW, the diagnostic package allows finer grained control.<br>
<br>
3) In my experience with the MITgcm,<br>
to avoid to have jobs running for extremely long periods of time<br>
(sometimes doing the wrong thing),<br>
and to promote more flexibility to make changes such<br>
as the one you mention,<br>
one can handle these things in at least two ways, both laborious:<br>
<br>
A) If you're running on a cluster with a job scheduler, you could<br>
write a wrapper job submission script<br>
(in bash, tcsh, perl, python, whatever) that:<br>
A1) updates nIter0 (in data) after each model run ends (after mpirun<br>
./mitgcmuv)<br>
A2) then resubmits the job (with "qsub $0",<br>
or equivalent command for the job scheduler you use).<br>
Most job schedulers are configured to allow jobs to resubmit themselves<br>
this way.<br>
<br>
B) If you're running on a workstation or standalone machine without<br>
a job scheduler, you could do the steps in item A) manually.<br>
<br>
This way you can gracefully break your run in two phases (spinup and<br>
"steady state"), and adjust the various data* namelists according to<br>
your needs.<br>
<br>
This also has the advantage of allowing you to QC partial results along<br>
the way, avoiding long and costly runs with mistaken setups.<br>
<br>
I hope this helps,<br>
Gus Correa<br>
<br>
<br>
On 03/26/2018 12:07 PM, Jody Klymak wrote:<br>
> Hi Taylor,<br>
><br>
> If you just use the state-dump files, I don?t know of a way to only<br>
> starts saving after N years, except to stop the run after saving a<br>
> pickup, and then restart with the pickup.<br>
><br>
> If you use the diagnostics package, then you can set the ?phase? of the<br>
> save to be many years in the future and the saves won?t start until that<br>
> point.<br>
><br>
> Cheers, ?Jody<br>
><br>
>> On Mar 26, 2018, at ?8:58 AM, Taylor Shropshire <<a href="mailto:tas14j@my.fsu.edu">tas14j@my.fsu.edu</a><br>
>> <mailto:<a href="mailto:tas14j@my.fsu.edu">tas14j@my.fsu.edu</a>>> wrote:<br>
>><br>
>> Hello,<br>
>><br>
>> I am attempting to save computational resources by having the model<br>
>> write output files only after a predetermined amount of spin up. I see<br>
>> variables like deltaTClock and dumpFreq?but can't seem to figure out<br>
>> how to make the model write output files only after a few years of a<br>
>> simulation have been ran.<br>
>><br>
>> Thanks,<br>
>> Taylor<br>
>><br>
>> *Taylor Shropshire<br>
>> PhD Candidate?- Oceanography*<br>
>> *Center for Ocean - Atmospheric Prediction Studies*<br>
>> *Florida State University*<br>
>> ______________________________<wbr>_________________<br>
>> MITgcm-support mailing list<br>
>> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a> <mailto:<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.<wbr>org</a>><br>
>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
><br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Tue, 27 Mar 2018 09:35:49 +0200<br>
From: Martin Losch <<a href="mailto:Martin.Losch@awi.de">Martin.Losch@awi.de</a>><br>
To: MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>><br>
Subject: Re: [MITgcm-support] SST relaxation with temporal variability<br>
Message-ID: <<a href="mailto:406CEEC4-FC4B-454C-8738-959C5042CF42@awi.de">406CEEC4-FC4B-454C-8738-<wbr>959C5042CF42@awi.de</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Yes,<br>
<br>
except that the relaxation time scale only affects the surface restoring fields (files thetaClimFile, saltClimFile), but I think you mean that.<br>
<br>
M.<br>
<br>
> On 26. Mar 2018, at 18:53, Uchida Takaya <<a href="mailto:tu2140@columbia.edu">tu2140@columbia.edu</a>> wrote:<br>
><br>
> Hi Martin,<br>
><br>
> Thank you very much for your reply and information.<br>
><br>
> If I understand correctly what you have mentioned, as long as the forcing data (e.g. `zonalWindFile`, `thetaClimeFile` etc.) all have the same relaxation time scale, cycle and period, the time varying surface forcing should work fine with just prescribing the forcing files along with the `tauThetaClimRelax`, `externForcingCycle` and `externForcingPeriod` parameters in the namelist file `data` without having to resort to the `rbcs` nor `exf` package. Is this understanding correct?<br>
><br>
> Best,<br>
> Takaya<br>
><br>
>> On Mar 26, 2018, at 4:37 AM, Martin Losch <<a href="mailto:Martin.Losch@awi.de">Martin.Losch@awi.de</a>> wrote:<br>
>><br>
>> Hi Takaya,<br>
>><br>
>> `externForcingCycle` and `externForcingPeriod` are used for all forcing data that is specified in the namelist file ?data?. That makes this way of reading and applying forcing data very unflexible.<br>
>><br>
>> You could use the rbcs package for surface restoring, although I would not recommend it as it?s a bit of an overkill. rbcs can be used for 3D restoring and involves a lot of extra storage.<br>
>><br>
>> The recommended surface forcing package is the exf-package. There you have different forcing periods etc for each forcing field, including surface restoring. The exf-package is more involved in terms of configuration, but in the long run it pays off to use it.<br>
>><br>
>> Martin<br>
>><br>
>>> On 20. Mar 2018, at 20:28, Uchida Takaya <<a href="mailto:tu2140@columbia.edu">tu2140@columbia.edu</a>> wrote:<br>
>>><br>
>>> Dear MITgcm support,<br>
>>><br>
>>><br>
>>> I am trying to include sea-surface temperature (SST) relaxation with temporal variability.<br>
>>><br>
>>> I found the documentation for which parameters to change for temporally varying wind stress (<a href="http://mitgcm.org/public/r2_manual/latest/online_documents/node104.html" rel="noreferrer" target="_blank">http://mitgcm.org/public/r2_<wbr>manual/latest/online_<wbr>documents/node104.html</a>) but the `externForcingCycle` and `externForcingPeriod` parameters seem only to be incorporated into the momentum equation.<br>
>>> Do the parameters also apply to the SST forcing prescribed by `thetaClimFile` in the `data` name file or would I need to use the `rbcs` package by prescribing `rbcsForcingPeriod`, `rbcsForcingCycle` and `tauRelaxT`?<br>
>>> In the latter case, should I simply comment out `tauThetaClimRelax` and `thetaClimFile` in the `data` name file?<br>
>>><br>
>>> Any clarification on this would be helpful. Thank you very much for your support in advance.<br>
>>><br>
>>><br>
>>> Best,<br>
>>> Takaya<br>
>>> ???????<br>
>>> PhD Candidate<br>
>>> Physical Oceanography<br>
>>> Columbia University in the City of New York<br>
>>> <a href="https://roxyboy.github.io/" rel="noreferrer" target="_blank">https://roxyboy.github.io/</a><br>
>>><br>
>>> ______________________________<wbr>_________________<br>
>>> MITgcm-support mailing list<br>
>>> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
>>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
>><br>
>> ______________________________<wbr>_________________<br>
>> MITgcm-support mailing list<br>
>> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Tue, 27 Mar 2018 11:16:40 +0000<br>
From: Andrew Twelves <<a href="mailto:agt_94@hotmail.co.uk">agt_94@hotmail.co.uk</a>><br>
To: "<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>" <<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>><br>
Subject: [MITgcm-support] Blocking of light by sea ice and plankton<br>
Message-ID:<br>
        <<a href="mailto:DB6PR0501MB266204BA9F0DD91D117BB072B2AC0@DB6PR0501MB2662.eurprd05.prod.outlook.com">DB6PR0501MB266204BA9F0DD91D11<wbr>7BB072B2AC0@DB6PR0501MB2662.<wbr>eurprd05.prod.outlook.com</a>><br>
<br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Morning all,<br>
<br>
<br>
<br>
I am about to start biogeochemical modelling using the BLING package, and have a couple of questions regarding the treatment of light availability:<br>
<br>
<br>
<br>
  1.  Does BLING by default set light penetration to zero underneath sea ice, or does it allow for limited penetration subject to a strong attenuation coefficient?  If the latter, does BLING take the value of effective sea ice thickness HEFF directly from the Sea Ice package, or does it do something more sophisticated (like a sub-grid scale parametrization of the thickness distribution)?<br>
  2.  Does BLING take the self-shading effect of phytoplankton themselves into account for light calculations?<br>
<br>
<br>
<br>
I appreciate that these are potentially broad questions but any guidance would be much appreciated.<br>
<br>
<br>
<br>
Many thanks in advance,<br>
<br>
<br>
<br>
Andrew Twelves<br>
<br>
University of Edinburgh<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20180327/ec6abe7a/attachment.html" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>pipermail/mitgcm-support/<wbr>attachments/20180327/ec6abe7a/<wbr>attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
______________________________<wbr>_________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-<wbr>support</a><br>
<br>
<br>
------------------------------<br>
<br>
End of MITgcm-support Digest, Vol 177, Issue 20<br>
******************************<wbr>*****************<br>
</blockquote></div><br></div></div></div>