Changeset 3820 for branches/2012

Timestamp:

2013-02-21T18:58:13+01:00 (11 years ago)

Author:

hliu

Message:

add the 2D-enhancement for loglayer bfr coefficient. see ticket 1058

File:

• branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90 (modified) (4 diffs)

branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90

@@ -87 +87 @@
       INTEGER  ::   ji, jj                       ! dummy loop indices
       INTEGER  ::   ikbu, ikbv                   ! local integers
-      REAL(wp) ::   zvu, zuv, zecu, zecv, ztmp   ! temporary scalars
+      REAL(wp) ::   zvu, zuv, zecu, zecv         ! temporary scalars
+      REAL(wp) ::   ztmp, ztmp1                  ! temporary scalars
       !!----------------------------------------------------------------------
       !
@@ -100 +101 @@
 
          IF(ln_loglayer) THEN       ! "log layer" bottom friction coefficient
+
+           ! add 2D-enhancement bottom friction
+           ztmp1 = 1._wp
+           IF(ABS(rn_bfri2) >= 1.e-10 ) THEN
+             ztmp1 = 1._wp / rn_bfri2
+           ELSE
+             CALL ctl_stop( 'rn_bfri2 must not be less than 1.e-10')
+           END IF
+
 #  if defined key_vectopt_loop
            DO jj = 1, 1
@@ -108 +118 @@
 #  endif
                 ztmp = 0.5_wp * fse3t(ji,jj,mbkt(ji,jj))
-                ztmp = max(ztmp, rn_bfrz0)
-                bfrcoef2d(ji,jj) = ( log( ztmp / rn_bfrz0 ) / vkarmn ) ** (-2)
+                ztmp = max(ztmp, rn_bfrz0 + 1.e-10)
+                bfrcoef2d(ji,jj) = bfrcoef2d(ji,jj) * ztmp1 * &
+                                 & ( log( ztmp / rn_bfrz0 ) / vkarmn ) ** (-2)
              END DO
            END DO
@@ -251 +262 @@
          !
       END SELECT
+
+      IF( nn_bfr /= 2 .AND. ln_loglayer ) THEN
+         IF(lwp) THEN
+            WRITE(numout,*)
+            WRITE(numout,*) 'Loglayer can only be by applied for quadratic bottom friction' 
+            WRITE(numout,*) 'but you have set: nn_bfr /= 2 and ln_loglayer=.true.!!!!'
+            WRITE(ctmp1,*)  'check nn_bfr and ln_loglayer (should be 2 and true)'
+            CALL ctl_stop( ctmp1 )
+         END IF
+      END IF
+
+
+
       IF(lwp) WRITE(numout,*) '      implicit bottom friction switch                ln_bfrimp  = ', ln_bfrimp
       !