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

Timestamp:

2019-08-08T16:02:49+02:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF_nemo : add UKESM ice sheet coupling method (ticket #2142)

File:

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

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

@@ -4 +4 @@
    USE isf
    USE phycst
+   USE in_out_manager
 
    IMPLICIT NONE
@@ -13 +14 @@
 CONTAINS
 
-   SUBROUTINE isf_hdiv( phdiv )
+   SUBROUTINE isf_hdiv( kt, phdiv )
       !!----------------------------------------------------------------------
       !!                  ***  SUBROUTINE isf_hdiv  ***
@@ -25 +26 @@
       REAL(wp), DIMENSION(:,:,:), INTENT( inout ) ::   phdiv   ! horizontal divergence
       !!----------------------------------------------------------------------
+      INTEGER, INTENT(in) :: kt
       !
       ! ice shelf cavity contribution
@@ -33 +35 @@
       !
       ! ice sheet coupling contribution (if conservation needed)
-      !IF ( ln_iscpl_hsb  ) CALL isf_hdiv_cpl(hdiv_iscpl, phdivn)
+      IF ( ll_isfcpl ) THEN
+         !
+         ! correct divergence only for the first time step
+         IF ( kt == nit000 ) CALL isf_hdiv_cpl(risfcpl_vol, phdiv)
+         !
+         ! correct divergence every time step to remove any trend due to coupling
+         IF ( ll_isfcpl_cons ) CALL isf_hdiv_cpl(risfcpl_cons_vol, phdiv)
+         !
+      END IF
       !
    END SUBROUTINE isf_hdiv
@@ -81 +91 @@
    END SUBROUTINE isf_hdiv_mlt
 
+   SUBROUTINE isf_hdiv_cpl(pqvol, phdiv)
+      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: phdiv
+      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in   ) :: pqvol
+   
+   INTEGER :: jk
+      
+      DO jk=1,jpk 
+         phdiv(:,:,jk) =  phdiv(:,:,jk) + pqvol(:,:,jk) * r1_e1e2t(:,:) / e3t_n(:,:,jk)
+      END DO
+
+   END SUBROUTINE
+
 END MODULE isfhdiv