[MITgcm-support] How to add a linear drag for U and V when temperature is lower than a value

[Martin Losch]

Hi Yanhong,

I think the routines apply_forcing_u/v can work because “THETA” is availalbe to them (via #include DYNVARS.h). So you should be able to just add this extra linear drag tendency to gU/V_arr at the end of each routine.

Martin

> On 23. Jun 2023, at 05:13, 赖燕红 <yhlai at pku.edu.cn> wrote:
> 
> Hello everyone!
> 
> I want to reduce u and v close to zero when the temperature is lower than a critical value 
> 
> (such as 2000 K, the melting temperature of silicates) by adding a linear drag in the equation in each layer.
> 
> 
> 
> I want to ask, how to realize it in model/src/apply_forcing.F, 
> 
> because the temperature is needed to determine whether to use a linear drag or not.
> 
> I'm not sure how to call temperature in this file.
> 
> 
> 
> Looking forward to your suggestions！ Thanks in advance.
> 
> Best,
> 
> Yanhong  
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230623/c5d901da/attachment.html>