[MITgcm-support] mitgcm seaice

[Michael A Spall]

Hi Jinlun and Martin,
Thanks for the feedback. I thought it might have something to do with stability. I will play around with it a little
to see what happens. Monitoring one of my runs it looks like it does hover around 0.25.
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 Martin Losch <Martin.Losch at awi.de>
Sent: Monday, November 28, 2022 3:06 AM
To: MITgcm Support
Subject: Re: [MITgcm-support] mitgcm seaice

Hi Mike,

the non-zero drag used to be necessary to have a non-singular system matrix for the solver also where there there is no ice.

In a more recent version of the code I added a runtime parameter SEAICEdWatMin, which defaults to 0.25. This value can be set to zero, and the model should still run. After different (older) changes it is not longer required that there’s a minimum drag.

I am not sure about checkpoint 67c, but you can most likely reduce 0.25 to a much smaller value, if necessary.

Martin


> On 23. Nov 2022, at 21:01, Jinlun Zhang <jlzhang at uw.edu> wrote:
>
> Hi Mike,
>
> The 0.25 value is used to prevent water drag from being too small in order to avoid possible instability. Water drag is often higher than 0.25, so the insensitivity might not come from this. Otherwise, the ice V or ocean V might be a bit too small.
>
> Cheers,
>
> Jinlun
>
> 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://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__https%3A%2F%2Fwww2.whoi.edu%2Fstaff%2Fmspall__%3B!!K-Hz7m0Vt54!nS_Dw_xQpzrotMZ2-NkZZ2MwbO4ort_7yY9DyiMC6iuck-GyACH2psYlBuH2mbOXELbNedtpD81SnnU%24&data=05%7C01%7Cmspall%40whoi.edu%7Cdf8eba7e3dc94929e10608dad1179534%7Cd44c5cc6d18c46cc8abd4fdf5b6e5944%7C0%7C0%7C638052196475135571%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=ekska85vDow87oVMrIpduf2j3PidMa5Au1r9wCe72AA%3D&reserved=0
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Furldefense.com%2Fv3%2F__http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support__%3B!!K-Hz7m0Vt54!nS_Dw_xQpzrotMZ2-NkZZ2MwbO4ort_7yY9DyiMC6iuck-GyACH2psYlBuH2mbOXELbNedtp8fziL7Q%24&data=05%7C01%7Cmspall%40whoi.edu%7Cdf8eba7e3dc94929e10608dad1179534%7Cd44c5cc6d18c46cc8abd4fdf5b6e5944%7C0%7C0%7C638052196475135571%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=0cbo0pzdo3lP0sOvkr5cB3zO95W20T%2FUUp9ukoUp8x4%3D&reserved=0
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C01%7Cmspall%40whoi.edu%7Cdf8eba7e3dc94929e10608dad1179534%7Cd44c5cc6d18c46cc8abd4fdf5b6e5944%7C0%7C0%7C638052196475135571%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=mnFnuLAC%2Fak%2F2ETH3539%2Ffni7wQz2jsmnqoP59HeIM8%3D&reserved=0

_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org
https://nam02.safelinks.protection.outlook.com/?url=http%3A%2F%2Fmailman.mitgcm.org%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C01%7Cmspall%40whoi.edu%7Cdf8eba7e3dc94929e10608dad1179534%7Cd44c5cc6d18c46cc8abd4fdf5b6e5944%7C0%7C0%7C638052196475135571%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=mnFnuLAC%2Fak%2F2ETH3539%2Ffni7wQz2jsmnqoP59HeIM8%3D&reserved=0