Changeset 11931 for NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcavmlt.F90

Timestamp:

2019-11-19T18:30:57+01:00 (4 years ago)

Author:

mathiot

Message:

ENHANCE-02_ISF_nemo: add comments, improve memory usage of ln_isfcpl_cons option, fix issue in ISOMIP+ configuration

File:

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

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

@@ -3 +3 @@
    !!                       ***  MODULE  isfcavmlt  ***
    !! ice shelf module :  update surface ocean boundary condition under ice
-   !!                   shelf
+   !!                   shelves
    !!======================================================================
    !! History :  4.0  !  2019-09  (P. Mathiot) Original code
@@ -9 +9 @@
 
    !!----------------------------------------------------------------------
-   !!   isfcav_mlt    : update surface ocean boundary condition under ice shelf
+   !!   isfcav_mlt    : compute or read ice shelf fwf/heat fluxes in the ice shelf cavity
    !!----------------------------------------------------------------------
 
@@ -23 +23 @@
    USE in_out_manager              ! I/O manager
    USE iom        , ONLY: iom_put  ! I/O library
-   USE fldread    , ONLY: fld_read !
+   USE fldread    , ONLY: fld_read, FLD, FLD_N !
    USE lib_fortran, ONLY: glob_sum !
    USE lib_mpp    , ONLY: ctl_stop !
@@ -51 +51 @@
       !! ** Purpose    : compute or read ice shelf fwf/heat fluxes in the ice shelf cavity
       !!
-      !!---------------------------------------------------------------------
       !!-------------------------- OUT -------------------------------------
       REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) :: pqhc, pqoce, pqfwf  ! heat and fwf fluxes
@@ -58 +57 @@
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pgt  , pgs    ! gamma t and gamma s
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pttbl, pstbl  ! top boundary layer tracer
-      !!---------------------------------------------------------------------
       !!---------------------------------------------------------------------
       !
@@ -110 +108 @@
       !!--------------------------------------------------------------------
       !
-      ! Calculate freezing temperature
+      ! Compute freezing temperature
       CALL eos_fzp( pstbl(:,:), ztfrz(:,:), risfdep(:,:) )
       !
@@ -131 +129 @@
       !!----------------------------------------------------------------------
       !!
-      !!                          ***  ROUTINE isfcav_mlt_spe  ***
+      !!                          ***  ROUTINE isfcav_mlt_2eq  ***
       !!
       !! ** Purpose    : Compute ice shelf fwf/heqt fluxes using ISOMIP formulation (Hunter et al., 2006)
@@ -144 +142 @@
       !!                 Tech.  Rep.  June,  Antarctic  Climate  &  Ecosystems  Cooperative  Research  Centre,  available  at:  
       !!                 http://staff.acecrc.org.au/~bkgalton/ISOMIP/test_cavities.pdf (last access: 21 July 2016), 2006.
-      !!---------------------------------------------------------------------
+      !!
       !!-------------------------- OUT -------------------------------------
       REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) :: pqhc, pqoce, pqfwf  ! hean content, ocean-ice heat and fwf fluxes
@@ -192 +190 @@
       !!                   MISMIP v. 3 (MISMIP +), ISOMIP v. 2 (ISOMIP +) and MISOMIP v. 1 (MISOMIP1), 
       !!                   Geosci. Model Dev., 9, 2471-2497, https://doi.org/10.5194/gmd-9-2471-2016, 2016. 
-      !!---------------------------------------------------------------------
+      !!
       !!-------------------------- OUT -------------------------------------
       REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) :: pqhc, pqoce, pqfwf  ! latent heat and fwf fluxes
@@ -308 +306 @@
       !
       CALL iom_put('isftfrz_cav', ztfrz * mskisf_cav(:,:) )
+      !
    END SUBROUTINE isfcav_mlt_oasis