Changeset 292 for codes

Timestamp:

10/22/14 13:19:35 (10 years ago)

Author:

millour

Message:

Added an explicit error message when nsplit_i and nsplit_j values are not compatible with number of MPI processes.
EM

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/domain_param.f90 (modified) (2 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/domain_param.f90

@@ -11 +11 @@
   
   SUBROUTINE init_domain_param
-  USE ioipsl
+  USE ioipsl, ONLY: getin
+  USE mpipara, ONLY: mpi_size, is_mpi_master
   IMPLICIT NONE
     nsplit_i=default_nsplit_i
@@ -17 +18 @@
     CALL getin('nsplit_i',nsplit_i)
     CALL getin('nsplit_j',nsplit_j)
+    ! Check that nsplit_i*nsplit_j*10>=mpi_size
+    IF (nsplit_i*nsplit_j*10<mpi_size) THEN
+      IF (is_mpi_master) THEN
+        WRITE(*,*) "Error: Number of MPI processes must be"
+        WRITE(*,*) "       <= 10*nsplit_i_nsplit_j"
+        WRITE(*,*) " but nsplit_i=",nsplit_i," , nsplit_j=",nsplit_j
+        WRITE(*,*) " while using ",mpi_size," MPI processes"
+      ENDIF
+      STOP
+    ENDIF
   END SUBROUTINE init_domain_param