<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Yanhong,<div><br></div><div>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.</div><div><br></div><div>Martin<br><div><br><blockquote type="cite"><div>On 23. Jun 2023, at 05:13, 赖燕红 <yhlai@pku.edu.cn> wrote:</div><br class="Apple-interchange-newline"><div><meta http-equiv="Content-Type" content="text/html; charset=utf-8"><p>
Hello everyone!
</p><p>
I want to reduce u and v close to zero when the temperature is lower than a critical value
</p><p>
(such as 2000 K, the melting temperature of silicates) by adding a linear drag in the equation in each layer.
</p><p>
<br>
</p><p>
I want to ask, how to realize it in model/src/apply_forcing.F,
</p><p>
because the temperature is needed to determine whether to use a linear drag or not.
</p><p>
I'm not sure how to call temperature in this file.
</p><p>
<br>
</p><p>
Looking forward to your suggestions！ Thanks in advance.
</p><p>
Best,
</p><p>
Yanhong
</p>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></body></html>