inital.f90
577903136c3egeom_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 fields_init
			call timera(0, "read_geom")
			call read_geom(lu_in, rnorm, splrz, Potinn, Potout)
			call timera(1, "read_geom")
			call init_mag
			CALL load_parts ! will call the localisation

			IF(mpisize .gt. 1) THEN
			CALL calc_Zbounds(partslist(1), Zbounds, femorder)
			END IF

			Do i=1,size(partslist)
			CALL keep_mpi_self_parts(partslist(i), Zbounds)
			END DO

			call fields_start
			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