Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (geom_test) Browse

fennecs/src/resume.f90
5fc7f41985a6geom_test

Last Change
Raw File

resume.f90
View Options

SUBROUTINE resume
!
! Resume from previous run
!
Use beam
Use basic
Use fields
Use sort
Use maxwsrce
Use geometry
Use neutcol
IMPLICIT NONE
!
! Local vars and arrays
INTEGER:: i
!________________________________________________________________________________
WRITE(*,'(a)') ' Resume from previous run'
!________________________________________________________________________________
!
! Open and read initial conditions from restart file
CALL chkrst(0)
! If we want to start a new run with a new file
IF (newres) THEN
! we start time and step from 0
cstep=0
time=0
END IF
call read_geom(lu_in, rnorm, rgrid, Potinn, Potout)
! Initialize the magnetic field and the electric field solver
call init_mag
CALL fields_init
call localisation(partslist(1))
IF(mpisize .gt. 1) THEN
CALL calc_Zbounds(partslist(1), Zbounds, femorder)
DO i=1,nbspecies
CALL keep_mpi_self_parts(partslist(i),Zbounds)
call collectparts(partslist(i))
END DO
END IF
! Add the requested test particles and read them from input files
IF (nbaddtestspecies .gt. 0) THEN
Do i=1,nbaddtestspecies
CALL load_part_file(partslist(nbspecies+i),addedtestspecfile(i))
END DO
nbspecies=nbspecies+nbaddtestspecies
END IF
CALL bound(partslist(1))
CALL localisation(partslist(1))
! Allocate the variables for MPI communications
CALL fields_comm_init(Zbounds)
! Solve Poisson for the initial conditions
CALL rhscon(partslist)
CALL poisson(splrz)
! Compute the fields at the particles positions
CALL EFieldscompatparts(partslist(1))
if(mpirank .eq. 0) WRITE(*,*) "Initial forces computed"
!
END SUBROUTINE resume
c4science · Help