Changeset 5202 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

Timestamp:

2015-04-07T17:40:16+02:00 (9 years ago)

Author:

clem

Message:

LIM3: important bug fix to avoid crashes. See ticket #1508

File:

• trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90 (modified) (3 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

@@ -89 +89 @@
       REAL(wp) :: zfric_u, zqld, zqfr
       REAL(wp) :: zvi_b, zsmv_b, zei_b, zfs_b, zfw_b, zft_b 
-      REAL(wp), PARAMETER :: zfric_umin = 0._wp        ! lower bound for the friction velocity (cice value=5.e-04)
-      REAL(wp), PARAMETER :: zch        = 0.0057_wp    ! heat transfer coefficient
+      REAL(wp), PARAMETER :: zfric_umin = 0._wp           ! lower bound for the friction velocity (cice value=5.e-04)
+      REAL(wp), PARAMETER :: zch        = 0.0057_wp       ! heat transfer coefficient
       !
       REAL(wp), POINTER, DIMENSION(:,:) ::  zqsr, zqns
@@ -353 +353 @@
       END DO
  
-      !------------------------
-      ! Ice natural aging              
-      !------------------------
-      oa_i(:,:,:) = oa_i(:,:,:) + a_i(:,:,:) * rdt_ice /rday
-
       !----------------------------------
       ! Change thickness to volume
@@ -365 +360 @@
       smv_i(:,:,:) = sm_i(:,:,:) * v_i(:,:,:)
 
+      ! update ice age (in case a_i changed, i.e. becomes 0 or lateral melting in monocat)
+      DO jl  = 1, jpl
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               rswitch = MAX( 0._wp , SIGN( 1._wp, a_i_b(ji,jj,jl) - epsi10 ) )
+               oa_i(ji,jj,jl) = rswitch * oa_i(ji,jj,jl) * a_i(ji,jj,jl) / MAX( a_i_b(ji,jj,jl), epsi10 )
+            END DO
+         END DO
+      END DO
+
       CALL lim_var_zapsmall
+
       !--------------------------------------------
       ! Diagnostic thermodynamic growth rates