boris_update_velocity Subroutine

public pure subroutine boris_update_velocity(v, q, m, dt, e, b, v_new)

電場半ステップ加速と磁場回転を組み合わせ、1タイムステップ後の速度を計算する。

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: v(3)
real(kind=dp), intent(in) :: q

粒子1個あたりの電荷 [C]。
real(kind=dp), intent(in) :: m

粒子1個あたりの電荷 [C]。 粒子1個あたりの質量 [kg]。
real(kind=dp), intent(in) :: dt

粒子1個あたりの電荷 [C]。 粒子1個あたりの質量 [kg]。 時間刻み幅 [s]。
real(kind=dp), intent(in) :: e(3)

粒子1個あたりの電荷 [C]。 粒子1個あたりの質量 [kg]。 時間刻み幅 [s]。
real(kind=dp), intent(in) :: b(3)

粒子1個あたりの電荷 [C]。 粒子1個あたりの質量 [kg]。 時間刻み幅 [s]。
real(kind=dp), intent(out) :: v_new(3)

Calls

proc~~boris_update_velocity~~CallsGraph proc~boris_update_velocity boris_update_velocity proc~cross~2 cross proc~boris_update_velocity->proc~cross~2

Called by

proc~~boris_update_velocity~~CalledByGraph proc~boris_update_velocity boris_update_velocity proc~boris_push boris_push proc~boris_push->proc~boris_update_velocity proc~build_particle_step_candidate build_particle_step_candidate proc~build_particle_step_candidate->proc~boris_push proc~advance_particle_step advance_particle_step proc~advance_particle_step->proc~build_particle_step_candidate proc~advance_particle_step_upper_panel_fourier advance_particle_step_upper_panel_fourier proc~advance_particle_step_upper_panel_fourier->proc~build_particle_step_candidate proc~resolve_particle_boundary_candidate resolve_particle_boundary_candidate proc~resolve_particle_boundary_candidate->proc~build_particle_step_candidate

Source Code

  pure subroutine boris_update_velocity(v, q, m, dt, e, b, v_new)
    real(dp), intent(in) :: v(3), q, m, dt, e(3), b(3)
    real(dp), intent(out) :: v_new(3)
    real(dp) :: qm, v_minus(3), t(3), s(3), v_prime(3), v_plus(3), t2

    qm = q/m
    v_minus = v + qm*e*(0.5d0*dt)
    t = qm*b*(0.5d0*dt)
    t2 = sum(t*t)
    s = 2.0d0*t/(1.0d0 + t2)
    v_prime = v_minus + cross(v_minus, t)
    v_plus = v_minus + cross(v_prime, s)
    v_new = v_plus + qm*e*(0.5d0*dt)
  end subroutine boris_update_velocity