<p style="font-family:宋体, arial, Verdana, sans-serif;font-size:14px;white-space:normal;">
        Hi Martin,
</p>
<p style="font-family:宋体, arial, Verdana, sans-serif;font-size:14px;white-space:normal;">
        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. 
</p>
<p style="font-family:宋体, arial, Verdana, sans-serif;font-size:14px;white-space:normal;">
        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? 
</p>
<p style="font-family:宋体, arial, Verdana, sans-serif;font-size:14px;white-space:normal;">
        <br>
</p>
<p style="font-family:宋体, arial, Verdana, sans-serif;font-size:14px;white-space:normal;">
        Thanks in advance!
</p>
<p style="font-family:宋体, arial, Verdana, sans-serif;font-size:14px;white-space:normal;">
        <br>
</p>
<p style="font-family:宋体, arial, Verdana, sans-serif;font-size:14px;white-space:normal;">
        Yanhong 
</p>