get_unit_vector_spherical Subroutine

public subroutine get_unit_vector_spherical(x_sph, unitv_r, unitv_theta, unitv_phi)

Arguments

Type IntentOptional Attributes Name
double precision, intent(in) :: x_sph(1:3)
double precision, intent(inout) :: unitv_r(1:3)
double precision, intent(inout) :: unitv_theta(1:3)
double precision, intent(inout) :: unitv_phi(1:3)

Calls

proc~~get_unit_vector_spherical~~CallsGraph proc~get_unit_vector_spherical get_unit_vector_spherical dcos dcos proc~get_unit_vector_spherical->dcos dsin dsin proc~get_unit_vector_spherical->dsin

Called by

proc~~get_unit_vector_spherical~~CalledByGraph proc~get_unit_vector_spherical get_unit_vector_spherical proc~integrate_emission_spherical integrate_emission_spherical proc~integrate_emission_spherical->proc~get_unit_vector_spherical proc~integrate_whitelight_spherical integrate_whitelight_spherical proc~integrate_whitelight_spherical->proc~get_unit_vector_spherical proc~get_image get_image proc~get_image->proc~integrate_emission_spherical proc~get_image->proc~integrate_whitelight_spherical proc~get_euv_image get_EUV_image proc~get_euv_image->proc~get_image proc~get_sxr_image get_SXR_image proc~get_sxr_image->proc~get_image proc~get_whitelight_image get_whitelight_image proc~get_whitelight_image->proc~get_image

Source Code

    subroutine get_unit_vector_spherical(x_sph,unitv_r,unitv_theta,unitv_phi)
      double precision, intent(in) :: x_sph(1:3)
      double precision, intent(inout) :: unitv_r(1:3),unitv_theta(1:3),&
         unitv_phi(1:3)

      unitv_r(1)=dsin(x_sph(2))*dcos(x_sph(3))
      unitv_r(2)=dsin(x_sph(2))*dsin(x_sph(3))
      unitv_r(3)=dcos(x_sph(2))
      unitv_theta(1)=dcos(x_sph(2))*dcos(x_sph(3))
      unitv_theta(2)=dcos(x_sph(2))*dsin(x_sph(3))
      unitv_theta(3)=-dsin(x_sph(2))
      unitv_phi(1)=-dsin(x_sph(3))
      unitv_phi(2)=dcos(x_sph(3))
      unitv_phi(3)=zero

    end subroutine get_unit_vector_spherical