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

Timestamp:

2009-07-29T16:03:14+02:00 (15 years ago)

Author:

ctlod

Message:

only cosmetic changes

File:

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

trunk/NEMO/OPA_SRC/ZDF/zdfmxl.F90

@@ -4 +4 @@
    !! Ocean physics: mixed layer depth 
    !!======================================================================
-   !! History :
-   !!   9.0  !  03-08  (G. Madec)  OpenMP/autotasking optimization
+   !! History :  1.0  ! 2003-08  (G. Madec)  original code
+   !!            3.2  ! 2009-07  (S. Masson, G. Madec)  IOM + merge of DO-loop
    !!----------------------------------------------------------------------
    !!   zdf_mxl      : Compute the turbocline and mixed layer depths.
    !!----------------------------------------------------------------------
-   !! * Modules used
    USE oce             ! ocean dynamics and tracers variables
    USE dom_oce         ! ocean space and time domain variables
@@ -20 +19 @@
    PRIVATE
 
-   !! * Routine accessibility
-   PUBLIC zdf_mxl           ! called by step.F90
+   PUBLIC   zdf_mxl    ! called by step.F90
 
-   !! * Shared module variables
-   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)
-      hmlp ,             &  !: mixed layer depth  (rho=rho0+zdcrit) (m)
-      hmlpt                 !: mixed layer depth at t-points (m)
+   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]
 
-   !! * module variables
-   REAL(wp) ::   &
-      avt_c = 5.e-4_wp,  &  ! Kz criterion for the turbocline depth
-      rho_c = 0.01_wp       ! density criterion for mixed layer depth
+   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
 #  include "domzgr_substitute.h90"
    !!----------------------------------------------------------------------
-   !!   OPA 9.0 , LOCEAN-IPSL (2005) 
+   !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 
    !! $Id$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
@@ -51 +44 @@
       !!                   
       !! ** Purpose :   Compute the turbocline depth and the mixed layer depth
-      !!      with density criteria.
+      !!              with density criteria.
       !!
       !! ** Method  :   The turbocline depth is the depth at which the vertical
@@ -58 +51 @@
       !!      value defined locally (avt_c here taken equal to 5 cm/s2)
       !!
-      !! ** Action  :
+      !! ** Action  :   nmln, hmld, hmlp, hmlpt
+      !!----------------------------------------------------------------------
+      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
       !!
-      !! History :
-      !!        !  94-11  (M. Imbard)  Original code
-      !!   8.0  !  96-01  (E. Guilyardi)  sub mixed layer temp.
-      !!   8.1  !  97-07  (G. Madec)  optimization
-      !!   8.5  !  02-06  (G. Madec)  F90: Free form and module
-      !!----------------------------------------------------------------------
-      !! * Arguments
-      INTEGER, INTENT( in ) ::   kt         ! ocean time-step index
-
-      !! * Local declarations
       INTEGER ::   ji, jj, jk     ! dummy loop indices
-      INTEGER ::   ik             ! temporary integer
-      INTEGER, DIMENSION(jpi,jpj) ::   &
-         imld                     ! temporary workspace
+      INTEGER ::   iki, ikn       ! temporary integer
+      INTEGER, DIMENSION(jpi,jpj) ::   imld   ! temporary workspace
       !!----------------------------------------------------------------------
 
@@ -82 +66 @@
       ENDIF
 
-
-      ! 1. Turbocline depth
-      ! -------------------
-      ! last w-level at which avt<avt_c (starting from the bottom jk=jpk)
-      ! (since avt(.,.,jpk)=0, we have jpk=< imld =< 2 )
-      DO jk = jpk, 2, -1
+      ! 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
          DO jj = 1, jpj
             DO ji = 1, jpi
-               IF( avt(ji,jj,jk) < avt_c ) imld(ji,jj) = jk 
+               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
             END DO
          END DO
       END DO
-
-      ! Turbocline depth and sub-turbocline temperature
+      ! depth of the mixing and mixed layers
       DO jj = 1, jpj
          DO ji = 1, jpi
-            ik = imld(ji,jj)
-            hmld (ji,jj) = fsdepw(ji,jj,ik) * tmask(ji,jj,1)
+            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
          END DO
       END DO
       CALL iom_put( "mldturb", hmld )   ! turbocline depth
-
-!!gm idea
-!!   
-!!gm  DO jk = jpk, 2, -1
-!!gm     DO jj = 1, jpj
-!!gm        DO ji = 1, jpi
-!!gm           IF( avt(ji,jj,jk) < avt_c ) hmld(ji,jj) = fsdepw(ji,jj,jk) * tmask(ji,jj,1)
-!!gm        END DO
-!!gm     END DO
-!!gm  END DO
-!!gm
-
-      ! 2. Mixed layer depth
-      ! --------------------
-      ! Initialization to the number of w ocean point mbathy
-      nmln(:,:) = mbathy(:,:)
-
-      ! Last w-level at which rhop>=rho surf+rho_c (starting from jpk-1)
-      ! (rhop defined at t-point, thus jk-1 for w-level just above)
-      DO jk = jpkm1, 2, -1
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-               IF( rhop(ji,jj,jk) > rhop(ji,jj,1) + rho_c )   nmln(ji,jj) = jk
-            END DO
-         END DO
-      END DO
-
-      ! Mixed layer depth
-      DO jj = 1, jpj
-         DO ji = 1, jpi
-            ik = nmln(ji,jj)
-            hmlp (ji,jj) = fsdepw(ji,jj,ik) * tmask(ji,jj,1)
-            hmlpt(ji,jj) = fsdept(ji,jj,ik-1)
-         END DO
-      END DO
-      CALL iom_put( "mld010", hmlp )   ! mixed layer depth
+      CALL iom_put( "mld010" , hmlp )   ! mixed layer depth
       
       IF(ln_ctl)   CALL prt_ctl( tab2d_1=REAL(nmln,wp), clinfo1=' nmln : ', tab2d_2=hmld, clinfo2=' hmld : ', ovlap=1 )
-
+      !
    END SUBROUTINE zdf_mxl