function sumLL=tjpc(L)
LL=dms2rad(mat2dms(L))
A=[1,1,1];
w=sum(LL(:))-pi;
w=dms2mat(rad2dms(w));
P=eye(3);
Naa=A*inv(P)*A';
Ka=-inv(Naa)*w;
V=A'*Ka;
L1=L+V;
LL=dms2rad(mat2dms(L1));
sumLL=sum(LL);