Changeset 12062 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfhdiv.F90

Timestamp:

2019-12-05T10:57:46+01:00 (4 years ago)

Author:

mathiot

Message:

changes required by N.J. review

File:

• NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfhdiv.F90 (modified) (1 diff)

NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfhdiv.F90

@@ -51 +51 @@
          IF ( ln_isfcpl .AND. kt /= 0 ) THEN
             !
-            ! correct divergence only for the first time step
+            ! Dynamical stability at start up after change in under ice shelf cavity geometry is achieve by correcting the divergence.
+            ! This is achieved by applying a volume flux in order to keep the horizontal divergence after remapping 
+            ! the same as at the end of the latest time step. So correction need to be apply at nit000 (euler time step) and
+            ! half of it at nit000+1 (leap frog time step).
             IF ( kt == nit000   ) CALL isf_hdiv_cpl(risfcpl_vol       , phdiv)
-            IF ( kt == nit000+1 ) CALL isf_hdiv_cpl(risfcpl_vol*0.5_wp, phdiv)
+            IF ( kt == nit000+1 ) CALL isf_hdiv_cpl(risfcpl_vol*0.5_wp, phdiv) ! risfcpl_vol*0.5_wp = 0.5 * (risfcpl_vol_b + risfcpl_vol_n)
             !
             ! correct divergence every time step to remove any trend due to coupling