Changeset 610 for trunk/NEMO/OPA_SRC/lib_mpp.F90

Timestamp:

2007-02-21T10:22:21+01:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_bugfix_016 : CT : when using key_mpp_mpi, fill extra-hallow bands with a computed value, i.e. inside the domain not with value 1.

File:

• trunk/NEMO/OPA_SRC/lib_mpp.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/lib_mpp.F90

@@ -635 +635 @@
 
       !! * Local variables
-      INTEGER ::   ji, jk, jl   ! dummy loop indices
+      INTEGER ::   ji, jj, jk, jl                        ! dummy loop indices
       INTEGER ::   imigr, iihom, ijhom, iloc, ijt, iju   ! temporary integers
       INTEGER ::   ml_req1, ml_req2, ml_err     ! for key_mpi_isend
@@ -645 +645 @@
 
       IF( PRESENT( cd_mpp ) ) THEN
-         ! only fill extra allows with 1.
-         ptab(     1:nlci, nlcj+1:jpj, :) = 1.e0
-         ptab(nlci+1:jpi ,       :   , :) = 1.e0
+         DO jj = nlcj+1, jpj   ! only fill extra allows last line
+            ptab(1:nlci, jj, :) = ptab(1:nlci, nlcj, :)
+         END DO
+         DO ji = nlci+1, jpi   ! only fill extra allows last column
+            ptab(ji    , : , :) = ptab(nlci  , :   , :)
+         END DO
       ELSE      
 
@@ -1095 +1098 @@
       ! ------------------------------
       IF (PRESENT(cd_mpp)) THEN
-         ! only fill extra allows with 1.
-         pt2d(     1:nlci, nlcj+1:jpj) = 1.e0
-         pt2d(nlci+1:jpi ,       :   ) = 1.e0
-     
+         DO jj = nlcj+1, jpj   ! only fill extra allows last line
+            pt2d(1:nlci, jj) = pt2d(1:nlci, nlcj)
+         END DO
+         DO ji = nlci+1, jpi   ! only fill extra allows last column
+            pt2d(ji    , : ) = pt2d(nlci  , :   )
+         END DO     
       ELSE