Changeset 422

Timestamp:

2006-03-21T10:16:38+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_042 : CT : use avt0 defined in the namelist to compute the background eddy diffusivity profil avtb in the case navtb /= 0

File:

• trunk/NEMO/OPA_SRC/ZDF/zdftke.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

@@ -21 +21 @@
    USE taumod          ! surface stress
    USE prtctl          ! Print control
+   USE traadv_ctl      ! advection scheme control
 
    IMPLICIT NONE
@@ -791 +792 @@
          !   Background profile of avt (fit a theoretical/observational profile (Krauss 1990) 
          avmb(:) = avm0
-         avtb(:) = 1.e-5 + 2.8e-8 * gdepw(:)   ! m2/s
+         avtb(:) = avt0 + ( 3.0e-4 - 2 * avt0 ) * 1.0e-4 * gdepw(:)   ! m2/s
       ENDIF
 
-      !   Increase the background in the surface layers
-      avmb(1) = 10.  * avmb(1)      ;      avtb(1) = 10.  * avtb(1)
-      avmb(2) = 10.  * avmb(2)      ;      avtb(2) = 10.  * avtb(2)
-      avmb(3) =  5.  * avmb(3)      ;      avtb(3) =  5.  * avtb(3)
-      avmb(4) =  2.5 * avmb(4)      ;      avtb(4) =  2.5 * avtb(4)
+      IF( cp_cfg == "orca"  .AND. jp_cfg == 2 .AND. ln_traadv_cen2 )  THEN
+         !   Increase the background in the surface layers
+         avmb(1) = 10.  * avmb(1)      ;      avtb(1) = 10.  * avtb(1)
+         avmb(2) = 10.  * avmb(2)      ;      avtb(2) = 10.  * avtb(2)
+         avmb(3) =  5.  * avmb(3)      ;      avtb(3) =  5.  * avtb(3)
+         avmb(4) =  2.5 * avmb(4)      ;      avtb(4) =  2.5 * avtb(4)
+      ENDIF