ESCOMP / CISM

Community Ice Sheet Model
GNU Lesser General Public License v3.0
6 stars 11 forks source link

Replacing value of minthck_ice_sheet in glissade_masks.F90 #69

Open mpetrini-norce opened 1 month ago

mpetrini-norce commented 1 month ago

minthck_ice_sheet is the thickness threshold (in meters) under which new ice sheet cells will not be initialized. A value of 2000 meters was originally selected to prevent the formation of ice caps in northern Greenland during the CESM2/CISM2 coupled spinup (https://doi.org/10.1029/2019MS001984). However, in the coupled 1pctCO2 simulation this value led to massive chunks of ice to be immediately discharged to the ocean as soon they would disconnect from the ice sheet. This resulted in unrealistic, very high spikes in the calving rate. To avoid this, a SourceMods change was introduced setting this constant to a safe value of 100 meters. To avoid the use of SourceMods, I suggest that this value (100 m) should be the standard value (at least for the CISM version included in NorESM2.5_alpha05 and subsequent releases.

https://github.com/ESCOMP/CISM/blob/3eeda3e85f98cfc137c6bdc968a1894a78af6096/libglissade/glissade_masks.F90#L628

@hgoelzer @mvertens