[MITgcm-support] Re: mitgcm seaice
ESTANISLAO-GAVILAN-PASCUAL-AHU
e.gavilan at hhu.edu.cn
Mon Nov 28 21:04:44 EST 2022
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
>**********************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20221129/263346f8/attachment.html>
More information about the MITgcm-support
mailing list