Grid sizes larger than ndim^3=2,147,483,647 (4-bytes) not allowed

[bgriffen]

When I compile the code with a grid size of 1024 I get the following compile error:

mpif90  -c  -fopenmp -DMY_OPENMP  -DMPI  mat_ini_cubep3m.F90
mat_ini_cubep3m.F90:548.55:

       call MPI_BCAST(temperature_grid,mesh(1)*mesh(2)*mesh(3)*3,MPI_REAL,0,...

As Garrelt has noted via private communication this is because ndim^3 is above the 4-byte maximum allowed. Since my runs have logical,parameter :: isothermal=.true. I can comment out:

!#ifdef MPI       
!       ! Distribute the input parameters to the other nodes
!       call MPI_BCAST(temperature_grid,mesh(1)*mesh(2)*mesh(3)*3,MPI_REAL,0,&
!            MPI_COMM_NEW,mympierror)
!#endif

It now compiles but needs to be fixed for future.