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

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

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

@@ -16 +16 @@
    USE isf            ! Ice shelf variable
    USE isfutils       !
-   USE iscplini       ! Ice sheet coupling
    !
    USE in_out_manager ! I/O manager
@@ -56 +55 @@
       !
       ! ice sheet coupling case
-      !CALL tra_isf_cpl()
+      IF ( ll_isfcpl .AND. kt == nit000 ) CALL tra_isf_cpl(risfcpl_tsc, tsa)
+      !
+      ! ice sheet coupling case (
+      IF ( ll_isfcpl_cons ) CALL tra_isf_cpl(risfcpl_cons_tsc, tsa)
       !
       IF( ln_timing )   CALL timing_stop('tra_isf')
@@ -76 +78 @@
       REAL(wp), DIMENSION(jpi,jpj)     , INTENT(in   ) :: phtbl, pfrac
       REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in   ) :: ptsc , ptsc_b
-      REAL(wp), DIMENSION(jpi,jpj) :: ztc
       !!----------------------------------------------------------------------
       INTEGER                      :: ji,jj,jk  ! loop index   
@@ -107 +108 @@
    END SUBROUTINE tra_isf_mlt
    !
-   !SUBROUTINE tra_isf_cpl
-       !
-!      !----------------------------------------
-!      !        Ice Sheet coupling imbalance correction to have conservation
-!      !----------------------------------------
-!      !
-!      IF( ln_iscpl .AND. ln_iscpl_hsb) THEN         ! input of heat and salt due to river runoff 
-!         DO jk = 1,jpk
-!            DO jj = 2, jpj 
-!               DO ji = fs_2, fs_jpim1
-!                  zdep = 1._wp / e3t_n(ji,jj,jk) 
-!                  tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - htsc_iscpl(ji,jj,jk,jp_tem) * zdep
-!                  tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - htsc_iscpl(ji,jj,jk,jp_sal) * zdep  
-!               END DO  
-!            END DO  
-!         END DO
-!      ENDIF
-!      !
-!   END SUBROUTINE tra_isf_cpl
+   SUBROUTINE tra_isf_cpl( ptsc, ptsa )
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE tra_isf_cpl  ***
+      !!
+      !! *** Action :: Update tsa with the ice shelf coupling trend 
+      !!
+      !!----------------------------------------------------------------------
+      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(inout) :: ptsa
+      !!----------------------------------------------------------------------
+      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) :: ptsc  !( >0 out )
+      !!----------------------------------------------------------------------
+      INTEGER :: jk
+      !!----------------------------------------------------------------------
+      !
+      DO jk = 1,jpk
+         ptsa(:,:,jk,jp_tem) = ptsa(:,:,jk,jp_tem) - ptsc(:,:,jk,jp_tem) * r1_e1e2t(:,:) / e3t_n(:,:,jk)
+         ptsa(:,:,jk,jp_sal) = ptsa(:,:,jk,jp_sal) - ptsc(:,:,jk,jp_sal) * r1_e1e2t(:,:) / e3t_n(:,:,jk)
+      END DO
+      !
+   END SUBROUTINE tra_isf_cpl
    !
 END MODULE traisf