Changeset 9397

Timestamp:

2018-03-14T14:27:22+01:00 (7 years ago)

Author:

clem

Message:

dev_CNRS_2017: change ice routines so that reproducibility issue is solved for the ice (but there is still repro issue between r8996 and r9012 after 268 time steps of orca2-lim3-pisces)

Location:

branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3

Files:

• iceitd.F90 (modified) (2 diffs)
• icethd_ent.F90 (modified) (1 diff)

branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/iceitd.F90

@@ -119 +119 @@
             ! Compute thickness change in each ice category
             DO ji = 1, npti
-               zdhice(ji,jl) = h_i_2d(ji,jl) - h_ib_2d(ji,jl)
+               IF( a_i_2d(ji,jl) >  epsi10) zdhice(ji,jl) = h_i_2d(ji,jl) - h_ib_2d(ji,jl)
             END DO
          END DO
@@ -690 +690 @@
          zalpha = 0.05_wp
          zhmax  = 3._wp * rn_himean
-         hi_max(0) = 0._wp
+         hi_max(0) = 0._wp 
          DO jl = 1, jpl
             znum = jpl * ( zhmax+1 )**zalpha

branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_ent.F90

@@ -103 +103 @@
       END DO
 
-      zeh_cum1(1:npti,0) = 0._wp 
+      zeh_cum1(1:npti,0:nlay_i) = 0._wp 
       ! new cumulative q*h => linear interpolation
       DO jk0 = 1, nlay_i+1