<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:"\@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:11.0pt;
font-family:"Calibri",sans-serif;}
span.EmailStyle18
{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 Stanislav, <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I just want to add that I have encountered this problem and Martin’s solution helped. Sea ice piling up at the boundaries will happen, when the sea ice velocity near the boundaries does not match the boundary values. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>You can try<o:p></o:p></p><p class=MsoNormal>#define OBCS_UVICE_OLD<o:p></o:p></p><p class=MsoNormal>#define OBCS_SEAICE_AVOID_CONVERGENCE<o:p></o:p></p><p class=MsoNormal>in OBCS_OPTION.h<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The details of this implimentation are in obcs_adjust_uvice.F. The downside of this method is that sea ice velocity will not exactly be what has been previously prescribed. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Best,<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='font-size:12.0pt;color:black'>From: </span></b><span style='font-size:12.0pt;color:black'>MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Stanislav Martyanov <martyanov.sd@gmail.com><br><b>Reply-To: </b><mitgcm-support@mitgcm.org><br><b>Date: </b>Monday, May 18, 2020 at 11:39<br><b>To: </b><MITgcm-support@mitgcm.org><br><b>Subject: </b>[MITgcm-support] SEAICE, very thick ice<o:p></o:p></span></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>Hello everyone!</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>I get very thick sea ice at several locations in the high-res Kara Sea model domain (horizontal resolution is 500-1200 m). Interestingly, the general pattern of sea ice concentration and thickness is OK compared to MERCATOR ocean reanalysis, but in some regions the modeled HEFF is reported by the monitor to be about 15 m after 5 months of integration (starting from realistic initial conditions). I have attached the figure to clarify what I mean and where it happens.</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>I suspect the ridging process near the open boundaries to be involved in this situation, but I do not know what can be done with it.</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>Such HEFF also produce problems in vertical layers when useRealFreshWaterFlux is TRUE, but I checked that such too thick ice is also produced when useRealFreshWaterFlux is OFF, so it is not a cause (at least, not the main cause).</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>Unfortunately, there exist too many options to just iterate through all of them.</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>May be you can advise what can be done to mitigate such overestimated HEFF? My preliminary suggestions are as follows, but I am not sure about them:</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>- to switch off the uice and vice at the OBCS?</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>- to set SEAICE_no_slip=TRUE to reduce the sea ice velocities near the coasts to prevent enhanced ridging?</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>- to use sponge at the open boundaries?</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>- there was a CPP key OBCS_SEAICE_AVOID_CONVERGENCE, but it is only for OBCS_UVICE_OLD which, in turn, is undefined by default and even may be not implemented, as stated.</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>Any advice is very welcome!</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>Stanislav</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'> </span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif;color:black;background:#FCFCFC'>PS: I use a SEAICE package configuration which is very close to the default one. What I have changed is:</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>#define SEAICE_VARIABLE_SALINITY</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>SEAICEadvScheme = 33,</span><o:p></o:p></p><p class=MsoNormal style='line-height:115%'><span style='font-size:12.0pt;line-height:115%;font-family:"Times New Roman",serif'>SEAICE_multDim = 1,</span><o:p></o:p></p></div><p class=MsoNormal>_______________________________________________ MITgcm-support mailing list MITgcm-support@mitgcm.org http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <o:p></o:p></p></div></body></html>