;+
;FUNCTION: pangle,theta,phi,b_theta,b_phi
;PURPOSE:
; Computes pitch angle given two sets of theta and phi
;INPUT:
; theta,phi: double (array or scaler) first directions
; b_theta,b_phi : double (array or scaler) second directions
;RETURNS: pitch angle (array or scaler) same dimensions as theta and phi
;
;CREATED BY: Davin Larson
;LAST MODIFICATION: @(#)pangle.pro 1.4 95/11/28
;
;-
function pangle, theta,phi,b_theta,b_phi,vec = vec
if n_elements(vec) eq 3 then xyz_to_polar,vec,theta=b_theta,phi=b_phi
sphere_to_cart,1.,b_theta,b_phi,bx,by,bz
sphere_to_cart,1., theta, phi,sx,sy,sz
dot = (bx*sx + by*sy + bz*sz)
pa = 180./!dpi*acos(dot)
return, pa
end