stepon.f90
0b252c046024master

stepon.f90
View Options

			!---------------------------------------------------------------------------
			! DESCRIPTION:
			! Advance one time-step
			!---------------------------------------------------------------------------
			SUBROUTINE stepon

			USE basic
			USE constants
			USE fields
			USE beam
			IMPLICIT NONE

			! Compute the positions of the particles at t = n+1/2
			CALL comp_halfpos(parts)

			! Localize and remove the particles that are outside the domain
			CALL boundary_loss(parts)

			! Deposit charge on the grid
			CALL deposit_charge(parts)

			! Compute the self consistent electric field on the grid at t = n+1/2
			CALL poisson_solver

			! Compute the electric field at the particle position at t = n+1/2
			CALL comp_fields_p(parts)

			! Compute the velocity of the particles at t = n+1/2
			CALL comp_velocity(parts)

			! Compute the positions of the particles at t = n+1
			CALL push(parts)

			! Localize and remove the particles that are outside the domain
			CALL boundary_loss(parts)

			! Compute the main physical quantities
			!CALL part_diagnostics

			! Save the data
			CALL diagnose(step)
			END SUBROUTINE stepon