resume.f90
22b60c541a80geom_testmumps

resume.f90
View Options

			SUBROUTINE resume
			!
			! Resume from previous run
			!
			Use beam
			Use basic
			Use fields
			Use sort
			Use maxwellsource
			Use geometry
			IMPLICIT NONE

			!
			! Local vars and arrays
			INTEGER:: i
			!________________________________________________________________________________
			WRITE(*,'(a)') ' Resume from previous run'
			!________________________________________________________________________________
			!
			! Open and read initial conditions from restart file
			CALL chkrst(0)

			! TODO: change subroutine to take into account multiple species

			IF (newres) THEN
			cstep=0
			time=0
			END IF
			call read_geom(lu_in, rnorm, rgrid, Potinn, Potout)
			! Compute Self consistent electric field
			CALL fields_init
			IF(mpisize .gt. 1) THEN
			CALL distribpartsonprocs(partslist(1))
			DO i=1,nbspecies
			CALL keep_mpi_self_parts(partslist(i),Zbounds)
			END DO
			END IF
			CALL bound(partslist(1))
			CALL localisation(partslist(1))
			! Init Electric and Magnetic Fields
			CALL fields_comm_init
			WRITE(,) "Fields initialized"
			CALL rhscon(partslist(1))
			WRITE(,) "Initial rhs computed"
			CALL poisson(splrz)
			WRITE(,) "Initial field solved"
			CALL EFieldscompatparts(partslist(1))
			WRITE(,) "Initial forces computed"
			IF (newres) THEN
			CALL diagnostics
			WRITE(,) "Initial beam%diagnostics done!"
			END IF
			IF (nlmaxwellsource) CALL maxwellsource_init(lu_in, time)
			!
			END SUBROUTINE resume