a = [15600;11760;11610;15170];
b = [7540;2750;14630;610];
c = [10380;16190;7680;13320];
rho = 26570;
phi = asin(c/rho);
theta = asin((b./cos(phi))/rho);
R=(a.^2+b.^2+(c-6370).^2).^(1/2);
t = .0001 + R./299792.458;
tPhiRef = rand()*pi/2;
tightPhi = rand(1,4)*.05*tPhiRef;
tThetaRef = rand()*2*pi;
tightTheta = rand(1,4)*.05*tThetaRef;
tIghta = rho.*cos(tightPhi).*cos(tightTheta);
tIghtb = rho.*cos(tightPhi).*sin(tightTheta);
tIghtc = rho.*sin(tightPhi);
tightR=((tIghta+19.053).^2+(tIghtb-11.318).^2+(tIghtc-6370).^2).^(1/2);
tIghtt = .0001 + tightR./299792.458;