only 2 and 3 dimension can be 1

This commit is contained in:
Martin Diehl 2021-04-09 08:25:30 +02:00
parent 1b89032086
commit 5f608ed572
1 changed files with 2 additions and 0 deletions

View File

@ -74,6 +74,8 @@ subroutine discretization_grid_init(restart)
allocate(materialAt_global(0)) ! needed for IntelMPI
endif
if (grid(1) < 2) call IO_error(844, ext_msg='cells(1) must be larger than 1')
call MPI_Bcast(grid,3,MPI_INTEGER,0,PETSC_COMM_WORLD, ierr)
if (ierr /= 0) error stop 'MPI error'
call MPI_Bcast(geomSize,3,MPI_DOUBLE,0,PETSC_COMM_WORLD, ierr)