Changeset 12811 for NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_pnd.F90

Timestamp:

2020-04-24T17:20:27+02:00 (4 years ago)

Author:

clem

Message:

debug a restartability issue. All sette tests passed now

File:

• NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_pnd.F90 (modified) (2 diffs)

NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_pnd.F90

@@ -184 +184 @@
             h_ip_1d(ji)      = 0._wp
             h_il_1d(ji)      = 0._wp
-            !
-            ! clem: problem with conservation or not ?
             !                                                         !--------------------------------!
          ELSE                                                         ! Case ice thickness >= rn_himin !
@@ -267 +265 @@
                zhp = ( h_i_1d(ji) * ( rau0 - rhoi ) + h_ip_1d(ji) * ( rau0 - rhow * a_ip_1d(ji) / a_i_1d(ji) ) ) * r1_rau0
 
-               ! Calculate the permeability of the ice (Assur 1958)
+               ! Calculate the permeability of the ice (Assur 1958, see Flocco 2010)
                DO jk = 1, nlay_i
                   zsbr = - 1.2_wp                                  &
                      &   - 21.8_wp    * ( t_i_1d(ji,jk) - rt0 )    &
                      &   - 0.919_wp   * ( t_i_1d(ji,jk) - rt0 )**2 &
-                     &   - 0.0178_wp  * ( t_i_1d(ji,jk) - rt0 )**3 ! clem: error here the factor was 0.01878 instead of 0.0178 (cf Flocco 2010)
+                     &   - 0.0178_wp  * ( t_i_1d(ji,jk) - rt0 )**3
                   ztmp(jk) = sz_i_1d(ji,jk) / zsbr
                END DO