Changeset 2462 for branches/DEV_r2460_v3_3beta_NOL/NEMOGCM/NEMO/OPA_SRC/LBC/mppini_2.h90

Timestamp:

2010-12-07T18:56:19+01:00 (13 years ago)

Author:

acc

Message:

DEV_r2460_v3_3beta_NOL. Changes to v3.3beta to improve assignment of ocean processes to ioserver processes when discarding land-only regions. See ticket #776

File:

• branches/DEV_r2460_v3_3beta_NOL/NEMOGCM/NEMO/OPA_SRC/LBC/mppini_2.h90 (modified) (3 diffs)

branches/DEV_r2460_v3_3beta_NOL/NEMOGCM/NEMO/OPA_SRC/LBC/mppini_2.h90

@@ -41 +41 @@
       USE in_out_manager  ! I/O Manager
       USE iom
+#if defined key_iomput
+      USE mod_event_client
+# endif
       !! 
       INTEGER :: ji, jj, jn, jproc, jarea     ! dummy loop indices
@@ -288 +291 @@
             iin(icont+1) = ii
             ijn(icont+1) = ij
+            nregproc_out(icont+1) = jarea - 1
          ENDIF
       END DO
@@ -461 +465 @@
       IF (lwp) THEN
          CALL ctl_opn( inum, 'layout.dat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE., narea )
-         WRITE(inum,'(6i8)') jpnij,jpi,jpj,jpk,jpiglo,jpjglo
-         WRITE(inum,'(a)') 'NAREA nlci nlcj nldi nldj nlei nlej nimpp njmpp'
+         WRITE(inum,'(8i8)') jpnij,jpi,jpj,jpk,jpiglo,jpjglo,jpni,jpnj
+         WRITE(inum,'(a)') 'NAREA nlci nlcj nldi nldj nlei nlej nimpp njmpp nregproc_out'
 
         DO  jproc = 1, jpnij
-         WRITE(inum,'(9i5)') jproc, nlcit(jproc), nlcjt(jproc), &
+         WRITE(inum,'(10i5)') jproc, nlcit(jproc), nlcjt(jproc), &
                                       nldit(jproc), nldjt(jproc), &
                                       nleit(jproc), nlejt(jproc), &
-                                      nimppt(jproc), njmppt(jproc)
+                                      nimppt(jproc), njmppt(jproc), nregproc_out(jproc)
         END DO
         CLOSE(inum)   
+#ifdef key_iomput
+     ! Check nregproc_in matches nregproc_out and abort if it does not.
+     ! This should only occur on the first run of a new domain decomposition
+     ! when jpni*jpnj /= jpnij. If ioservers are not being used, Init_parallel
+     ! will not have been called and nregproc_in will not have been allocated.
+     ! In such a case the check is unnecessary so only check for a match if
+     ! nregproc_in has been allocated.
+       IF (ALLOCATED(nregproc_in)) THEN
+        DO  jn = 1, jpnij
+         IF ( nregproc_in(jn) /= nregproc_out(jn) ) THEN
+          WRITE(numout,*) ' nregproc_in and nregproc_out do not match.'
+          WRITE(numout,*) ' This is expected when starting a new domain docomposition.'
+          WRITE(numout,*) ' The layout.dat file has been updated and the model may be rerun.'
+          CALL ctl_stop( ' mpp_init: IO server domain colour assignment mismatch ')
+         END IF
+        END DO
+       END IF
+#endif
       END IF