
<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Fatemeh,<div><br></div><div>have a look at the RBCS package, pkg/rbcs, <a href="https://mitgcm.readthedocs.io/en/latest/phys_pkgs/rbcs.html">https://mitgcm.readthedocs.io/en/latest/phys_pkgs/rbcs.html</a></div><div><br></div><div>With this package you can restore 3D fields with arbitrary 3D masks. Via the mask you can also use regionally varying restoring timescales, i.e. the mask can have values between 0 (no restoring) to 1 (strongest restoring)</div><div><br></div><div>Setting the restoring timescale to the timestep is almost like prescribing the restoring temperature.</div><div><br></div><div>Hope that helps,</div><div><br></div><div>M<br id="lineBreakAtBeginningOfMessage"><div><br><blockquote type="cite"><div>On 4. Jul 2025, at 14:55, Fateme Sharifi <fateme.sharifi1391@gmail.com> wrote:</div><br class="Apple-interchange-newline"><div><meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr"><div><p>Dear MITgcm Support Team,</p><p>I am currently working on an idealized model setup using MITgcm to simulate <strong>seasonal ice-covered lakes</strong>, with a focus on <strong>moat formation</strong>—ice-free lateral regions that remain at approximately 4 °C throughout the running period.</p><p>In my setup, I represent an idealized freshwater lake with a <strong>Gaussian-shaped bathymetry</strong>. The lake is initialized with a <strong>vertically stratified temperature profile</strong>, ranging from 0 °C at the surface to 2 °C at the bottom in the interior. To simulate moats, I initialize a lateral zone along the lake boundary with a uniform temperature of <strong>4 °C</strong>, representing maximum freshwater density.</p><p>My goal is to maintain this <strong>lateral 4 °C moat temperature</strong> fixed throughout the simulation period, while allowing the interior to evolve freely according to the model dynamics.</p><p>I would appreciate any guidance on the following:</p>

<ul>

<li><p>Is there a way to <strong>prescribe a constant temperature</strong> (e.g., via a boundary condition or restoring term) in specific lateral grid cells (e.g., within a ~500 m wide margin)?</p>

</li>

<li><p>Would you recommend using <strong>sponge layers</strong>, <strong>restoring zones</strong>, or another method for this purpose in the MITgcm framework?</p>

</li>

</ul><p>Thank you in advance for your support. </p><p>Best regards,</p></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Fatemeh Sadat Sharifi<br>Ph.D. Candidate at<br>The Leibniz Institute of Freshwater<div> Ecology and Inland Fisheries (IGB)<br></div></div></div></div></div>

_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></body></html>