Changeset 15573

Timestamp:

2021-12-03T13:16:28+01:00 (2 years ago)

Author:

edblockley

Message:

Back-porting 20m supercool ice thicknes cutoff

File:

• NEMO/branches/UKMO/NEMO_4.0.4_supercool_cutoff/src/ICE/icethd.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.4_supercool_cutoff/src/ICE/icethd.F90

@@ -171 +171 @@
             qsb_ice_bot(ji,jj) = rswitch * MIN( qsb_ice_bot(ji,jj), - zqfr_neg * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) )
 
+            ! If conditions are always supercooled (such as at the mouth of ice-shelves), then ice grows continuously
+            ! ==> stop ice formation by artificially setting up the turbulent fluxes to 0 when volume > 20m (arbitrary)
+            IF( ( t_bo(ji,jj) - ( sst_m(ji,jj) + rt0 ) ) > 0._wp .AND. vt_i(ji,jj) >= 20._wp ) THEN
+               zqfr               = 0._wp
+               zqfr_pos           = 0._wp
+               qsb_ice_bot(ji,jj) = 0._wp
+            ENDIF
+
             ! --- Energy Budget of the leads (qlead, J.m-2) --- !
             !     qlead is the energy received from the atm. in the leads.