<div dir="ltr"><div>Hi Paul,</div><div><br></div><div>Thank you very much for the information!<br></div><div>I'm using a static ice shelf, so I think I will continue applying the digging method.<br></div><div><br></div><div>Kind regards,</div><div><br></div><div>Enrico<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 18 nov 2020 alle ore 20:33 Paul Holland - UKRI BAS <<a href="mailto:pahol@bas.ac.uk">pahol@bas.ac.uk</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Enrico<br>
<br>
Briefly, I always choose to 'dig' out the cavity in order to make sure that all of the floating ice has some ocean under it.  My logic is that we know the grounding line very well from satellite interferometry, but we do not know the seabed/ice bed very well at all, and so I always choose to sacrifice the seabed in order to preserve the grounding line.<br>
<br>
When doing the 'digging', we usually choose to do it in such a way that the velocity points on the C grid have two cells open in the vertical.  This would in theory make sure that all ocean columns can be connected to each other with a very under-resolved overturning circulation.  This gets a bit complicated with partial cells.  You can see in Kaitlin Naughten's github code you refer to exactly how she is handling this.  Basically, you have to over-dig the tracer point bathymetry in order to make sure that the velocity points are open, however you want them to be.<br>
<br>
The above logic of digging the seabed only applies to simulations with static ice.  If you had a coupled ice sheet model, it would be a bit strange to dig the seabed in different places as the model advances or retreats the ice.  In that case, we have adopted other strategies, including digging the ice instead of the seabed (i.e. shallowing the ice in order to preserve a modellable water column thickness) and other trickery with the streamice package.<br>
<br>
I hope that helps a bit...<br>
<br>
Cheers<br>
<br>
Paul<br>
<br>
-------------------<br>
Message: 2<br>
Date: Wed, 18 Nov 2020 06:55:51 -0800<br>
From: Dimitris Menemenlis <<a href="mailto:menemenlis@jpl.nasa.gov" target="_blank">menemenlis@jpl.nasa.gov</a>><br>
To: MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>>, Yoshihiro Nakayama<br>
<<a href="mailto:Yoshihiro.Nakayama@lowtem.hokudai.ac.jp" target="_blank">Yoshihiro.Nakayama@lowtem.hokudai.ac.jp</a>>, "Michael P. Schodlok"<br>
<<a href="mailto:Michael.P.Schodlok@jpl.nasa.gov" target="_blank">Michael.P.Schodlok@jpl.nasa.gov</a>><br>
Subject: Re: [MITgcm-support] [EXTERNAL] Treatment of Antarctic<br>
grounding zone regions<br>
Message-ID: <<a href="mailto:C8C7DE4E-7A33-48F0-8803-32F35C681DBF@jpl.nasa.gov" target="_blank">C8C7DE4E-7A33-48F0-8803-32F35C681DBF@jpl.nasa.gov</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Michael and Yoshi, do you have any advice to give for question below?<br>
What decision did you make for your hi-res simulations, e.g., for:<br>
<a href="https://www.nature.com/articles/s41598-019-53190-6" rel="noreferrer" target="_blank">https://www.nature.com/articles/s41598-019-53190-6</a><br>
<br>
Thanks, Dimitris<br>
<br>
> On Nov 18, 2020, at 6:47 AM, Pochini, Enrico <<a href="mailto:epochini@inogs.it" target="_blank">epochini@inogs.it</a>> wrote:<br>
><br>
> Dear community,<br>
><br>
> I have a question for the polar oceanographers:<br>
><br>
> I have noticed that in bathymetric reconstructions available for Antarctica (e.g. Bedmap2, RTopo-2) there are regions corresponding to the grounding zone of the biggest ice shelves where the water column thickness (resulting from the subtraction draft - bathy) equates exactly 1m.<br>
><br>
> For my simulations I work in the Ross Sea at 5km resolution and vertical layers of a few tens of metres at the depth of the grounding zone. I am wondering whether this region should be cut out directly as non-resolvable in my vertical discretization, or whether, in a less conservative approach, the bathymetry should be deepened in order to resolve the region in at least one layer.<br>
><br>
> Employing a preprocessing algorithm (<a href="https://github.com/knaughten/mitgcm_python" rel="noreferrer" target="_blank">https://github.com/knaughten/mitgcm_python</a> <<a href="https://urldefense.us/v3/__https://github.com/knaughten/mitgcm_python__;!!PvBDto6Hs4WbVuu7!bWCWtDMcJoVTIrR5IAjyeN1eGS1eXJVRuH3eFS2l_JFJ3IhWUV5BBvLC7nnfv5WjG2XalrUIZ7A$" rel="noreferrer" target="_blank">https://urldefense.us/v3/__https://github.com/knaughten/mitgcm_python__;!!PvBDto6Hs4WbVuu7!bWCWtDMcJoVTIrR5IAjyeN1eGS1eXJVRuH3eFS2l_JFJ3IhWUV5BBvLC7nnfv5WjG2XalrUIZ7A$</a>>) I obtain a grounding-zone water-thickness of ~ 50-80 m, enough to fit ~1-2 layers.<br>
><br>
> I am uncertain about which way to take, so I'm curious to hear how you treat such thin cavities in high resolution models in other Antarctic seas, whether you cut them or keep them, and how you make them resolvable without affecting too much the overall setup.<br>
><br>
> Thanks!<br>
><br>
> Enrico P.<br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/" rel="noreferrer" target="_blank">https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/</a><br>
> mitgcm-support__;!!PvBDto6Hs4WbVuu7!bWCWtDMcJoVTIrR5IAjyeN1eGS1eXJVRuH<br>
> 3eFS2l_JFJ3IhWUV5BBvLC7nnfv5WjG2XakE4dloY$<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. UK Research and Innovation (UKRI) 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. UKRI does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses. Opinions, conclusions or other information in this message and attachments that are not related directly to UKRI business are solely those of the author and do not represent the views of UKRI.<br>
_______________________________________________<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>