[MITgcm-support] How to add a new viscosity scheme in MITgcm?

[Martin Losch]

Yanhong,

the 3D viscosity coefficients (momentum equations) are computed in model/src/calc_viscosity.F and you can see how the more complicated packages do that, in case you want to do something similar, but you can just modify kappaRU/RV directly.

In analogy, the 3D diffusivity coefficients (tracer equations) are computed in calc_3d_diffusivity.F

Martin

> On 31. Dec 2021, at 06:27, 赖燕红 <yhlai at pku.edu.cn> wrote:
> 
> Hi Martin,
> 
> I want to add a new viscosity scheme in MITgcm, but I'm not clearly which files or how many files in models/ I need to modify. I want to make the vertical viscosity coefficient be dependent on temperature, like it is equal to the input parameter viscAr when the temperature is relatively higher, but varies with temperature following a kind of relationship when the temperature is low. And this scheme is to distinguish the solid phase and molten phase of silicates when simulating the magma ocean. 
> 
> I know that I need to modify the file models/src/calc_viscosity.F, but maybe there are other files I need to modify at the same time, do you know that? 
> 
> 
> 
> Thanks in advance!
> 
> 
> 
> Yanhong 
> 
> _______________________________________________
> 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/20220103/57644229/attachment.p7s>