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

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

@@ -18 +18 @@
    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 !
@@ -49 +49 @@
       !!                        1 : Specified melt flux
       !!                        2 : Beckmann & Goose parameterization
-      !!----------------------------------------------------------------------
+      !!
       !!-------------------------- OUT -------------------------------------
       REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pqfwf, pqoce, pqhc  ! fresh water, ice-ocean heat and heat content fluxes
@@ -81 +81 @@
       !!              data read into a forcing files.
       !!
-      !!----------------------------------------------------------------------
       !!-------------------------- OUT -------------------------------------
       REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pqhc, pqfwf, pqoce  ! fresh water and ice-ocean heat fluxes
@@ -105 +104 @@
       pqoce(:,:) =   pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux
       pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux 
+      !
       CALL iom_put('isftfrz_par', ztfrz )
       !
@@ -155 +155 @@
       ! output thermal driving
       CALL iom_put('isfthermald_par',( ztfrz(:,:) - ztavg(:,:) ) * mskisf_par(:,:))
+      !
+      ! output freezing point used to define the thermal driving and heat content fluxes
       CALL iom_put('isftfrz_par', ztfrz )
       !
@@ -206 +208 @@
       zfwf(:,:) = zfwf(:,:) * zfwf_oasis / zfwf_fld
       ! 
-      ! i3. -----------Define fwf and qoce
+      ! 3. -----------Define fwf and qoce
       ! ocean heat flux is assume to be equal to the latent heat
       pqfwf(:,:) =   zfwf(:,:)                         ! fwf                ( >0 out )