Changeset 14412

Timestamp:

2021-02-05T13:46:51+01:00 (4 years ago)

Author:

agn

Message:

revert to simple hmle=zmld

File:

• NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser_4.0/src/OCE/ZDF/zdfosm.F90 (modified) (4 diffs)

NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser_4.0/src/OCE/ZDF/zdfosm.F90

@@ -640 +640 @@
          CALL zdf_osm_external_gradients( mld_prof, zdtdz_mle_ext, zdsdz_mle_ext, zdbdz_mle_ext )
          CALL zdf_osm_osbl_state_fk( lpyc, lflux, lmle, zwb_fk )
-         CALL zdf_osm_mle_parameters( mld_prof, hmle, zhmle, zvel_mle, zdiff_mle )
+         CALL zdf_osm_mle_parameters( zmld, mld_prof, hmle, zhmle, zvel_mle, zdiff_mle )
       ELSE    ! ln_osm_mle
 ! FK not selected, Boundary Layer only.
@@ -2578 +2578 @@
       
  END SUBROUTINE zdf_osm_zmld_horizontal_gradients
-  SUBROUTINE zdf_osm_mle_parameters( mld_prof, hmle, zhmle, zvel_mle, zdiff_mle )
+  SUBROUTINE zdf_osm_mle_parameters( zmld,  mld_prof, hmle, zhmle, zvel_mle, zdiff_mle )
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE zdf_osm_mle_parameters  ***
@@ -2589 +2589 @@
       !!             Fox-Kemper and Ferrari, JPO, 38, 1166-1179, 2008
 
+      REAL(wp), DIMENSION(jpi,jpj)     :: zmld ! ==  estimated FK BLD used for MLE horiz gradients  == !
       INTEGER, DIMENSION(jpi,jpj)      :: mld_prof
       REAL(wp), DIMENSION(jpi,jpj)     :: hmle, zhmle, zwb_fk, zvel_mle, zdiff_mle
@@ -2631 +2632 @@
            hmle(ji,jj) = MAX(MIN(hmle(ji,jj), ht_n(ji,jj)),  gdepw_n(ji,jj,4))
            IF(ln_osm_hmle_limit) hmle(ji,jj) = MIN(hmle(ji,jj), rn_osm_hmle_limit*hbl(ji,jj) )
+           hmle(ji,jj) = zmld(ji,jj)
         END DO
       END DO