Changeset 9415 for branches/2017/dev_merge_2017/NEMOGCM/NEMO

Timestamp:

2018-03-21T16:51:02+01:00 (6 years ago)

Author:

clem

Message:

debug ice diffusion: avoid too small ice thicknesses by including a threshold (rn_himin) after ice dynamics. This should prevent the ice-ocean conduction flux from being abnormally strong

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icecor.F90 (modified) (1 diff)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icecor.F90

@@ -66 +66 @@
       ENDIF
       !
-      IF( kn == 2 ) THEN
-         !                          !-----------------------------------------------------
-         !                          !  thickness of the smallest category above himin    !
-         !                          !-----------------------------------------------------
-         WHERE( a_i(:,:,1) >= epsi20 )   ;   h_i(:,:,1) = v_i (:,:,1) / a_i(:,:,1)
-         ELSEWHERE                       ;   h_i(:,:,1) = 0._wp
-         END WHERE
-         WHERE( h_i(:,:,1) < rn_himin )      a_i(:,:,1) = a_i (:,:,1) * h_i(:,:,1) / rn_himin
-         !
-      ENDIF
+      !                             !-----------------------------------------------------
+      !                             !  ice thickness must exceed himin (for ice diff)    !
+      !                             !-----------------------------------------------------
+      WHERE( a_i(:,:,:) >= epsi20 )   ;   h_i(:,:,:) = v_i(:,:,:) / a_i(:,:,:)
+      ELSEWHERE                       ;   h_i(:,:,:) = 0._wp
+      END WHERE
+      WHERE( h_i(:,:,:) < rn_himin )      a_i(:,:,:) = a_i(:,:,:) * h_i(:,:,:) / rn_himin
+      !
       !                             !-----------------------------------------------------
       !                             !  ice concentration should not exceed amax          !