Changeset 13959 for NEMO/branches/2020/SI3_martin_ponds/src/ICE/icethd.F90

Timestamp:

2020-12-01T23:47:44+01:00 (4 years ago)

Author:

clem

Message:

multiple adds refering to tickets #2573 #2575 #2576. It concerns small bugs corrections for having a perfectly conservative sea-ice system, plus the removal of restriction on snow layers (one can have several layers in the snow now), plus a bug fix in very peculiar situation where ocean is always supercool

File:

• NEMO/branches/2020/SI3_martin_ponds/src/ICE/icethd.F90 (modified) (1 diff)

NEMO/branches/2020/SI3_martin_ponds/src/ICE/icethd.F90

@@ -166 +166 @@
          qsb_ice_bot(ji,jj) = rswitch * MIN( qsb_ice_bot(ji,jj), - zqfr_neg * r1_Dt_ice / 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.