bvectortgrid Subroutine

private pure subroutine bvectortgrid(invR, nth, gradrptot, dpsidrp, bi, bvector)

FIXME : Add documentation

Arguments

Type IntentOptional Attributes Name
real, intent(in), dimension(-ntgrid: ) :: invR
integer, intent(in) :: nth
real, intent(in), dimension(-ntgrid:,:) :: gradrptot
real, intent(in) :: dpsidrp
real, intent(in) :: bi
real, intent(out), dimension(-ntgrid:,:) :: bvector

Contents

Source Code

bvectortgrid

Source Code

  pure subroutine bvectortgrid(invR, nth, gradrptot, dpsidrp, bi, bvector)
    integer, intent (in) :: nth
    real, dimension(-ntgrid:  ), intent (in) :: invR
    real, dimension(-ntgrid:,:), intent (in) :: gradrptot
    real, dimension(-ntgrid:,:), intent (out) :: bvector
    real, intent (in) :: dpsidrp, bi
    integer :: i
    do i = -nth, nth
       bvector(i,1) = -dpsidrp * gradrptot(i,2) * invR(i)
       bvector(i,2) = dpsidrp * gradrptot(i,1) * invR(i)
       bvector(i,3) = bi * invR(i)
    end do
  end subroutine bvectortgrid