Hi Martin,<div><br></div><div><div>I was a little unrealistic in changing tAlpha to get the density I wanted. At first, I just thought of it as a constant, and I did not think of it in relation to stability. Now the problem has been solved in other ways.</div><div><br></div><div>Thank you for your answer, sincerely!</div><div><br></div><div>-best，</div><div>Soong</div><br><br><blockquote name="replyContent" class="ReferenceQuote" style="padding-left:5px;margin-left:5px;border-left:#b6b6b6 2px solid;margin-right:0">-----原始邮件-----<br>
<b>发件人:</b><span id="rc_from">"Martin Losch" <Martin.Losch@awi.de></span><br>
<b>发送时间:</b><span id="rc_senttime">2021-12-23 21:28:10 (星期四)</span><br>
<b>收件人:</b> "MITgcm Support" <mitgcm-support@mitgcm.org><br>
<b>抄送:</b> <br>
<b>主题:</b> Re: [MITgcm-support] How to select the thermal expansion coefficient<br><br>Hi Soong,<br><br>my guess is that with a larger tAlpha you change the stability  (drho/dr = - rhoConst*talpha d(theta-Tref)/dz) which could make the model unstable.<br><br>Did you try with a shorter timestep deltaT?<br><br>Martin<br><br>> On 23. Dec 2021, at 13:24, 万松 <wansong@stu.ouc.edu.cn> wrote:<br>> <br>> Hi,<br>> <br>> I'm sorry to bother you with my question! I recently modified the thermal expansion coefficient 'tAlpha' to get the density I want when running the baroclinic model, but the puzzle is that if I use a higher value than the default value ( 2.0E-4 ), the model will crash, can anyone share the solution?<br>> <br>> -best,<br>> Soong<br>> _______________________________________________<br>> MITgcm-support mailing list<br>> MITgcm-support@mitgcm.org<br>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br><br></blockquote></div>