C	PROGRAM URSAMAJOR
C	TRANSFORMS POSITIONAL, PROPER MOTION AND RADIAL VELOCITY DATA
C	FOR URSA MAJOR CLUSTER STARS TO CARTESIAN COORDINATE GRID
	REAL  X,Y,XX,XY,Z,X0,Y0,Z0,SUMVR,SUMV0,AVEVR,AVEV0
	REAL  SEVR,SEV0,DV02,DVR2
	INTEGER  I,J
	DIMENSION  A(50),B(50),BB(50),PA(50),PB(50),VR(50),C1(50),AB(50)
	DIMENSION  C2(50),AA(50),C3(50),C4(50),C5(50),APA(50),APB(50)
	DIMENSION  C6(50),C7(50),C8(50),VRR(50),AAA(50),AAB(50)
	DATA  X0,Y0,Z0,I / 204.2057834,46.56422222,57.2957795131,0 /
	OPEN (UNIT=7, FILE='uma.dat', STATUS='OLD')
	SUMVR = 0.0
10	READ (7,*,END=12) X, Y, XX, XY, Z
	 I = I + 1
	 A(I) = (X - X0)/Z0
	 B(I) = Y/Z0
	 BB(I) = (Y - Y0)/Z0
	 PA(I) = XX
	 PB(I) = XY
	 VR(I) = Z
	 SUMVR = SUMVR + VR(I)
	 GOTO 10
12	AVEVR = SUMVR/REAL(I)
	DVR2 = 0.0
	DO 13 J=1,I
	 C1(J) = SIN(B(J))*COSD(Y0) - COS(B(J))*SIND(Y0)*COS(A(J))
	 AB(J) = ASIN(C1(J))
	 C2(J) = SIN(A(J))*COS(B(J))/COS(AB(J))
	 AA(J) = ASIN(C2(J))
	 C3(J) = A(J)*COS(B(J))*TAND(Y0)
	 C4(J) = 0.5*BB(J)*BB(J)
	 C5(J) = ASIN(C3(J))
	 APA(J) = PA(J)*COS(B(J))*(COS(C5(J)) + C4(J)) - PB(J)*SIN(C5(J))
	 APB(J) = PB(J)*COS(C5(J)) + PA(J)*COS(B(J))*SIN(C5(J))
13	 DVR2 = DVR2 + (VR(J) - AVEVR)**2
	SEVR = SQRT(DVR2/(REAL(I)*REAL(I-1)))
	SUMV0 = 0.0
	DO 14, J=1,I
	 C6(J) = VR(J)/(COS(AA(J))*COS(AB(J)))
	 C7(J) = 23.59*4.74*0.068*TAN(AA(J))
	 C8(J) = 23.59*4.74*(-0.030)*TAN(AB(J))/COS(AA(J))
	 VRR(J) = C6(J) - C7(J) - C8(J)
	 SUMV0 = SUMV0 + VRR(J)
	 AAA(J) = AA(J)*Z0
14	 AAB(J) = AB(J)*Z0
	AVEV0 = SUMV0/REAL(I)
	DV02 = 0.00
	DO 16, J=1,I
	DV02 = DV02 + (VRR(J) - AVEV0)**2
	WRITE (6,15) AAA(J), AAB(J), APA(J), APB(J), VRR(J)
15	FORMAT (1X,5F15.5)
16	CONTINUE
	SEV0 = SQRT(DV02/(REAL(I)*REAL(I-1)))
	WRITE (6,17) AVEVR, SEVR, AVEV0, SEV0
17	FORMAT (2X,4F15.5)
	STOP
	END