inital.f90
2e034541f079geom_testmumps

inital.f90
View Options

			SUBROUTINE inital
			!
			USE basic
			USE beam
			USE fields
			USE mpihelper
			USE maxwsrce
			Use geometry
			Use neutcol
			!
			! Set initial conditions
			!
			IMPLICIT NONE
			INTEGER:: i
			!
			!________________________________________________________________________________
			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 init_mag
			CALL load_parts ! will call the localisation

			IF(mpisize .gt. 1) THEN
			CALL calc_Zbounds(partslist(1), Zbounds, femorder)
			Do i=1,size(partslist)
			CALL keep_mpi_self_parts(partslist(i), Zbounds)
			END DO
			END IF

			CALL fields_init
			CALL fields_comm_init(Zbounds)
			WRITE(,) "Fields initialized"
			CALL rhscon(partslist)
			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"





			!________________________________________________________________________________
			END SUBROUTINE inital