leapfrog_final Subroutine

public subroutine leapfrog_final(v0, a0, dt, n)

速度后退半步长，使其与位移、加速度同步

Arguments

Type IntentOptional Attributes Name
real(kind=rk), intent(inout), dimension(*) :: v0
real(kind=rk), intent(in), dimension(*) :: a0
real(kind=rk), intent(in) :: dt
integer, intent(in) :: n

Contents

Source Code

leapfrog_final

Source Code

    subroutine leapfrog_final(v0, a0, dt, n)
        real(kind=rk), intent(inout), dimension(*) :: v0
        real(kind=rk), intent(in), dimension(*) :: a0
        real(kind=rk), intent(in) :: dt
        integer, intent(in) :: n

        v0(:n) = v0(:n) - a0(:n)*dt/2

    end subroutine leapfrog_final