Changeset 1577 for trunk/NEMO/OPA_SRC/ZDF/zdfmxl.F90

Timestamp:

2009-08-04T16:56:59+02:00 (15 years ago)

Author:

smasson

Message:

update diahth and zdfmxl, see ticket:468

File:

• trunk/NEMO/OPA_SRC/ZDF/zdfmxl.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdfmxl.F90

@@ -5 +5 @@
    !!======================================================================
    !! History :  1.0  ! 2003-08  (G. Madec)  original code
-   !!            3.2  ! 2009-07  (S. Masson, G. Madec)  IOM + merge of DO-loop
+   !!            3.2  ! 2009-07  (S. Masson, G. Madec)  IOM + 10m ref.
    !!----------------------------------------------------------------------
-   !!   zdf_mxl      : Compute the turbocline and mixed layer depths.
+   !!   zdf_mxl       : Compute mixed layer depth and level
    !!----------------------------------------------------------------------
    USE oce             ! ocean dynamics and tracers variables
    USE dom_oce         ! ocean space and time domain variables
-   USE zdf_oce         ! ocean vertical physics
    USE in_out_manager  ! I/O manager
    USE prtctl          ! Print control
@@ -22 +21 @@
 
    INTEGER , PUBLIC, DIMENSION(jpi,jpj) ::   nmln    !: number of level in the mixed layer
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   hmld    !: mixing layer depth (turbocline)      [m]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   hmlp    !: mixed layer depth  (rho=rho0+zdcrit) [m]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   hmlpt   !: mixed layer depth at t-points        [m]
-
-   REAL(wp) ::   avt_c = 5.e-4_wp   ! Kz criterion for the turbocline depth
-   REAL(wp) ::   rho_c = 0.01_wp    ! density criterion for mixed layer depth
 
    !! * Substitutions
@@ -43 +38 @@
       !!                  ***  ROUTINE zdfmxl  ***
       !!                   
-      !! ** Purpose :   Compute the turbocline depth and the mixed layer depth
-      !!              with density criteria.
+      !! ** Purpose :   the mixed layer depth with density criteria.
       !!
-      !! ** Method  :   The turbocline depth is the depth at which the vertical
-      !!      eddy diffusivity coefficient (resulting from the vertical physics
-      !!      alone, not the isopycnal part, see trazdf.F) fall below a given
-      !!      value defined locally (avt_c here taken equal to 5 cm/s2)
+      !! ** 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
       !!
-      !! ** Action  :   nmln, hmld, hmlp, hmlpt
+      !! ** Action  :   nmln, hmlp, hmlpt
       !!----------------------------------------------------------------------
       INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
       !!
-      INTEGER ::   ji, jj, jk     ! dummy loop indices
-      INTEGER ::   iki, ikn       ! temporary integer
-      INTEGER, DIMENSION(jpi,jpj) ::   imld   ! temporary workspace
+      INTEGER  ::   ji, jj, jk         ! dummy loop indices
+      INTEGER  ::   iikn               ! temporary integer within a do loop
+      REAL(wp) ::   zrho_c = 0.01_wp   ! density criterion for mixed layer depth
       !!----------------------------------------------------------------------
 
@@ -67 +60 @@
 
       ! w-level of the mixing and mixed layers
-      nmln(:,:) = mbathy(:,:)      ! Initialization to the number of w ocean point mbathy
-      imld(:,:) = mbathy(:,:)
-      DO jk = jpkm1, 2, -1         ! from the bottom to the surface
+      nmln(:,:) = mbathy(:,:)          ! Initialization to the number of w ocean point mbathy
+      DO jk = jpkm1, nlb10, -1         ! from the bottom to nlb10
          DO jj = 1, jpj
             DO ji = 1, jpi
-               IF( avt (ji,jj,jk) < avt_c                 )   imld(ji,jj) = jk      ! Turbocline 
-               IF( rhop(ji,jj,jk) > rhop(ji,jj,1) + rho_c )   nmln(ji,jj) = jk      ! Mixed layer
+               IF( rhop(ji,jj,jk) > rhop(ji,jj,nla10) + zrho_c )   nmln(ji,jj) = jk      ! Mixed layer
             END DO
          END DO
@@ -80 +71 @@
       DO jj = 1, jpj
          DO ji = 1, jpi
-            iki = imld(ji,jj)
-            ikn = nmln(ji,jj)
-            hmld (ji,jj) = fsdepw(ji,jj,iki) * tmask(ji,jj,1)      ! Turbocline depth 
-            hmlp (ji,jj) = fsdepw(ji,jj,ikn) * tmask(ji,jj,1)      ! Mixed layer depth
-            hmlpt(ji,jj) = fsdept(ji,jj,ikn-1)                     ! depth of the last T-point inside the mixed layer
+            iikn = nmln(ji,jj)
+            hmlp (ji,jj) = fsdepw(ji,jj,iikn  ) * tmask(ji,jj,1)    ! Mixed layer depth
+            hmlpt(ji,jj) = fsdept(ji,jj,iikn-1)                     ! depth of the last T-point inside the mixed layer
          END DO
       END DO
-      CALL iom_put( "mldturb", hmld )   ! turbocline depth
-      CALL iom_put( "mld010" , hmlp )   ! mixed layer depth
+      CALL iom_put( "mldr10_1" , hmlp )   ! mixed layer depth
       
-      IF(ln_ctl)   CALL prt_ctl( tab2d_1=REAL(nmln,wp), clinfo1=' nmln : ', tab2d_2=hmld, clinfo2=' hmld : ', ovlap=1 )
+      IF(ln_ctl)   CALL prt_ctl( tab2d_1=REAL(nmln,wp), clinfo1=' nmln : ', tab2d_2=hmlp, clinfo2=' hmlp : ', ovlap=1 )
       !
    END SUBROUTINE zdf_mxl