Changeset 8905 for branches/UKMO

Timestamp:

2017-12-05T18:17:16+01:00 (6 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

Location:

branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO

Files:

• LIM_SRC_3/icethd_zdf.F90 (modified) (1 diff)
• OPA_SRC/SBC/sbccpl.F90 (modified) (1 diff)
• OPA_SRC/SBC/sbcmod.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

branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -260 +260 @@
       !
       CALL wrk_alloc( jpi,jpj,   zacs, zaos )
+
+      IF (ln_meto_cpl) THEN
+         tsfc_ice(:,:,:) = 0.0
+         a_ip(:,:,:) = 0.0
+         v_ip(:,:,:) = 0.0
+         t1_ice(:,:,:) = rt0
+         cnd_ice(:,:,:) = 0.0
+         sstfrz(:,:) = 0.0
+      ENDIF
 
       ! ================================ !

branches/UKMO/dev_r8183_ICEMODEL_svn_removed/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

@@ -137 +137 @@
          WRITE(numout,*) '         mixed forced-coupled     formulation       ln_mixcpl     = ', ln_mixcpl
 !!gm  lk_oasis is controlled by key_oasis3  ===>>>  It shoud be removed from the namelist 
-         WRITE(numout,*) '         Met Office coupling formulation            ln_mixcpl     = ', ln_meto_cpl
+         WRITE(numout,*) '         Met Office coupling formulation (sfc exchange in atmos)      ln_meto_cpl     = ', ln_meto_cpl
          WRITE(numout,*) '         OASIS coupling (with atm or sas)           lk_oasis      = ', lk_oasis
          WRITE(numout,*) '         components of your executable              nn_components = ', nn_components