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

Timestamp:

2017-12-05T18:17:16+01:00 (7 years ago)

Author:

alexwestmohc

Message:

This branch version enables running of coupled NEMO-LIM at the Met Office for at
least ~200 timesteps (until killed by wallclock).

Fixes relative to previous version:

icethd_zdf -

Minimum snow, and ice thickness thresholds implemented, for calculation of
effective conductivity

sbccpl -

Sea ice coupling fields initialised to 0. in sbc_cpl_init, if ln_meto_cpl is
set. This may not be the right place to do this but I think these do need to
be set to zero before the initial coupling exchange and it's not clear that
this would happen otherwise (as LIM is not called till after this exchange

sbcmod -

Corrected an output message which reported the wrong logical

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

@@ -911 +911 @@
       ! effective conductivity and 1st layer temperature (Jules coupling)
       DO ji = 1, npti
-         cnd_ice_1d(ji) = 2._wp * ( isnow(ji) * zkappa_s(ji,0) + ( 1._wp - isnow(ji) ) * zkappa_i(ji,0) )
+         IF (h_s_1d(ji) > 0.1 ) THEN 
+             cnd_ice_1d(ji) = 2._wp * zkappa_s(ji,0)
+         ELSE
+             IF (h_i_1d(ji) > 0.1 ) THEN
+                 cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0)
+             ELSE
+                 cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) / 0.1
+             ENDIF
+         ENDIF
          t1_ice_1d (ji) =         ( isnow(ji) * t_s_1d  (ji,1) + ( 1._wp - isnow(ji) ) * t_i_1d  (ji,1) )
       END DO