MODULE domain_param

INTEGER :: nsplit_i
INTEGER :: nsplit_j
INTEGER :: halo=1

INTEGER, PARAMETER :: default_nsplit_i=1
INTEGER, PARAMETER :: default_nsplit_j=1

CONTAINS
  
  SUBROUTINE init_domain_param
  USE ioipsl, ONLY: getin
  USE mpipara, ONLY: mpi_size, is_mpi_master
  IMPLICIT NONE
    nsplit_i=default_nsplit_i
    nsplit_j=default_nsplit_j
    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
  
END MODULE domain_param