<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><br><div>Hi Mike,</div><div><br></div><div>You are using an old version of the model. The new version has totally changed the code of that subroutine. Your problem could be related with the units of the water drag. In the older versions, the water drag is to be input in terms of density (e.g. about 5.5 in the namelist). So If you are testing normal values such as 0.005, the model will always meet that condition. </div><div><br></div><div>Kind regards</div><pre><div>
</div><div> Estanislao</div><div><br></div><div><br></div>>On 11/23/2022 9:02 AM, Michael A Spall wrote:
>> Hi All,
>> I am doing a problem with the MITgcm and the seaice package. The results are nearly insensitive to the value of SEAICE_waterDrag.
>> I see in the code that if a certain condition is met then the drag is set to 0.25, which I suspect is being met in my calculations (I am
>> checking this).  Can someone explain the reasoning behind this limit?
>>
>> from MITgcm_c67c, seaice_oceandrag_coeffs.F:
>>
>>             IF ( SEAICE_waterDrag.LE.0. ) THEN
>>              CwatC(I,J,bi,bj)=0.
>>             ELSEIF ( TEMPVAR.LE.(0.25 _d 0/SEAICE_waterDrag)**2 ) THEN
>>              CwatC(I,J,bi,bj)=0.25 _d 0
>>             ELSE
>>              CwatC(I,J,bi,bj)=SEAICE_waterDrag*SQRT(TEMPVAR)
>>             ENDIF
>>
>> Thanks,
>> Mike
>>
>> ==================================================================
>>
>> Michael A. Spall
>> Senior Scientist
>> Clark 323A
>> Department of Physical Oceanography
>> Woods Hole Oceanographic Institution
>> 360 Woods Hole Road MS #21
>> Woods Hole, MA 02543
>>
>> 508-289-3342
>> mspall@whoi.edu
>> https://urldefense.com/v3/__https://www2.whoi.edu/staff/mspall__;!!K-Hz7m0Vt54!nS_Dw_xQpzrotMZ2-NkZZ2MwbO4ort_7yY9DyiMC6iuck-GyACH2psYlBuH2mbOXELbNedtpD81SnnU$
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support@mitgcm.org
>> https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!K-Hz7m0Vt54!nS_Dw_xQpzrotMZ2-NkZZ2MwbO4ort_7yY9DyiMC6iuck-GyACH2psYlBuH2mbOXELbNedtp8fziL7Q$
>
>
>------------------------------
>
>Subject: Digest Footer
>
>_______________________________________________
>MITgcm-support mailing list
>MITgcm-support@mitgcm.org
>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
>------------------------------
>
>End of MITgcm-support Digest, Vol 233, Issue 8
>**********************************************
</pre></div><br>