Changeset 2493

Timestamp:

2010-12-20T14:47:12+01:00 (13 years ago)

Author:

gm

Message:

v3.3beta: #777 backward compatibility with the setting of the minimum ocean depth: set mbathy to 0 over lands

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (2 diffs)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -805 +805 @@
       INTEGER  ::   ji, jj, jk       ! dummy loop indices
       INTEGER  ::   ik, it           ! temporary integers
-      INTEGER  ::   nlbelow          ! temporary integer
       LOGICAL  ::   ll_print         ! Allow  control print for debugging
       REAL(wp) ::   ze3tp , ze3wp    ! Last ocean level thickness at T- and W-points
@@ -833 +832 @@
       ! ===================
       zmax = gdepw_0(jpk) + e3t_0(jpk)          ! maximum depth (i.e. the last ocean level thickness <= 2*e3t_0(jpkm1) )
-      mbathy(:,:) = jpkm1                       ! initialize mbathy to the maximum ocean level available
       bathy(:,:) = MIN( zmax ,  bathy(:,:) )    ! bounded value of bathy (min already set at the end of zgr_bat)
+      WHERE( bathy(:,:) == 0._wp )   ;   mbathy(:,:) = 0       ! land  : set mbathy to 0
+      ELSE WHERE                     ;   mbathy(:,:) = jpkm1   ! ocean : initialize mbathy to the max ocean level
+      END WHERE
 
       ! Compute mbathy for ocean points (i.e. the number of ocean levels)