start.f90
d24dd0989d0emaster

start.f90
View Options

			SUBROUTINE start
			!
			! Start or restart the run
			!
			USE basic,ONLY: mpirank, nlmaxwellsource, time, lu_in, Zbounds, nbspecies, cstep, nlres
			USE maxwsrce
			USE neutcol
			USE geometry
			USE beam
			USE fields
			USE mpi
			USE psupply

			IMPLICIT NONE
			INTEGER:: i, nbbounds

			IF(mpirank.eq.0) WRITE(*,'(a/)') '=== Start or restart the run ==='

			! Particles diagnostics
			IF (cstep .eq. 0) THEN
			CALL partdiagnostics
			DO i=1,nbspecies
			IF(.not. partslist(i)%calc_moments) CYCLE
			CALL momentsdiag(partslist(i))
			END DO
			END IF

			! Initialize electron neutral collisions
			CALL neutcol_init(lu_in, partslist(1))

			! Initialize the external power supply
			IF(nlres .and. the_ps%active) THEN
			WRITE(,) "the_ps is active"
			CALL psupply_init(lu_in, cstep, nbbounds, neutdens, the_ps%bias)
			ELSE
			WRITE(,) "use default psupply init"
			CALL psupply_init(lu_in, cstep, nbbounds, neutdens)
			END IF

			! Activate the source
			IF (nlmaxwellsource) CALL maxwsrce_init(lu_in, time, Zbounds)

			END SUBROUTINE start