Changeset 13169 for NEMO/branches

Timestamp:

2020-06-27T19:10:19+02:00 (4 years ago)

Author:

clem

Message:

adapt conductivities for met-office. Here a minimum of 1mm of snow is used, but we may need to increase it to 1cm (only for met-office specific fields)

File:

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

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

@@ -133 +133 @@
       REAL(wp), DIMENSION(jpij)            ::   za_s_fra    ! ice fraction covered by snow 
       REAL(wp), DIMENSION(jpij)            ::   isnow       ! snow presence (1) or not (0) 
+      REAL(wp), DIMENSION(jpij)            ::   isnow_comb  ! snow presence for met-office 
       !
       ! Mono-category
@@ -182 +183 @@
             isnow (ji) = 0._wp
          ENDIF
+         ! for Met-Office
+         IF( h_s_1d(ji) < zh_min ) THEN
+            isnow_comb(ji) = h_s_1d(ji) / zh_min
+         ELSE
+            isnow_comb(ji) = 1._wp
+         ENDIF
       END DO
       ! clem: we should apply correction on snow thickness to take into account snow fraction
@@ -331 +338 @@
          DO ji = 1, npti   ! Snow-ice interface
             IF ( .NOT. l_T_converged(ji) ) THEN
-               ! Calculate combined surface snow and ice conductivity to pass through the coupler
-               zkappa_comb(ji) = isnow(ji) * zkappa_s(ji,0) + ( 1._wp - isnow(ji) ) * zkappa_i(ji,0)
+               ! Calculate combined surface snow and ice conductivity to pass through the coupler (met-office)
+               zkappa_comb(ji) = isnow_comb(ji) * zkappa_s(ji,0) + ( 1._wp - isnow_comb(ji) ) * zkappa_i(ji,0)
                ! If there is snow then use the same snow-ice interface conductivity for the top layer of ice
                IF( h_s_1d(ji) > 0._wp )   zkappa_i(ji,0) = zkappa_s(ji,nlay_s)
@@ -896 +903 @@
       DO ji = 1, npti
          IF( h_i_1d(ji) >= zhi_ssl ) THEN
-!!clem metO            cnd_ice_1d(ji) = 2._wp * zkappa_comb(ji)
-            cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0)
+            cnd_ice_1d(ji) = 2._wp * zkappa_comb(ji)
+            !!cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0)
          ELSE
             cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) / zhi_ssl ! cnd_ice is capped by: cond_i/zhi_ssl