Homec4science
Diffusion FENNECS lite (3Dmpimode) Browse

fennecslite/src/inital.f90
6e6d693148f13Dmpimode

Last Change
Raw File

inital.f90
View Options

SUBROUTINE inital
!
USE basic
USE beam
USE fields
USE mpihelper
USE maxwsrce
USE neutcol
Use geometry
USE splinebound
!
! Set initial conditions
!
IMPLICIT NONE
INTEGER:: i, m, nbbounds
!
!________________________________________________________________________________
IF(mpirank .eq. 0) WRITE(*,'(a/)') '=== Set initial conditions ==='
!________________________________________________________________________________
!
! Init Electric and Magnetic Fields
ALLOCATE(partslist(nbspecies))
ALLOCATE(partslist_towrite(nbspecies))
CALL fields_init
call timera(0, "read_geom")
call read_geom(lu_in, rnorm, splrz(1), Potinn, Potout)
call timera(1, "read_geom")
call mag_init
! resize nblost array to adapt for correct number of boundaries
nbbounds=2
if(the_domain%nbsplines .gt. 0) nbbounds=the_domain%nbsplines
Do i=1,nbspecies
if( allocated(partslist(i)%nblost)) deallocate(partslist(i)%nblost)
allocate(partslist(i)%nblost(4+nbbounds))
partslist(i)%nblost=0
end do
call fields_start
CALL load_parts ! will call the localisation
qnorm=abs(partslist(1)%weight*partslist(1)%q)
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)
call collectparts(partslist(i),partslist_towrite(i))
END DO
CALL fields_comm_init(Zbounds)
CALL rhscon(partslist)
CALL poisson
!$OMP PARALLEL
CALL Update_phi
!$OMP END PARALLEL
!CALL EFieldscompatparts(partslist(1))
!call comp_mag_p(partslist(1))
!CALL adapt_vinit(partslist(1))
!________________________________________________________________________________
END SUBROUTINE inital
c4science · Help