<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>