Changeset 1268

Timestamp:

2009-01-14T15:08:06+01:00 (15 years ago)

Author:

ctlod

Message:

include the missing Von Karman constant in surface mixing length computation and set to 0.1 the default value for the rn_lmin parameter, see ticket: #291

Location:

trunk/NEMO/OPA_SRC/ZDF

Files:

• zdftke.F90 (modified) (2 diffs)
• zdftke2.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

@@ -85 +85 @@
    INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1/2)
    REAL(wp) ::   rn_lmin0 = 0.4_wp          ! surface  min value of mixing length
-   REAL(wp) ::   rn_lmin  = 0.4_wp          ! interior min value of mixing length
+   REAL(wp) ::   rn_lmin  = 0.1_wp          ! interior min value of mixing length
    REAL(wp) ::   rn_efr   = 1.0_wp          ! fraction of TKE surface value which penetrates in the ocean
    REAL(wp) ::   rn_lc    = 0.15_wp         ! coef to compute vertical velocity of Langmuir cells
@@ -200 +200 @@
       ! Buoyancy length scale: l=sqrt(2*e/n**2)
       ! ---------------------
-      IF( ln_mxl0 ) THEN         ! surface mixing length = F(stress) : l=2.e5*sqrt(utau^2 + vtau^2)/(rau0*g)
+      IF( ln_mxl0 ) THEN         ! surface mixing length = F(stress) : l=vkarmn*2.e5*sqrt(utau^2 + vtau^2)/(rau0*g)
 !!gm  this should be useless
          zmxlm(:,:,1) = 0.e0
 !!gm end
-         zraug = 0.5 * 2.e5 / ( rau0 * grav )
+         zraug = 0.5 * vkarmn * 2.e5 / ( rau0 * grav )
          DO jj = 2, jpjm1
 !CDIR NOVERRCHK

trunk/NEMO/OPA_SRC/ZDF/zdftke2.F90

@@ -83 +83 @@
    INTEGER  ::   nn_htau  = 0               ! type of tke profile of penetration (=0/1/2)
    REAL(wp) ::   rn_lmin0 = 0.4_wp          ! surface  min value of mixing length
-   REAL(wp) ::   rn_lmin  = 0.4_wp          ! interior min value of mixing length
+   REAL(wp) ::   rn_lmin  = 0.1_wp          ! interior min value of mixing length
    REAL(wp) ::   rn_efr   = 1.0_wp          ! fraction of TKE surface value which penetrates in the ocean
    REAL(wp) ::   rn_lc    = 0.15_wp         ! coef to compute vertical velocity of Langmuir cells
@@ -451 +451 @@
       ! Buoyancy length scale: l=sqrt(2*e/n**2)
       ! ---------------------
-      IF( ln_mxl0 ) THEN         ! surface mixing length = F(stress) : l=2.e5*sqrt(utau^2 + vtau^2)/(rau0*g)
+      IF( ln_mxl0 ) THEN         ! surface mixing length = F(stress) : l=vkarmn*2.e5*sqrt(utau^2 + vtau^2)/(rau0*g)
 !!gm  this should be useless
          zmxlm(:,:,1) = 0.e0
 !!gm end
-         zraug = 0.5 * 2.e5 / ( rau0 * grav )
+         zraug = 0.5 * vkarmn * 2.e5 / ( rau0 * grav )
          DO jj = 2, jpjm1
 !CDIR NOVERRCHK