<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
Hello Gordan,
<div><br>
</div>
<div>question 1:</div>
<div><br>
</div>
<div>We use the penalty method to constrain the ice temperature below pressure melting point (Tpmp).</div>
<div>If the ice exceeds Tpmp, we add a penalty to the stiffness matrix and load vector:</div>
<div><br>
</div>
<div>kappa x T = kappa x Tpmp</div>
<div><br>
</div>
<div>where kappa = kmax x 10^penalty_factor is a large number.</div>
<div>This forces T to be equal to Tpmp.</div>
<div>We can then retrieve the melting rate by computing how much energy was necessary to constrain T:</div>
<div><br>
</div>
<div>kappa (T-Tpmp) = L / c Mb</div>
<div><br>
</div>
<div>where L is the latent heat, c is the heat capacity and Mb the melting rate.</div>
<div>The problem with this method is that the melting rate is VERY noisy and can reach values that are not physically realistic. This is due to the fact that no regularization is applied (Mb is not in H1 but L2, if you are familiar with Hilbert and Lebesgue
spaces).</div>
<div><br>
</div>
<div>One alternative is to use the enthalpy model, which does not use penalties. It is still under development but you could give it a try. Just set md.thermal.isenthalpy=1;</div>
<div>Let us know if it works better.</div>
<div><br>
</div>
<div>question 2:</div>
<div><br>
</div>
<div>FS is much more sensitive than SIA numerically, and requires to have a mesh that is refined enough, boundary conditions that are consistent etc.</div>
<div>What we do generally is we don't jump from SIA to FS. You should try SSA (MacAyeal) and HO (Blatter/Pattyn) first. Also, most of the time, HO does a great job and it is not necessary to use a FS model.</div>
<div><br>
</div>
<div>Do not hesitate if you have more questions</div>
<div><br>
</div>
<div>Regards,</div>
<div>Mathieu<br>
<div>
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
<div>--</div>
<div>Mathieu Morlighem</div>
<div><br>
</div>
<div>Assistant Project Scientist</div>
<div>Department of Earth System Science</div>
<div>University of California, Irvine</div>
<div>Croul Hall, Irvine, CA 92697-3100</div>
<div>(949) 824-1353 (UCI office)</div>
<div>(818) 354-4134 (JPL office)</div>
<div>(626) 429-5780 (cell)</div>
<div><a href="mailto:Mathieu.Morlighem@uci.edu">Mathieu.Morlighem@uci.edu</a></div>
</div>
</div>
<br>
<div>
<div>On Jan 21, 2013, at 10:52 AM, Gordan Stuhne <<a href="mailto:gordan@atmosp.physics.utoronto.ca">gordan@atmosp.physics.utoronto.ca</a>> wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">Hello,<br>
<br>
I am quite interested in your model, and have been working with it<br>
using some of my own grids and data relating to Greenland.<br>
My SIA (Hutter) results basically look rational for 100yr prognostic<br>
runs (using thermal penalty_lock=10 for stability), but I have a few<br>
questions. Firstly, the basal melting rates reported in the<br>
output files are very high (on the order of m/yr). This is<br>
presumably computed from latent heat fluxes where the<br>
basal ice becomes temperate, but the details don't seem<br>
to be documented in the code and I am not sure if I am<br>
interpreting it right.<br>
<br>
Secondly, I am having trouble trying to run the same case in<br>
FS mode. I've tried adjusting time-step, viscosity overshoot,<br>
and other parameters with no success. Is there some<br>
useful rule of thumb for preparing the initial conditions?<br>
<br>
Your insight would be greatly appreciated,<br>
Thanks in advance,<br>
Gordan Stuhne,<br>
Dept. of Physics,<br>
University of Toronto<br>
<br>
_______________________________________________<br>
issm-support mailing list<br>
<a href="mailto:issm-support@issm.ess.uci.edu">issm-support@issm.ess.uci.edu</a><br>
http://issm.ess.uci.edu/cgi-bin/mailman/listinfo/issm-support<br>
</blockquote>
</div>
<br>
</div>
</body>
</html>