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

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

@@ -14 +14 @@
    !!   isfpar       : compute ice shelf melt using a prametrisation of ice shelf cavities
    !!----------------------------------------------------------------------
-   !USE oce            ! ocean dynamics and tracers
    USE isf            ! ice shelf
    !
@@ -23 +22 @@
    USE isfdiags , ONLY: isf_diags_flx  ! ice shelf diags subroutine
    !
-   USE dom_oce  , ONLY: jpi,jpj, bathy ! ocean space and time domain
+   USE dom_oce  , ONLY: bathy          ! ocean space and time domain
+   USE par_oce  , ONLY: jpi,jpj        ! ocean space and time domain
    USE phycst   , ONLY: r1_rau0_rcp    ! physical constants
    !
@@ -62 +62 @@
       CALL isfpar_mlt( kt, zqhc, zqoce, pqfwf  )
       !
-      ! compute heat and water flux  (change signe directly in the melt subroutine)
+      ! compute heat and water flux ( > 0 out )
       pqfwf(:,:) = pqfwf(:,:) * mskisf_par(:,:)
       zqoce(:,:) = zqoce(:,:) * mskisf_par(:,:)
       zqhc (:,:) = zqhc(:,:)  * mskisf_par(:,:)
       !
-      ! compute heat content flux
+      ! compute heat content flux ( > 0 out )
       zqlat(:,:) = pqfwf(:,:) * rLfusisf    ! 2d latent heat flux (W/m2)
       !
-      ! total heat flux
+      ! total heat flux ( > 0 out )
       zqh(:,:) = ( zqhc (:,:) + zqoce(:,:) )
       !
@@ -107 +107 @@
       CALL isf_alloc_par()
       !
-      ! par
+      ! initialisation
       misfkt_par(:,:)     = 1         ; misfkb_par(:,:)       = 1         
       rhisf_tbl_par(:,:)  = 1e-20     ; rfrac_tbl_par(:,:)    = 0.0_wp