Changeset 10994 for NEMO/trunk/src/ICE/icecor.F90

Timestamp:

2019-05-17T15:08:34+02:00 (5 years ago)

Author:

clem

Message:

couple of commits to 1) deal with ice concentration reaching 1 (so, no more limitations in the max ice concentration imposed in the namelist), 2) deal with very small negative values that can occur due to roundoff errors

File:

• NEMO/trunk/src/ICE/icecor.F90 (modified) (2 diffs)

NEMO/trunk/src/ICE/icecor.F90

@@ -66 +66 @@
          WRITE(numout,*) '~~~~~~~'
       ENDIF
-      !
       !                             !-----------------------------------------------------
-      !                             !  ice thickness must exceed himin (for ice diff)    !
+      !                             !  ice thickness must exceed himin (for temp. diff.) !
       !                             !-----------------------------------------------------
       WHERE( a_i(:,:,:) >= epsi20 )   ;   h_i(:,:,:) = v_i(:,:,:) / a_i(:,:,:)
@@ -97 +96 @@
          END DO
       ENDIF
-
       !                             !-----------------------------------------------------
       !                             !  Rebin categories with thickness out of bounds     !