Changeset 3563 for trunk

Timestamp:

2012-11-15T18:03:44+01:00 (11 years ago)

Author:

acc

Message:

#831 and #962: quick fix to initialise ghost rows of the envelope bathymetry in zgr_sco (domzgr.F90). This avoids undefined e3 errors later on but does not solve the generic problem of unset ghost rows with some mpp decompositions.

File:

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

trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -1245 +1245 @@
          END DO
          !
-         ! Apply lateral boundary condition   CAUTION: kept the value when the lbc field is zero
+         ! Apply lateral boundary condition   CAUTION: keep the value when the lbc field is zero
          ztmp(:,:) = zenv(:,:)   ;   CALL lbc_lnk( zenv, 'T', 1._wp )
          DO jj = 1, nlcj
@@ -1256 +1256 @@
       !                                                     ! ================ !
       !
-      !                                        ! envelop bathymetry saved in hbatt
+      ! Fill ghost rows with appropriate values to avoid undefined e3 values with some mpp decompositions
+      DO ji = nlci+1, jpi 
+         zenv(ji,1:nlcj) = zenv(nlci,1:nlcj)
+      END DO
+      !
+      DO jj = nlcj+1, jpj
+         zenv(:,jj) = zenv(:,nlcj)
+      END DO
+      !
+      ! Envelope bathymetry saved in hbatt
       hbatt(:,:) = zenv(:,:) 
       IF( MINVAL( gphit(:,:) ) * MAXVAL( gphit(:,:) ) <= 0._wp ) THEN