Hi, the variable `nstages` is an integer, but you are passing `MPI_DOUBLE_PRECISION` to the MPI call in `fortress_smc_t.f90`: ``` call mpi_bcast(self%temp%nstages, 1, MPI_DOUBLE_PRECISION, & 0, MPI_COMM_WORLD, mpierror) ``` This is clearly a bug.