[MITgcm-support] mitgcm seaice

[Michael A Spall]

Hi Estanislao,
Thanks for the suggestion but I have been using 5.5, so accounting for the units correctly.
I have downloaded the current version of the code anyway because the old one seems to crash
at random times, even after successfully running for a long time. If I restart it will continue to 
run beyond the time where it previously crashed. Hoping the new code is more stable.

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 at whoi.edu
https://www2.whoi.edu/staff/mspall

________________________________________
From: MITgcm-support <mitgcm-support-bounces at mitgcm.org> on behalf of ESTANISLAO-GAVILAN-PASCUAL-AHU <e.gavilan at hhu.edu.cn>
Sent: Monday, November 28, 2022 9:04 PM
To: mitgcm-support at mitgcm.org
Subject: [MITgcm-support]  Re: mitgcm seaice


Hi Mike,

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.

Kind regards



 Estanislao


>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 at 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 at 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 at mitgcm.org
>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
>------------------------------
>
>End of MITgcm-support Digest, Vol 233, Issue 8
>**********************************************