Changeset 255 for trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90

Timestamp:

2005-09-02T17:06:49+02:00 (19 years ago)

Author:

opalod

Message:

nemo_v1_update_002 : CT : Integration of the KPP turbulent closure scheme

File:

• trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90 (modified) (9 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90

@@ -15 +15 @@
    USE dom_oce         ! ocean space and time domain variables
    USE zdf_oce         ! ocean vertical physics variables
+   USE zdfkpp          ! KPP vertical mixing
    USE in_out_manager  ! I/O manager
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
@@ -24 +25 @@
    PUBLIC zdf_evd      ! called by step.F90
 
+   !! * Substitutions
+#  include "domzgr_substitute.h90"
    !!----------------------------------------------------------------------
    !!   OPA 9.0 , LOCEAN-IPSL (2005) 
@@ -44 +47 @@
       !!
       !! ** Action  :   Update avt, avmu, avmv in statically instable cases
-      !!
+      !!                and avt_evd which is avt due to convection
       !! References :
       !!      Lazar, A., these de l'universite Paris VI, France, 1997
@@ -50 +53 @@
       !!   7.0  !  97-06  (G. Madec, A. Lazar)  Original code
       !!   8.5  !  02-06  (G. Madec)  F90: Free form and module
+      !!   9.0  !  05-06  (C. Ethe) KPP parameterization
       !!----------------------------------------------------------------------
       !! * Arguments
@@ -65 +69 @@
       ENDIF
 
+      ! Initialisation of avt_evd (vertical diffusion due to convection) to avt and avmu_evd to avmu
+      avt_evd  (:,:,:) = avt(:,:,:) 
+      avmu_evd (:,:,:) = avmu(:,:,:) 
+
       SELECT CASE ( nevdm )
  
@@ -75 +83 @@
             jj = 1                     ! big loop forced
             DO ji = jpi+2, jpij   
+#   if defined key_zdfkpp
+!! no implicit mixing in the boundary layer with KPP
+               IF( ( rn2(ji,jj,jk) <= -1.e-12 ) .AND. ( fsdepw(ji,jj,jk) > hkpp(ji,jj) ) ) THEN
+#   else
                IF( rn2(ji,jj,jk) <= -1.e-12 ) THEN
+#   endif
                   avt (ji  ,jj  ,jk) = avevd * tmask(ji  ,jj  ,jk)
                   avmu(ji  ,jj  ,jk) = avevd * umask(ji  ,jj  ,jk)
@@ -86 +99 @@
             DO jj = 2, jpj             ! no vector opt.
                DO ji = 2, jpi
-                  IF( rn2(ji,jj,jk) <= -1.e-12 ) THEN
+#   if defined key_zdfkpp
+!! no implicit mixing in the boundary layer with KPP
+               IF( ( rn2(ji,jj,jk) <= -1.e-12 ) .AND. ( fsdepw(ji,jj,jk) > hkpp(ji,jj) ) ) THEN
+#   else
+               IF( rn2(ji,jj,jk) <= -1.e-12 ) THEN
+#   endif
                      avt (ji  ,jj  ,jk) = avevd * tmask(ji  ,jj  ,jk)
                      avmu(ji  ,jj  ,jk) = avevd * umask(ji  ,jj  ,jk)
@@ -114 +132 @@
             jj = 1                     ! big loop forced
             DO ji = 1, jpij   
+#   if defined key_zdfkpp
+!! no implicit mixing in the boundary layer with KPP
+               IF( ( rn2(ji,jj,jk) <= -1.e-12 ) .AND. ( fsdepw(ji,jj,jk) > hkpp(ji,jj) ) ) &              
+                  avt(ji,jj,jk) = avevd * tmask(ji,jj,jk)
+#   else
                IF( rn2(ji,jj,jk) <= -1.e-12 )   avt(ji,jj,jk) = avevd * tmask(ji,jj,jk)
+#   endif
             END DO
 #   else
             DO jj = 1, jpj             ! loop over the whole domain (no lbc_lnk call)
                DO ji = 1, jpi
+#   if defined key_zdfkpp
+!! no implicit mixing in the boundary layer with KPP
+               IF( ( rn2(ji,jj,jk) <= -1.e-12 ) .AND. ( fsdepw(ji,jj,jk) > hkpp(ji,jj) ) ) &          
+                  avt(ji,jj,jk) = avevd * tmask(ji,jj,jk)
+#   else
                   IF( rn2(ji,jj,jk) <= -1.e-12 )   avt(ji,jj,jk) = avevd * tmask(ji,jj,jk)
+#   endif
                END DO
             END DO
@@ -128 +158 @@
       END SELECT 
 
+      ! update of avt_evd and avmu_evd
+      avt_evd  (:,:,:) = avt (:,:,:)  - avt_evd  (:,:,:) 
+      avmu_evd (:,:,:) = avmu(:,:,:)  - avmu_evd (:,:,:) 
+
    END SUBROUTINE zdf_evd