[MITgcm-support] How to select the thermal expansion coefficient

[Martin Losch]

Hi Soong,

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.

Did you try with a shorter timestep deltaT?

Martin

> On 23. Dec 2021, at 13:24, 万松 <wansong at stu.ouc.edu.cn> wrote:
> 
> Hi,
> 
> 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?
> 
> -best,
> Soong
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1665 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20211223/867455d4/attachment.p7s>