<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body>
<div dir="ltr">Hi Sara
<div><br>
</div>
<div>Thanks for the interest and sounds like a very interesting project -- i don't have anything to add regarding the calculation of tendency or vertical grid resolution, thanks for that Ken.</div>
<div><br>
</div>
<div>Just to add that if you are interested in resolving vertical movement of the ice-ocean interface on the scale of 10s to hundreds of meters, shelfice remeshing  would be a good tool to use, though if you are just interested in melt alone and a static ice-ocean
 interface there is no need. this paper (<a href="https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2017JC013251">https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2017JC013251</a>) is a more detailed reference for remeshing than the docs (and
 note that you do not need to use a dynamic ice-shelf model to make use of it). If you have any queries about its use please contact me. </div>
<div><br>
</div>
<div>Best</div>
<div>Dan</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Wed, Jun 2, 2021 at 10:26 PM Ken Hughes <<a href="mailto:kenneth.hughes@oregonstate.edu">kenneth.hughes@oregonstate.edu</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>
<p>Hi Sara,</p>
<p><br>
</p>
<p>With the conventional approach, melting/freezing doesn't really occur. In each horizontal grid cell, the ice thickness remains constant. A melt rate is calculated, but only so that the correct amount of heat and salt are added into the ocean beneath the
 ice. "Melting", for example, removes heat from grid cell below the ice shelf and effectively removes salt (essentially equivalent to diluting with freshwater).
<br>
</p>
<p><br>
</p>
<p>This addition/removal of heat/salt is what I think you're getting at when mentioning the tendency. If you know the temperature and salinity in the cell vertically adjacent to the ice, together with its velocity (to parameterize turbulent transfers), then
 you have everything you need to calculate the boundary. It can seem like a chicken–egg problem, but there are standard ways to deal with the heat flux into the ice to ensure the problem is solvable, e.g., the so-called three-equation formulation.<br>
</p>
<p><br>
</p>
<p>You mention you want a dynamic ice shell, so you might need to use the remeshing add on for the SHELFICE package. I haven't personally used it, but it is worth a look.</p>
<p><br>
</p>
<p>As for why we don't use a dense grid near the boundary, that's because the MITgcm works on a fixed vertical grid. Unless the ice shelf draft is everywhere the same, you can't always have a highly resolved grid below the ice.</p>
<p><br>
</p>
<p>If you haven't done so already, work through the Losch (2008) paper, which describes the physics behind the package.<br>
</p>
<p><br>
</p>
<p>Ken</p>
<p><br>
</p>
<div>On 2021-06-02 1:20 p.m., Miller, Sara G wrote:<br>
</div>
<blockquote type="cite">
<p><span style="color:rgb(215,63,9)">[This email originated from outside of OSU. Use caution with links and attachments.]</span></p>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hello all,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I am trying to understand the treatment of the ice-ocean boundary layer in the SHELFICE package. My goal is to use MITgcm to model a dynamic, global ice shell over an ocean on Jupiter's moon Europa. I believe the SHELFICE package will be useful for me, but
 I do not fully understand the description given in the user manual regarding boundary layer modeling. </div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">What is the total tendency and how is it used in the calculation of the boundary layer temperature at a future time step? The boundary layer temperature and freezing/melting
 information feels like a chicken-and-egg problem. It would make sense to me that you need to know the temperature in the BL to determine if you're in the freezing or melting regime, but MITgcm shows boundary layer temperature as a function of that total tendency
 term which I think captures the freezing/melting information. </span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">More fundamentally, how do you melt the ice in MITgcm? Does it slowly melt until it passes the threshold to step up to the next grid, or does its pressure just change?
 That is, the ice base is always the first grid cell, and as it melts upwards the pressure term just changes with it. I am new to modeling, but </span><span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">I don’t understand
 why we’re working in fractions of a grid cell instead of just making the grid more densely-spaced in the boundary layer region.</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">I apologize for the volume of questions in a single email! I appreciate any and all wisdom/advice.</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt"><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Best regards,</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt">Sara</span></div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
---</div>
<div id="gmail-m_-419345353792182637Signature">
<div>
<div id="gmail-m_-419345353792182637divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255);font-family:Calibri,Arial,Helvetica,sans-serif">
Sara G. Miller (she/her)
<div>PhD Student </div>
<div><span style="font-family:Calibri,Arial,Helvetica,sans-serif;background-color:rgb(255,255,255);display:inline">Georgia Institute of Technology</span><br>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset></fieldset>
<pre>_______________________________________________
MITgcm-support mailing list
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a>
</pre>
</blockquote>
</div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
-- <br>
<div dir="ltr" class="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div>--- PLEASE NOTE THAT I AM CURRENTLY WORKING FROM HOME AS A MEASURE OF SOCIAL DISTANCING DURING THE COVID-19 PANDEMIC ---</div>
<div dir="ltr"><br style="font-size:12.8px">
<span style="font-size:12.8px">Daniel Goldberg, PhD</span><br style="font-size:12.8px">
<span style="font-size:12.8px">Reader in Glaciology</span>
<div style="font-size:12.8px">School of Geosciences, University of Edinburgh<br>
<span style="color:rgb(0,0,0);font-family:arial,helvetica,clean,sans-serif;line-height:15px">Geography Building, Drummond Street, Edinburgh EH8 9XP</span><br>
<br>
<br>
em: <a href="mailto:dan.goldberg@ed.ac.uk" target="_blank">dan.goldberg@ed.ac.uk</a><br>
web: <a href="https://www.geos.ed.ac.uk/homes/dgoldber" target="_blank">https://www.geos.ed.ac.uk/homes/dgoldber</a></div>
</div>
</div>
</div>
</div>
</div>
</div>
The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. Is e buidheann carthannais a th’ ann an Oilthigh Dhùn Èideann, clàraichte an Alba, àireamh clàraidh SC005336.
</body>
</html>