Create initial particles
subroutine particles_create() use mod_global_parameters select case(physics_type_particles) case('advect') call advect_create_particles() case('Lorentz') ! The Vay mover can use the same routine call Lorentz_create_particles() case('GCA') call gca_create_particles() case('sample') call sample_create_particles() case default if (mype == 0) then print *, "Unknown physics_type_particles", & trim(physics_type_particles) print *, "Options are: advect, Lorentz, GCA, sample" call mpistop("Unknown physics_type_particles") end if end select end subroutine particles_create