Changeset 11541 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isftbl.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/isftbl.F90 (modified) (8 diffs)

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

@@ -21 +21 @@
    PRIVATE
 
-   PUBLIC isf_tbl, isf_tbl_avg, isf_tbl_lvl, isftbl_ktop, isftbl_kbot
+   PUBLIC isf_tbl, isf_tbl_avg, isf_tbl_lvl, isf_tbl_ktop, isf_tbl_kbot
 
 CONTAINS
@@ -30 +30 @@
       !!
       !! ** Purpose : compute mean T/S/U/V in the boundary layer at T- point
+      !!
+      !! ** Method : Average properties over a specific thickness
+      !!
+      !! ** Reference : inspired from : Losch, Modeling ice shelf cavities in a z coordinate ocean general circulation model
+      !!                https://doi.org/10.1029/2007JC004368 , 2008
       !!
       !!--------------------------------------------------------------------
@@ -103 +108 @@
    SUBROUTINE isf_tbl_avg( ktop, kbot, phtbl, pfrac, pe3, pvarin, pvarout )
       !!--------------------------------------------------------------------
-      !!                  ***  ROUTINE isf_tbl_lvl  ***
+      !!                  ***  ROUTINE isf_tbl_avg  ***
       !!
       !! ** Purpose : compute mean property in the boundary layer
@@ -145 +150 @@
       !!                  ***  ROUTINE isf_tbl_lvl  ***
       !!
-      !! ** Purpose : - compute bottom level fully included in the top boundary layer
+      !! ** Purpose : - compute bottom level off the top boundary layer
       !!              - thickness of the top boundary layer
+      !!              - fraction of the bottom level affected by the tbl
       !!
       !!---------------------------------------------------------------------
@@ -180 +186 @@
       !
       ! get ktbl
-      CALL isftbl_kbot(ktop, phtbl, pe3, kbot)
+      CALL isf_tbl_kbot(ktop, phtbl, pe3, kbot)
       !
       ! get pfrac
@@ -197 +203 @@
    END SUBROUTINE isf_tbl_lvl
    !
-   SUBROUTINE isftbl_kbot(ktop, phtbl, pe3, kbot)
-      !!--------------------------------------------------------------------
-      !!                  ***  ROUTINE isf_tbl_lvl  ***
+   SUBROUTINE isf_tbl_kbot(ktop, phtbl, pe3, kbot)
+      !!--------------------------------------------------------------------
+      !!                  ***  ROUTINE isf_tbl_bot  ***
       !!
       !! ** Purpose : compute bottom level of the isf top boundary layer
@@ -232 +238 @@
       END DO
       !
-   END SUBROUTINE isftbl_kbot
-      !
-   SUBROUTINE isftbl_ktop(pdep, ktop)
-      !!--------------------------------------------------------------------
-      !!                  ***  ROUTINE isf_tbl_lvl  ***
+   END SUBROUTINE isf_tbl_kbot
+      !
+   SUBROUTINE isf_tbl_ktop(pdep, ktop)
+      !!--------------------------------------------------------------------
+      !!                  ***  ROUTINE isf_tbl_top  ***
       !!
       !! ** Purpose : compute top level of the isf top boundary layer in case of an ice shelf parametrisation
@@ -263 +269 @@
       END DO
       !
-   END SUBROUTINE isftbl_ktop
+   END SUBROUTINE isf_tbl_ktop
 
 END MODULE isftbl