<div dir="ltr">Hi Ou,<div><br></div><div>I forgot to mention in my previous post that there is no atmospheric forcing in all the three experiments. This is just a spindown problem of a mixed layer front.</div><div><br></div><div>Thanks,</div><div>Kalyan</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 13, 2020 at 9:52 AM Wang, Ou (US 329B) <<a href="mailto:ou.wang@jpl.nasa.gov">ou.wang@jpl.nasa.gov</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 lang="EN-US">
<div class="gmail-m_3434250195537659335WordSection1">
<p class="MsoNormal">Kalyan,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Do the three experiments have atmospheric forcing? If so, isn’t in Case A the ocean will feel the wind stress and other atmospheric forcing directly since there is no sea-ice (sea-ice concentration c=0), while Cases B and C don’t because
 the ocean is fully covered by sea-ice (c=1). <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I don’t have an explanation why Cases B and C have exactly similar results. They should be somewhat different.
<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Regards,<u></u><u></u></p>
<p class="MsoNormal">Ou<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(181,196,223);padding:3pt 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12pt;color:black">From: </span></b><span style="font-size:12pt;color:black">MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" target="_blank">mitgcm-support-bounces@mitgcm.org</a>> on behalf of Kalyan Shrestha <<a href="mailto:kalyansh@uw.edu" target="_blank">kalyansh@uw.edu</a>><br>
<b>Reply-To: </b>"<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>" <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>><br>
<b>Date: </b>Wednesday, May 13, 2020 at 9:28 AM<br>
<b>To: </b>"<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>" <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>><br>
<b>Subject: </b>[EXTERNAL] [MITgcm-support] Seaice water drag<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">Hello everyone, <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal">I am trying to explore baroclinic instability under sea ice. However, in the MITgcm sea ice model, I noticed that changing the drag doesn't seem to affect anything, which means that there is something weird happening with the sea ice-water
 drag. Let me elaborate this with an example - please refer the related figure in the dropbox link below. 
<u></u><u></u></p>
<div>
<p class="MsoNormal"><i><a href="https://www.dropbox.com/s/ipcnf2p95w00itp/SeaiceWaterDrag.pdf?dl=0" target="_blank">https://www.dropbox.com/s/ipcnf2p95w00itp/SeaiceWaterDrag.pdf?dl=0</a></i><u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">For this, I have designed three numerical experiments, which are as follows. The only difference in the following three cases is in sea ice concentration and sea ice water drag, while all other parameters are identical.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Case A-- sea ice concentration, c =0; and Sea ice-water drag, Cd = 5.35E-03<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Case B-- c = 1; Cd = 0<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Case C-- c = 1; Cd = 1.32E-03<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Physically, Case A and Case B should give the same result because both have zero sea ice-water stress. On the other hand, Case B and Case C should behave differently as their drag coefficients have changed. However, the results are surprising.
 Contrary to what is expected, Case A and Case B give unidentical results, whereas Case B and Case C give exactly similar results. Apart from total kinetic energy vs time plot and salinity contours in the attached figure, I also checked the energetics and other
 diagnostics, and the results do not make physical sense.  <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Could you please help me to understand this?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Kalyan Shrestha<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</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>