Changeset 3959 for branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

Timestamp:

2013-07-08T18:46:13+02:00 (11 years ago)

Author:

agn

Message:

Gurvan's FK + separate rn_rho_c_mle

File:

• branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90 (modified) (4 diffs)

branches/2013/dev_r3948_NOC_FK/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

@@ -24 +24 @@
 
    PUBLIC   zdf_mxl       ! called by step.F90
+
+   REAL(wp), PUBLIC ::   rho_c = 0.01_wp    ! density criterion for mixed layer depth
+   REAL(wp), PUBLIC ::   avt_c = 5.e-4_wp   ! Kz criterion for the turbocline depth
 
    INTEGER , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   nmln    !: number of level in the mixed layer (used by TOP)
@@ -63 +66 @@
       !! ** Method  :   The mixed layer depth is the shallowest W depth with 
       !!      the density of the corresponding T point (just bellow) bellow a
-      !!      given value defined locally as rho(10m) + zrho_c
+      !!      given value defined locally as rho(10m) + rho_c
       !!               The turbocline depth is the depth at which the vertical
       !!      eddy diffusivity coefficient (resulting from the vertical physics
@@ -76 +79 @@
       INTEGER  ::   iikn, iiki          ! temporary integer within a do loop
       INTEGER, POINTER, DIMENSION(:,:) ::   imld                ! temporary workspace
-      REAL(wp) ::   zrho_c = 0.01_wp    ! density criterion for mixed layer depth
-      REAL(wp) ::   zavt_c = 5.e-4_wp   ! Kz criterion for the turbocline depth
       !!----------------------------------------------------------------------
       !
@@ -98 +99 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               IF( rhop(ji,jj,jk) > rhop(ji,jj,nla10) + zrho_c )   nmln(ji,jj) = jk      ! Mixed layer
-               IF( avt (ji,jj,jk) < zavt_c                     )   imld(ji,jj) = jk      ! Turbocline 
+               IF( rhop(ji,jj,jk) > rhop(ji,jj,nla10) + rho_c )   nmln(ji,jj) = jk      ! Mixed layer
+               IF( avt (ji,jj,jk) < avt_c                     )   imld(ji,jj) = jk      ! Turbocline 
             END DO
          END DO