<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"Times New Roman \(Body CS\)";
panose-1:2 2 6 3 5 4 5 2 3 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.EmailStyle18
{mso-style-type:personal;
font-family:"Calibri",sans-serif;
color:windowtext;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style></head><body lang=EN-US link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Hi, <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Let me update my previous question to make it more specific. For ISOMIP exp1, based on the model setup, outputs and Losch (2008) paper, I assume that the original setup of the boundary conditions in the verification/isomip/input is: (Correct me if anything is wrong)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Top (if there is open water area): Free surface, set up by rigidLid=.FALSE., implicitFreeSurface=.TRUE., in data file<o:p></o:p></p><p class=MsoNormal>Bottom: Free-slip, set up by no_slip_bottom=.FALSE., in data file<o:p></o:p></p><p class=MsoNormal>The lateral B.C.s: Closed walls and free-slip, set up by no_slip_sides=.FALSE., in data file<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>However, I notice that in gendata.m file, the south and west lateral boundaries have additional treatment as follows: <o:p></o:p></p><p class=MsoNormal><i>bathy = ones(nx,ny)*H; <o:p></o:p></i></p><p class=MsoNormal><i>bathy(1,:) = 0; <o:p></o:p></i></p><p class=MsoNormal><i>bathy(:,1) = 0;<o:p></o:p></i></p><p class=MsoNormal>Hence, the south and west boundaries are marked as land or solid boundaries by setting bathymetry=0, while the north and east do not have this treatment. <o:p></o:p></p><p class=MsoNormal>Initially, I thought this was unnecessary because the lateral boundaries are already closed and free-slip, but if I remove this treatment, the code blows up. <o:p></o:p></p><p class=MsoNormal>Therefore, I am wondering why we need this special treatment for the south and west B.C.s, which turns out to be very important. Meanwhile, appreciate it if anyone could help me to check my statement of B.C. setup. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Thanks,<o:p></o:p></p><p class=MsoNormal>Yilang<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span style='color:black'>From: </span></b><span style='color:black'>MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Yilang Xu <yxu@whoi.edu><br><b>Reply-To: </b><mitgcm-support@mitgcm.org><br><b>Date: </b>Thursday, July 26, 2018 at 11:40<br><b>To: </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br><b>Subject: </b>[MITgcm-support] ISOMIP default open water boundary condition<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:11.0pt'><o:p> </o:p></span></p></div><p class=MsoNormal>Hi, <o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>I try to add an open water area near the ice shelf front in SHELFICE package - ISOMIP exp1. I am wondering, if I just extend the domain further north to create the open water, without using OBCS package, what the default boundary condition would be for this open water area. <o:p></o:p></p><p class=MsoNormal>Moreover, for ISOMIP example in verification directory, what are the boundary conditions in the eastern and northern sides? They are not specified in gendata.m file. <o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal>Thanks,<o:p></o:p></p><p class=MsoNormal>Yilang<o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>_______________________________________________ MITgcm-support mailing list MITgcm-support@mitgcm.org http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <o:p></o:p></span></p></div></body></html>