ppinit.f90
cd1286e2d410improved_sugama

ppinit.f90
View Options

			SUBROUTINE ppinit
			!
			! Parallel envrionment. Get the number of requested MPI processes
			!
			!
			USE prec_const
			USE basic
			USE basic_mpi
			!
			IMPLICIT NONE
			!
			INTEGER :: nargs
			INTEGER, PARAMETER :: ndims =3
			CHARACTER(len=32) :: str
			INTEGER,DIMENSION(ndims) :: dims
			!
			INTEGER :: i
			!________________________________________
			!
			! Setup parallel envrionement

			CALL MPI_INIT(ierr)
			CALL MPI_COMM_SIZE(MPI_COMM_WORLD,nprocs,ierr)
			!
			! Get myrank
			CALL MPI_COMM_RANK(MPI_COMM_WORLD,me,ierr)
			!
			WRITE(me_str,'(I3.3)') me
			!
			!
			! check the requested number of proccesses
			!
			nargs = command_argument_count()
			!
			IF( nargs .NE. 0 .AND. nargs .NE. ndims ) THEN
			WRITE(*, '(a,i4,a)') 'Number of arguments not equal to NDIMS =', ndims, '!'
			CALL mpi_abort(MPI_COMM_WORLD, -1, ierr)
			END IF
			!
			IF( nargs .NE. 0) THEN
			DO i=1,nargs
			CALL get_command_argument(i,str,ierr)
			READ(str(1:3),'(i3)') dims(i)
			ENDDO
			! number of processors per collision, i.e. electron-electron electron-ion ion-electron
			Cenproc = dims(1)
			Cinproc = dims(2)
			Caanproc = dims(3)
			ENDIF
			!
			IF(nargs .ne. 0) THEN
			IF( sum(dims) .ne. nprocs) THEN
			IF(me .eq. 0) WRITE(,) 'Number of processors per coll. operator do not match ',nprocs,': # of allocated MPI processes = ',sum(dims)
			CALL MPI_Abort(MPI_COMM_WORLD,-2,ierr)
			ENDIF
			ENDIF
			!
			IF(me .eq. 0) WRITE(,) 'Number of processors =',nprocs
			!
			!
			END SUBROUTINE ppinit