New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3563 – NEMO

Changeset 3563


Ignore:
Timestamp:
2012-11-15T18:03:44+01:00 (12 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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r3421 r3563  
    12451245         END DO 
    12461246         ! 
    1247          ! Apply lateral boundary condition   CAUTION: kept the value when the lbc field is zero 
     1247         ! Apply lateral boundary condition   CAUTION: keep the value when the lbc field is zero 
    12481248         ztmp(:,:) = zenv(:,:)   ;   CALL lbc_lnk( zenv, 'T', 1._wp ) 
    12491249         DO jj = 1, nlcj 
     
    12561256      !                                                     ! ================ ! 
    12571257      ! 
    1258       !                                        ! envelop bathymetry saved in hbatt 
     1258      ! Fill ghost rows with appropriate values to avoid undefined e3 values with some mpp decompositions 
     1259      DO ji = nlci+1, jpi  
     1260         zenv(ji,1:nlcj) = zenv(nlci,1:nlcj) 
     1261      END DO 
     1262      ! 
     1263      DO jj = nlcj+1, jpj 
     1264         zenv(:,jj) = zenv(:,nlcj) 
     1265      END DO 
     1266      ! 
     1267      ! Envelope bathymetry saved in hbatt 
    12591268      hbatt(:,:) = zenv(:,:)  
    12601269      IF( MINVAL( gphit(:,:) ) * MAXVAL( gphit(:,:) ) <= 0._wp ) THEN 
Note: See TracChangeset for help on using the changeset viewer.