inital.f90
22b60c541a80geom_testmumps

inital.f90
View Options

			SUBROUTINE inital
			!
			USE basic
			USE beam
			USE fields
			USE mpihelper
			USE maxwellsource
			Use geometry
			!
			! Set initial conditions
			!
			IMPLICIT NONE
			!
			!________________________________________________________________________________
			IF(mpirank .eq. 0) WRITE(*,'(a/)') '=== Set initial conditions ==='
			!________________________________________________________________________________
			!
			! Init Electric and Magnetic Fields
			ALLOCATE(partslist(nbspecies))
			call read_geom(lu_in, rnorm, rgrid, Potinn, Potout)
			CALL load_parts ! will call the localisation
			CALL fields_init
			CALL fields_comm_init
			WRITE(,) "Fields initialized"
			IF (nlmaxwellsource) CALL maxwellsource_init(lu_in, time)
			CALL rhscon(partslist(1))
			WRITE(,) "Initial rhs computed"
			CALL poisson(splrz)
			WRITE(,) "Initial field solved"
			CALL EFieldscompatparts(partslist(1))
			WRITE(,) "Initial forces computed"
			CALL adapt_vinit(partslist(1))
			WRITE(,)"init velocity adapted"
			CALL diagnostics
			WRITE(,) "Initial beam%diagnostics calculated"

			!________________________________________________________________________________
			END SUBROUTINE inital