Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90

Timestamp:

2012-01-28T17:44:18+01:00 (12 years ago)

Author:

rblod

Message:

Merge of 3.4beta into the trunk

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90 (modified) (9 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90

@@ -22 +22 @@
    USE lib_mpp         ! distributed memory computing
    USE prtctl          ! Print control
+   USE timing          ! Timing
 
    IMPLICIT NONE
@@ -36 +37 @@
    REAL(wp) ::   rn_bfrien = 30._wp      ! local factor to enhance coefficient bfri
    LOGICAL  ::   ln_bfr2d  = .false.     ! logical switch for 2D enhancement
-   
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::  bfrcoef2d   ! 2D bottom drag coefficient
+   LOGICAL , PUBLIC                            ::  ln_bfrimp = .false.  ! logical switch for implicit bottom friction
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::  bfrcoef2d            ! 2D bottom drag coefficient
 
    !! * Substitutions
@@ -82 +83 @@
       REAL(wp) ::   zvu, zuv, zecu, zecv   ! temporary scalars
       !!----------------------------------------------------------------------
-
+      !
+      IF( nn_timing == 1 )  CALL timing_start('zdf_bfr')
+      !
       IF( nn_bfr == 2 ) THEN                 ! quadratic botton friction
          ! Calculate and store the quadratic bottom friction coefficient bfrua and bfrva
@@ -120 +123 @@
             &                       tab2d_2=bfrva, clinfo2=       ' v: ', mask2=vmask,ovlap=1 )
       ENDIF
+
+      !
+      IF( nn_timing == 1 )  CALL timing_stop('zdf_bfr')
       !
    END SUBROUTINE zdf_bfr
@@ -142 +148 @@
       REAL(wp) ::  zfru, zfrv         !    -         -
       !!
-      NAMELIST/nambfr/ nn_bfr, rn_bfri1, rn_bfri2, rn_bfeb2, ln_bfr2d, rn_bfrien
-      !!----------------------------------------------------------------------
-
+      NAMELIST/nambfr/ nn_bfr, rn_bfri1, rn_bfri2, rn_bfeb2, ln_bfr2d, rn_bfrien, ln_bfrimp
+      !!----------------------------------------------------------------------
+      !
+      IF( nn_timing == 1 )  CALL timing_start('zdf_bfr_init')
+      !
       REWIND ( numnam )               !* Read Namelist nam_bfr : bottom momentum boundary condition
       READ   ( numnam, nambfr )
@@ -156 +164 @@
       !                              ! allocate zdfbfr arrays
       IF( zdf_bfr_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'zdf_bfr_init : unable to allocate arrays' )
+
+      !                              ! Make sure ln_zdfexp=.false. when use implicit bfr
+      IF( ln_bfrimp .AND. ln_zdfexp ) THEN
+         IF(lwp) THEN
+            WRITE(numout,*)
+            WRITE(numout,*) 'Implicit bottom friction can only be used when ln_zdfexp=.false.'
+            WRITE(numout,*) '         but you set: ln_bfrimp=.true. and ln_zdfexp=.true.!!!!'
+            WRITE(ctmp1,*)  '         set either ln_zdfexp = .false or ln_bfrimp = .false.'
+            CALL ctl_stop( ctmp1 )
+         END IF
+      END IF
 
       SELECT CASE (nn_bfr)
@@ -207 +226 @@
          !
       END SELECT
+      IF(lwp) WRITE(numout,*) '      implicit bottom friction switch                ln_bfrimp  = ', ln_bfrimp
       !
       ! Basic stability check on bottom friction coefficient
@@ -253 +273 @@
          CALL mpp_max( zmaxbfr )
       ENDIF
+      IF( .NOT.ln_bfrimp) THEN
       IF( lwp .AND. ictu + ictv > 0 ) THEN
          WRITE(numout,*) ' Bottom friction stability check failed at ', ictu, ' U-points '
@@ -259 +280 @@
          WRITE(numout,*) ' Bottom friction coefficient will be reduced where necessary'
       ENDIF
+      ENDIF
+      !
+      IF( nn_timing == 1 )  CALL timing_stop('zdf_bfr_init')
       !
    END SUBROUTINE zdf_bfr_init