Changeset 11541 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcav.F90

Timestamp:

2019-09-12T18:41:17+02:00 (5 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF: simplify use of ln_isf, add extra comments + minor changes (ticket #2142)

File:

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

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

@@ -1 +1 @@
 MODULE isfcav
    !!======================================================================
-   !!                       ***  MODULE  sbcisf  ***
-   !! Surface module :  update surface ocean boundary condition under ice
-   !!                   shelf
+   !!                       ***  MODULE  isfcav  ***
+   !! Ice shelf cavity module :  update ice shelf melting under ice
+   !!                            shelf
    !!======================================================================
    !! History :  3.2  !  2011-02  (C.Harris  ) Original code isf cav
-   !!            X.X  !  2006-02  (C. Wang   ) Original code bg03
    !!            3.4  !  2013-03  (P. Mathiot) Merging + parametrization
+   !!            4.1  !  2019-09  (P. Mathiot) Split ice shelf cavity and ice shelf parametrisation
    !!----------------------------------------------------------------------
 
    !!----------------------------------------------------------------------
-   !!   sbc_isf       : update sbc under ice shelf
+   !!   isf_cav       : update ice shelf melting under ice shelf
    !!----------------------------------------------------------------------
    USE oce            ! ocean dynamics and tracers
-   USE isf            !
-   USE isftbl         !
-   USE isfcavmlt
-   USE isfcavgam
-   USE isfdiags
+   USE isf            ! ice shelf public variables
+   USE isftbl         ! ice shelf top boundary layer properties
+   USE isfcavmlt      ! ice shelf melt formulation
+   USE isfcavgam      ! ice shelf melt exchange coeficient
+   USE isfdiags       ! ice shelf diags
    USE dom_oce        ! ocean space and time domain
    USE phycst         ! physical constants
@@ -26 +26 @@
    USE iom            ! I/O library
    USE fldread        ! read input field at current time step
-   USE lbclnk         !
+   USE lbclnk         ! lbclnk
 
    IMPLICIT NONE
@@ -43 +43 @@
    SUBROUTINE isf_cav( kt, ptsc, pqfwf )
       !!---------------------------------------------------------------------
-      !!                     ***  ROUTINE sbc_isf_cav  ***
+      !!                     ***  ROUTINE isf_cav  ***
       !!
       !! ** Purpose :   handle surface boundary condition under ice shelf
       !!
-      !! ** Method  : -
+      !! ** Method  :   based on Mathiot et al. (2017)
       !!
-      !! ** Action  :   utau, vtau : remain unchanged
-      !!                taum, wndm : remain unchanged
-      !!                qns        : update heat flux below ice shelf
-      !!                emp, emps  : update freshwater flux below ice shelf
+      !! ** Action  :   - compute geometry of the Losch top bournary layer (see Losch et al. 2008)
+      !!                - depending on the chooses option
+      !!                   - compute temperature/salt in the tbl
+      !!                   - compute exchange coeficient
+      !!                   - compute heat and fwf fluxes
+      !!                   - output
       !!---------------------------------------------------------------------
       !!-------------------------- OUT --------------------------------------
@@ -68 +70 @@
       REAL(wp), DIMENSION(jpi,jpj) :: zttbl, zstbl
       !!---------------------------------------------------------------------
-      !
-      ! compute misfkb_par, rhisf_tbl
-      rhisf_tbl_cav(:,:) = rn_htbl * mskisf_cav(:,:)
-      CALL isf_tbl_lvl( ht_n, e3t_n, misfkt_cav, misfkb_cav, rhisf_tbl_cav, rfrac_tbl_cav )
       !
       ! compute T/S/U/V for the top boundary layer
@@ -118 +116 @@
             END IF
          END SELECT
-
+         !
       END DO
       !
@@ -145 +143 @@
    SUBROUTINE isf_cav_init
       !!---------------------------------------------------------------------
-      !!                  ***  ROUTINE isf_diags_2dto3d ***
+      !!                  ***  ROUTINE isf_cav_init ***
       !!
-      !! ** Purpose : 
+      !! ** Purpose : initialisation of variable needed to compute melt under an ice shelf
       !!
       !!----------------------------------------------------------------------