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

Timestamp:

2019-08-05T19:14:52+02:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF_nemo : add comments, renaming file (AGRIF), add isfload module (ticket #2142)

File:

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

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

@@ -1 +1 @@
 MODULE traisf
    !!==============================================================================
-   !!                       ***  MODULE  trasbc  ***
-   !! Ocean active tracers:  surface boundary condition
+   !!                       ***  MODULE  traisf  ***
+   !! Ocean active tracers:  ice shelf boundary condition
    !!==============================================================================
    !! History :    4.0  !  2019-09  (P. Mathiot) original file
@@ -38 +38 @@
    SUBROUTINE tra_isf ( kt )
       !!----------------------------------------------------------------------
-      !!                  ***  ROUTINE tra_sbc  ***
+      !!                  ***  ROUTINE tra_isf  ***
       !!                   
-      !! ** Purpose :   Compute the tracer surface boundary condition trend of
-      !!      (flux through the interface, concentration/dilution effect)
-      !!      and add it to the general trend of tracer equations.
+      !! ** Purpose :  Compute the temperature trend due to the ice shelf melting (qhoce + qhc)
       !!
-      !! ** Method :   The (air+ice)-sea flux has two components: 
-      !!      (1) Fext, external forcing (i.e. flux through the (air+ice)-sea interface); 
-      !!      (2) Fwe , tracer carried with the water that is exchanged with air+ice. 
-      !!               The input forcing fields (emp, rnf, sfx, isf) contain Fext+Fwe,
-      !!             they are simply added to the tracer trend (tsa).
-      !!               In linear free surface case (ln_linssh=T), the volume of the
-      !!             ocean does not change with the water exchanges at the (air+ice)-sea
-      !!             interface. Therefore another term has to be added, to mimic the
-      !!             concentration/dilution effect associated with water exchanges.
-      !!
-      !! ** Action  : - Update tsa with the surface boundary condition trend 
-      !!              - send trends to trdtra module for further diagnostics(l_trdtra=T)
+      !! ** Action  : - update tsa for cav, par and cpl case
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt   ! ocean time-step index
@@ -62 +49 @@
       IF( ln_timing )   CALL timing_start('tra_sbc')
       !
-      !----------------------------------------
-      !       Ice Shelf effects (ISF)
-      !     tbl treated as in Losh (2008) JGR
-      !----------------------------------------
+      ! cavity case
       IF ( ln_isfcav_mlt ) CALL tra_isf_mlt(misfkt_cav, misfkb_cav, rhisf_tbl_cav, rfrac_tbl_cav, risf_cav_tsc, risf_cav_tsc_b, tsa)
+      !
+      ! parametrisation case
       IF ( ln_isfpar_mlt ) CALL tra_isf_mlt(misfkt_par, misfkb_par, rhisf_tbl_par, rfrac_tbl_par, risf_par_tsc, risf_par_tsc_b, tsa)
       !
+      ! ice sheet coupling case
       !CALL tra_isf_cpl()
       !
@@ -77 +64 @@
    SUBROUTINE tra_isf_mlt(ktop, kbot, phtbl, pfrac, ptsc, ptsc_b, pts)
       !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE tra_isf_mlt  ***
+      !!
+      !! *** Purpose :  Compute the temperature trend due to the ice shelf melting (qhoce + qhc) for cav or par case
+      !!
+      !! *** Action :: Update tsa with the surface boundary condition trend 
+      !!
+      !!----------------------------------------------------------------------
       REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(inout) :: pts
       !!----------------------------------------------------------------------
@@ -82 +76 @@
       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   
+      INTEGER                      :: ikt, ikb  ! top and bottom level of the tbl
+      REAL(wp), DIMENSION(jpi,jpj) :: ztc       ! total ice shelf tracer trend
       !!----------------------------------------------------------------------
       !
-      INTEGER  ::   ji, jj, jk, jn              ! dummy loop indices  
-      INTEGER  ::   ikt, ikb                    ! local integers
-      REAL(wp) :: zeps = 1.e-20
-      REAL(wp), DIMENSION(jpi,jpj) :: ztc
-      !!----------------------------------------------------------------------
-      !
-      !----------------------------------------
-      !       Ice Shelf effects (ISF)
-      !     tbl treated as in Losh (2008) JGR
-      !----------------------------------------
-      ! 
       ! compute 2d total trend due to isf
       ztc(:,:) = 0.5_wp * ( ptsc(:,:,jp_tem) + ptsc_b(:,:,jp_tem) ) / phtbl(:,:)