Changeset 1719 for trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

Timestamp:

2009-11-06T15:01:56+01:00 (14 years ago)

Author:

rblod

Message:

Suppress bottom boundary condition on tke, see ticket #592

File:

• trunk/NEMO/OPA_SRC/ZDF/zdftke.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

@@ -213 +213 @@
       !                     !  Bottom boundary condition on tke
       !                     !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
-      !                     en(bot)   = 0.5*sqrt(u_botfr^2+v_botfr^2) (min value rn_emin)
-!CDIR NOVERRCHK
-      DO jj = 2, jpjm1
-!CDIR NOVERRCHK
-         DO ji = fs_2, fs_jpim1   ! vector opt.
-            ikbu = MIN( mbathy(ji+1,jj), mbathy(ji,jj) )
-            ikbv = MIN( mbathy(ji,jj+1), mbathy(ji,jj) )
-            ikbum1 = MAX( ikbu-1, 1 )
-            ikbvm1 = MAX( ikbv-1, 1 )
-            ikbu = MIN( mbathy(ji,jj), mbathy(ji-1,jj) )
-            ikbv = MIN( mbathy(ji,jj), mbathy(ji,jj-1) )
-            ikbumm1 = MAX( ikbu-1, 1 )
-            ikbvmm1 = MAX( ikbv-1, 1 )
-            ikbt = MAX( mbathy(ji,jj), 1 )
-            ztx2 = bfrua(ji-1,jj) * fse3u(ji-1,jj  ,ikbumm1) + &
-                   bfrua(ji  ,jj) * fse3u(ji  ,jj  ,ikbum1 )
-            zty2 = bfrva(ji,jj  ) * fse3v(ji  ,jj  ,ikbvm1) + &
-                   bfrva(ji,jj-1) * fse3v(ji  ,jj-1,ikbvmm1 )
-            zebot = 0.25_wp * SQRT( ztx2 * ztx2 + zty2 * zty2 )
-            en (ji,jj,ikbt) = MAX( zebot, rn_emin ) * tmask(ji,jj,1)
-         END DO
-      END DO
+      !
+      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+      ! Tests to date have found the bottom boundary condition on tke to have very little effect.
+      ! The condition is coded here for completion but commented out until there is proof that the
+      ! computational cost is justified
+      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+      !
+      !                     en(bot)   = (rn_ebb0/rau0)*0.5*sqrt(u_botfr^2+v_botfr^2) (min value rn_emin)
+!CDIR NOVERRCHK
+!!    DO jj = 2, jpjm1
+!CDIR NOVERRCHK
+!!       DO ji = fs_2, fs_jpim1   ! vector opt.
+!!          ikbu = MIN( mbathy(ji+1,jj), mbathy(ji,jj) )
+!!          ikbv = MIN( mbathy(ji,jj+1), mbathy(ji,jj) )
+!!          ikbum1 = MAX( ikbu-1, 1 )
+!!          ikbvm1 = MAX( ikbv-1, 1 )
+!!          ikbu = MIN( mbathy(ji,jj), mbathy(ji-1,jj) )
+!!          ikbv = MIN( mbathy(ji,jj), mbathy(ji,jj-1) )
+!!          ikbumm1 = MAX( ikbu-1, 1 )
+!!          ikbvmm1 = MAX( ikbv-1, 1 )
+!!          ikbt = MAX( mbathy(ji,jj), 1 )
+!!          ztx2 = bfrua(ji-1,jj) * ub(ji-1,jj,ikbumm1) + &
+!!                 bfrua(ji  ,jj) * ub(ji  ,jj  ,ikbum1 )
+!!          zty2 = bfrva(ji,jj  ) * vb(ji  ,jj  ,ikbvm1) + &
+!!                 bfrva(ji,jj-1) * vb(ji  ,jj-1,ikbvmm1 )
+!!          zebot = 0.001875_wp * SQRT( ztx2 * ztx2 + zty2 * zty2 )   !  where 0.001875 = (rn_ebb0/rau0) * 0.5 = 3.75*0.5/1000.
+!!          en (ji,jj,ikbt) = MAX( zebot, rn_emin ) * tmask(ji,jj,1)
+!!       END DO
+!!    END DO
       !
       !