Changeset 12077 for NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/TRA/traisf.F90

Timestamp:

2019-12-05T18:41:39+01:00 (4 years ago)

Author:

mathiot

Message:

include ENHANCE-02_ISF_nemo in UKMO merge branch

File:

• NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/TRA/traisf.F90 (modified) (2 diffs)

NEMO/branches/2019/UKMO_MERGE_2019/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 phycst         ! physical constant
-   USE eosbn2         ! Equation Of State
-   USE isf            ! Ice shelf variable
-   USE isfutils       !
-   !
-   USE in_out_manager ! I/O manager
-   USE iom            ! xIOS server
-   USE timing         ! Timing
+   USE isf_oce                                     ! Ice shelf variables
+   USE dom_oce , ONLY : e3t, 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
@@ -65 +60 @@
       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(Kmm, risfcpl_tsc       , pts(:,:,:,:,Krhs))
          IF ( kt == nit000+1) CALL tra_isf_cpl(Kmm, risfcpl_tsc*0.5_wp, pts(:,:,:,:,Krhs))