<div dir="ltr"><div><p>Dear MITgcm Community,</p>
<p>I am currently working on a MITgcm model of a lake with a Gaussian-shaped bathymetry and a 3D lake mask. I have encountered an issue where some wet points still have temperature (T) and salinity (S) values equal to zero.</p>
<p>In my setup, I have defined an RBCS mask overlaying the Gaussian-shaped shallow region of the lake. The RBCS temperature files currently have zero values outside the RBCS mask, and only the shallow RBCS areas are set to 4 °C. All wet points outside the RBCS region are left free to evolve, ensuring that the interior of the lake is not artificially constrained.</p>
<p>I would like to clarify whether it is necessary to define the RBCS mask over the entire lake, or if it is sufficient—and indeed better—to define it only in the shallow regions where restoring is desired (i.e., the moat/top slab) to maintain 4 °C.</p>
<p>Any guidance or recommendations on best practices for applying RBCS in such a lake setup would be greatly appreciated.</p>
<p>Thank you very much for your time and help.</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><br></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Fri, Jul 4, 2025 at 2:55 PM Fateme Sharifi <<a href="mailto:fateme.sharifi1391@gmail.com">fateme.sharifi1391@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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"><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>
</blockquote></div>