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