Changeset 6403 for trunk/NEMOGCM/NEMO/TOP_SRC

Timestamp:

2016-03-25T17:24:35+01:00 (8 years ago)

Author:

cetlod

Message:

trunk:new developments already included in 3.6 stable, see points 1, 2 and 4 of ticket #1678

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf.F90 (modified) (5 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf.F90

@@ -40 +40 @@
    REAL(wp), PUBLIC ::   rn_ahtrc_0          !:   laplacian diffusivity coefficient for passive tracer [m2/s]
    REAL(wp), PUBLIC ::   rn_bhtrc_0          !: bilaplacian      -          --     -       -   [m4/s]
+   REAL(wp), PUBLIC ::   rn_fact_lap         !: Enhanced zonal diffusivity coefficent in the equatorial domain
    !
    !                      !!: ** lateral mixing namelist (nam_trcldf) **
@@ -64 +65 @@
       INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
       !
-      INTEGER            :: jn
+      INTEGER            :: ji, jj, jk, jn
+      REAL(wp)           :: zdep
       CHARACTER (len=22) :: charout
       REAL(wp), POINTER, DIMENSION(:,:,:)   ::   zahu, zahv
@@ -76 +78 @@
          ztrtrd(:,:,:,:)  = tra(:,:,:,:)
       ENDIF
-      !
-      !                                        !* set the lateral diffusivity coef. for passive tracer      
+      !                                  !* set the lateral diffusivity coef. for passive tracer      
       CALL wrk_alloc( jpi,jpj,jpk,   zahu, zahv )
-      zahu(:,:,:) = rldf * ahtu(:,:,:)
+      zahu(:,:,:) = rldf * ahtu(:,:,:) 
       zahv(:,:,:) = rldf * ahtv(:,:,:)
-
+      !                                  !* Enhanced zonal diffusivity coefficent in the equatorial domain
+      DO jk= 1, jpk
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               IF( gdept_n(ji,jj,jk) > 200. .AND. gphit(ji,jj) < 5. .AND. gphit(ji,jj) > -5. ) THEN
+                  zdep = MAX( gdept_n(ji,jj,jk) - 1000., 0. ) / 1000.
+                  zahu(ji,jj,jk) = zahu(ji,jj,jk) * MAX( 1., rn_fact_lap * EXP( -zdep ) )
+               ENDIF
+            END DO
+         END DO
+      END DO
+      !
       SELECT CASE ( nldf )                     !* compute lateral mixing trend and add it to the general trend
       !
@@ -136 +148 @@
       NAMELIST/namtrc_ldf/ ln_trcldf_lap, ln_trcldf_blp,                                  &
          &                 ln_trcldf_lev, ln_trcldf_hor, ln_trcldf_iso, ln_trcldf_triad,  &
-         &                 rn_ahtrc_0   , rn_bhtrc_0
+         &                 rn_ahtrc_0   , rn_bhtrc_0, rn_fact_lap  
       !!----------------------------------------------------------------------
       !
@@ -164 +176 @@
          WRITE(numout,*) '           laplacian                 rn_ahtrc_0      = ', rn_ahtrc_0
          WRITE(numout,*) '         bilaplacian                 rn_bhtrc_0      = ', rn_bhtrc_0
+         WRITE(numout,*) '      enhanced zonal diffusivity     rn_fact_lap     = ', rn_fact_lap
+
       ENDIF
       !