Changeset 8771 for branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/icethd_zdf.F90

Timestamp:

2017-11-21T17:28:51+01:00 (6 years ago)

Author:

clem

Message:

add 1st layer temperature for Jules

File:

• branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/icethd_zdf.F90 (modified) (1 diff)

branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/LIM_SRC_3/icethd_zdf.F90

@@ -909 +909 @@
       ! 11) Jules coupling: reset inner snow and ice temperatures, update conduction fluxes
       !---------------------------------------------------------------------------------------
-      ! effective conductivity (10cm treshold to avoid cnd_ice to be too big)
+      ! effective conductivity and 1st layer temperature (Jules coupling)
       DO ji = 1, npti
-         IF( zh_s(ji) >= 0.1_wp ) THEN   ;   cnd_ice_1d(ji) = zkappa_s(ji,0) * 2._wp
-         ELSE                            ;   cnd_ice_1d(ji) = zkappa_i(ji,0) * 2._wp
-         ENDIF
+         cnd_ice_1d(ji) = 2._wp * ( isnow(ji) * zkappa_s(ji,0) + ( 1._wp - isnow(ji) ) * zkappa_i(ji,0) )
+         t1_ice_1d (ji) =         ( isnow(ji) * t_s_1d  (ji,1) + ( 1._wp - isnow(ji) ) * t_i_1d  (ji,1) )
       END DO
       !