Changeset 2470 for trunk/NEMO/OPA_SRC/ZDF

Timestamp:

2010-12-08T16:08:16+01:00 (13 years ago)

Author:

cetlod

Message:

v3.2 : bug correction in zdfbrf.F90, dynbfr.F90 & dynspg_ts.F90 see ticket #767

File:

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

trunk/NEMO/OPA_SRC/ZDF/zdfbfr.F90

@@ -131 +131 @@
       USE iom   ! I/O module for ehanced bottom friction file
       !!
-      INTEGER ::   inum         ! logical unit for enhanced bottom friction file
-      INTEGER ::   ji, jj       ! dummy loop indexes
-      INTEGER ::   ikbu, ikbv   ! temporary integers
-      INTEGER ::   ictu, ictv   !    -          -
-      REAL(wp) ::  zminbfr, zmaxbfr   ! temporary scalars
-      REAL(wp) ::  zfru, zfrv         !    -         -
+      INTEGER ::   inum                         ! logical unit for enhanced bottom friction file
+      INTEGER ::   ji, jj                       ! dummy loop indexes
+      INTEGER ::   ikbu, ikbv, ikbum1, ikbvm1   ! temporary integers
+      INTEGER ::   ictu, ictv                   !    -          -
+      REAL(wp) ::  zminbfr, zmaxbfr             ! temporary scalars
+      REAL(wp) ::  zfru, zfrv                   !    -         -
       !!
       NAMELIST/nambfr/ nn_bfr, rn_bfri1, rn_bfri2, rn_bfeb2, ln_bfr2d, rn_bfrien
@@ -220 +220 @@
              ikbu = MIN( mbathy(ji+1,jj  ), mbathy(ji,jj) )
              ikbv = MIN( mbathy(ji  ,jj+1), mbathy(ji,jj) )
-             zfru = 0.5 * fse3u(ji,jj,ikbu) / rdt
-             zfrv = 0.5 * fse3v(ji,jj,ikbv) / rdt
+             ikbum1 = MAX( ikbu-1, 1 )
+             ikbvm1 = MAX( ikbv-1, 1 )
+             zfru = 0.5 * fse3u(ji,jj,ikbum1) / rdt
+             zfrv = 0.5 * fse3v(ji,jj,ikbvm1) / rdt
              IF( ABS( bfrcoef2d(ji,jj) ) > zfru ) THEN
                 IF( ln_ctl ) THEN