Changeset 12062 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/TRA

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/TRA/traisf.F90 (modified) (2 diffs)

NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/TRA/traisf.F90

@@ -10 +10 @@
    !!   tra_isf       : update the tracer trend at ocean surface
    !!----------------------------------------------------------------------
-   USE oce            ! ocean dynamics and active tracers
-   USE dom_oce        ! ocean space domain variables
-   USE isf_oce        ! Ice shelf variable
-   USE phycst         ! physical constant
-   USE eosbn2         ! Equation Of State
-   USE isfutils       !
-   !
-   USE in_out_manager ! I/O manager
-   USE iom            ! xIOS server
-   USE timing         ! Timing
+   USE isf_oce                                     ! Ice shelf variables
+   USE oce     , ONLY : tsa                        ! ocean dynamics and active tracers
+   USE dom_oce , ONLY : e3t_n, r1_e1e2t            ! ocean space domain variables
+   USE isfutils, ONLY : debug                      ! debug option
+   USE timing  , ONLY : timing_start, timing_stop  ! Timing
+   USE in_out_manager                              ! I/O manager
 
    IMPLICIT NONE
@@ -63 +59 @@
       IF ( ln_isfcpl ) THEN
          !
+         ! 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).
+         ! in accordance to this, the heat content flux due to injected water need to be added in the temperature and salt trend
+         ! at time step nit000 and nit000+1
          IF ( kt == nit000  ) CALL tra_isf_cpl(risfcpl_tsc       , tsa)
          IF ( kt == nit000+1) CALL tra_isf_cpl(risfcpl_tsc*0.5_wp, tsa)