- - . ( 150 ) , . GOSTRF.com - 1 . - . . ! , - !

  


||   ||   ||
- :
 


.


 

50-25645.325-89

 

 

1990

 

.



50-25645.325-89

01.07.90

() () , .

.

, , 1.

2 , .

3 -6.

4 , .

1.

1.1. :

0;

t;

t;

.

1.2. 0 FK5: 2000 ., 1, 12h UT1 ( 1,5). T0, J 2000.0, JD (T0) = 2451545.0.

t , .

1.3. T0 - 00z0, .

0 0.

Oz0 T0 .

0 00z0 .

1.4. t - ttzt, .

t t.

Ozt t .

Oyt ttzt .

1.5. t - z, .

t : 1, , T0 t t.

z t .

z .

1.6. - OXYZ, , .

OZ .

- .

OY OXYZ .

2.

2.1. : 2 .

() S (Scp) () : 3: 4 .

: 1 .

UT1.

2.2. S t

Scp = 6h41m50s,54841 + 236s,555367908 × d + 86400s × M + 0s,093104 × t2 -

- 6s,2 × 10-6 × t3 = 100,460618374 + 0,9856473662 × d + 360o × M +

+ 3,879333 × 10-4 × t2 - 2,583 × 10-8 × t3 = 1,7533685592 + 0,0172027918051 × d +

+ 6,2831853072 × M + 6,7707139 × 10-6 × t2 - 4,50876 × 10-10 × t3, (1)

h - ;

m - ;

s - ;

d - T0 t ,

d = JD(t) - 2451545,0; (2)

- UT1 , ;

t - 0 t 36525 ,

(3)

2.3. S

S = S + Na, (4)

Na - t.

2.4. S t

(5)

2.5. , UTC, AUT1 UTC UT1:

UT1 = UTC + DUT1, (6)

.

DUT1 . , .

3.

3.1. , 0 t , :

(7)

(8)

t, t, zt - , t;

0, 0, z0 - , T0;

N - t;

- , : 5( T0) : 6( t).

.

N 2.

3.2. , t , :

(9)

(10)

X, Y, Z - , ;

Rpol - , t ;

Rs - .

Rs Rpol 2.

3.3. , 0 t , :

(11)

(12)

, , z - , ;

Rm - .

Rm 2.

3.4. , t t , :

(13)

(14)

3.5. , t , :

(15)

(16)

- t Oz , .

2.

3.6. , 0 , :

(17)

(18)

(19)

(20)

1

, ,

*

1.

, -

2. UT1

,

3. UTC

, 31 () 30 , UT1 0,9 :

çUTC - UT1ç £ 0,s9

4.

, .

.

5.

, .

.

6. () : 3: 4

() ,

7.

, , (- ) ( )

8.

, ,

9. JD (t)

, 1 4713 . .. t.

. , 12h UT1 0- л

10.

1900 - 1905 ., ,

* . 15855-77.

2

,

1.

(21)

P11 = coszA × cosZA × cosqA - sinzA × sinZA;

P12 = -sinzA × cosZA × cosqA - coszA × sinZA;

P13 = -cosZA × sinqA;

P21 = coszA × sinZA × cosqA + sinzA × cosZA;

P22 = -sinzA × sinZA × cosqA + coszA × cosZA;

P23 = -sinZA × sinqA;

P31 = coszA × sinqA;

P32 = -sinzA × sinqA;

P33 = cosqA.

zA, ZA, qA :

zA = 2306²,2181 × t + 0²,30188 × t2 + 0²,017998 × t3 = 0,111808609 × 10-1 × t + + 0,146356 × 10-5 × t2 + 0,872 × 10-7 × t3;

ZA = 2306²,2181 × t + 1²,09468 × t2 + 0²,018203 × t3 = 0,111808609 × 10-1 × t + + 0,53072 × 10-5 × t2 + 0,883 × 10-7 × t3;

qA = 2004²,3109 × t - 0²,42665 × t2 - 0²,041833 × t3 = 0,97171735 × 10-2 × t - 0,20685 × 10-5 × t2 - 0,2028 × 10-6 × t3.

J 2000,0 t .

2. N

(22)

n11 = cosNf;

n12 = -sinNf × cose0;

n13 = -sinNf × sine0;

n21 = sinNf × cose;

n22 = cosNf × cose × cose0 + sine × sine0;

n23 = cosNf × cose × sine0 - sine × cose0;

n31 = sinNf × sine;

n32 = cosNf × sine × cose0 - cose × sine0;

n33 = cosNf × sine × sine0 + cose × cose0;

Nf - ;

e - ,

e = e0 + Ne;

e0 - ,

e0 = 84381²,448 - 46²,815 × t - 0²,00059 × t2 + 0²,001813 × t3 = 0,4090928042 - 0,2269655 × 10-3 × t - 0,29 × 10-8 × t2 + 0,88 × 10-8 × t3;

Ne - .

Nf Ne :

, , , - , Nf Ne;

, , , , - l, l¢, F, D, W ;

l - ;

l¢ - ;

F - ;

D - ;

W - .

:

l = 485866²,733 + 1717915922²,633 × t + 31²,310 × t2 + 0²,064 × t3 = 2,355548393 + 8328,69142288 × t + 1,517952 × 10-4 × t2 + 3,103 × 10-7 × t3;

l¢ = 1287099²,804 + 129596581²,224 × t - 0²,577 × t2 - 0²,012 × t3 = 6,24003594 + 628,30195602 × t - 2,7974 × 10-6 × t2 - 5,82 × 10-8 × t3

F = 335778²,877 + 1739527263²,137 × t - 13²,257 × t2 + 0²,011 × t3 = 1,62790193 + 8433,46615831 × t - 6,42717 × 10-5 × t2 + 5,33 × 10-8 × t3;

D = 1072261²,307 + 1602961601²,328 × t - 6²,891 × t2 + 0²,019 × t3 = 5,19846951 + 7771,37714617 × t - 3,34085 × 10-5 × t2 + 9,21 × 10-8 × t3;

W = 450160²,280 - 6962890²,539 × t + 7²,455 × t2 + 0²,008 × t3 = 2,182438624 - 33,757045936 × t + 3,61429 × 10-5 × t2 + 3,88 × 10-8 × t3.

, , , , , , , (i = 1, 2, 3, ..., 106) л ( 1986 .).

0²,1 ( 10-6) Nf Ne :

Nf = -17²,1996 × sinW + 0²,2062 × sin2W - 1²,3187 × sin2(F - D + W) +

+ 0²,1426 × sinl¢ - 0²,2274 × sin2(F + W) = -0,83386 × 10-4 × sinW +

+ 0,9997 × 10-6 × sin2W - 0,63932 × 10-5 - sin2(F - D + W) +

+ 0,6913 × 10-6 × sinl¢ - 0,11024 × 10-5 × sin2(F + W); (23)

Ne = 9²,2025 × cosW + 0²,5736 × cos2(F - D + W) + 0²,0977 × cos2(F + D) =

= 0,44615 × 10-4 × cosW + 0,27809 × 10-5 × cos2(F - D + W) +

+ 0,474 × 10-6 × cos2(F + W). (24)

10-8

(25)

3. Rs

(26)

S - t.

4.

(27)

- t.

5. Rpol

(28)

, - , .

t Opxpypzp, .

Opxp .

Opyp , 90 .

xp, yp .

6.

(29)

m - 0 t.

m

m = zA + ZA = 4612²,4362 × t + 1²,39656 × t2 + 0²,036201 × t3 = 0,2236172 × 10-1 × t + 0,67707 × 10-5 × t2 + 0,1755 × 10-6 × t3.

7. Na

Na = Nf × cose. (30)

3

-6

1. JULE (NDMG, T, JD, DJD, DM, D) 1900 2099 ., UT1 J2000.0 .

: NDMG - , , , ;

, 6 1988 . 06051988;

- UT1 , ,

;

, 13 35 12,79 133512,79.

: JD - ;

DJD - ;

D - J2000.0

;

DM - UT1 ,

.

2. SCP (D, DM, SC) .

: D - J2000.0

;

DM - UT1 ,

.

: SC - , .

3. SMOD (D, DM, SM) .

: D - J2000.0

;

DM - UT1 ,

.

: SM - ,

.

4. SICT (D, DM, NP1, SI) .

: D - J2000.0

;

DM - UT1 ,

.

NP1 - ,

0 1:

NP1 = 0 -

.

NP1 = 1 -

0².1 ( 1. - 6).

: SI - , .

5. (D, NP1, HPCI, , , ) , .

: D - J2000.0

;

NP1 - ,

0 1:

NP1 = 0 -

;

NP1 = 1 -

0².1 ( 1. - 6).

: HPCI - , ;

- , ;

- , ;

- , .

6. (D, NP1, NP2, ) .

: D - J2000.0

;

NP1 - ,

0 1:

NP1 = 0 -

;

NP1 = 1 -

0².1 ( 1. - 6).

NP2 - , 0 1:

NP2 = 0 - ;

NP2 = 1 - .

: - .

7. PREMAT (D,PM) .

: D - J2000.0

;

: - .

8. RMUMAT (D, NP1, RM) .

: D - J2000.0

;

NP1 - ,

0 1:

NP1 = 0 -

;

NP1 = 1 -

0².1 ( 1. - 6).

: RM - .

9.

SUBROUTINE JULE (NDMG, T, JD, DJD, DM, D)

NG = NDMG - (NDMG/10000)*10000

N = NDMG/1000000

NM = (NDMG - NG - N*1000000)/10000

JT = T/10000

JM = T/100 - JT*100

JC = T/100.

JC = JC*100

C = T - J

DM = JT/24. + JM/1440. + C/86400.

K = NM - 1

IF (K.EQ.0) GO TO 1

DO 1 I = 1, K

GO TO (2, 3, 2, 4, 2, 4, 2, 2, 4, 2, 4), I

2 N = N + 31

GO TO 1

3 N = N + 28

GO TO 1

4 N = N + 30

1 CONTINUE

L = NG - 1900

KK = L/4

J = KK*4

IF (NM.GT.2.AND.J.EQ.L.AND.L.NE.0) N = N + 1

J = (L - 1)/4

JD = N + J + L*365 + 2415019

IF (DM.LT.0.5) GO TO 10

DJD = DM - 0.5

JD = JD + 1

GO TO 11

10 DJD = DM + 0.5

11 CONTINUE

D = JD + DJD - 2451545.0

RETURN

END

C ******************************************************

SUBROUTINE SCP (D, DM, SC)

= D/36525

SC = 1.7533685592 + 1.72027918051*1.E - 2*D + 6.28318530718*DM

*+ 6.77071394*1. - 6** - 4.50876723*1.E - 10***

KSC = SC/6.28318530718

IF (KSC.LT.O) KSC = KSC - 1

SC = SC - KSC*6.28318530718

RETURN

END

C ********************************************************

SUBROUTINE SMOD (D, DM, SM)

= D/36525.

SM = 1.7533685592 + 0.01720I217957*D + 6.28318530718*DM

*-1.75958*1. - 7***

KSM = SM/6.28318530718

IF (KSM.LT.0) KSM = KSM - 1

SM = SM - KSM*6.28318530718

RETURN

END

C *****************************************************

SUBROUTINE SICT (D, DM, NP1, SI)

CALL HYTAC (D, NP1, HP, HE, , 0)

HALFA = HP*COS (E)

CALL SCP (D, DM, SC)

SI = SC + HALFA

RETURN

END

C *****************************************************

SUBBOUTINE (D, NP1, HPCI, , , 0)

DIMENSION KL (106), KL1 (106), KF (106), KD (106), (106)

DIMENSION APCI (106), BPCI (106), (106), (106)

DATA ((KL (I), I = 1, 22) = 0, 0, -2, 2, -2, 1, 0, 2, 0, 0, 0, 0, 0, 2,

*0, 0, 0, 0, 0, -2, 0, 2)

DATA ((KL(I), I = 23, 44) = 0, 1, 2, 0, 0, 0, -1, 0, 0, 1, 0, 1, 1, -1,

*0, 1, -1, -1, 1, 0, 2, 1)

DATA ((KL(I), I = 45, 66) = 2, 0, -1, -1, 1, -1, 1, 0, 0, 1, 1, 2, 0, 0, 1,

*0, 1, 2, 0, 1, 0, 1)

DATA ((KL(I), I = 67, 88) = 1, 1, -1, -2, 3, 0, 1, -1, 2, 1, 3, 0, -1, 1,

*-2, -1, 2, 1, 1, -2, -1, 1)

DATA ((KL (I), I = 89, 106) = 2, 2, 1, 0, 3, 1, 0, -1, 0, 0, 0, 1, 0, 1, 1,

*2, 0, 0)

DATA ((KL1 (I), I = 1, 22) = 0, 0, 0, 0, 0, -1, -2, 0, 0, 1, 1 -1, 0, 0, 0,

*2, 1, 2, -1, 0, -1, 0)

DATA ((KL1 (I), I = 23, 44) = 1, 0, 1, 0, 1, 1, 0, 1,

*0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

DATA ((KL1 (I), I = 45, 66) = 0, 0, 0, 0, 0, 0, 1, 1, -1,

*0, 0, 0, 0, 0, 0, 0, -1, 0, 1, 0, 0, 1)

DATA ((KL1 (I), I = 67, 88) = 0, -1, -1, 0, 0, -1, 1, 0, 0, 0, 0, 0,

*0, 0, 0, 0, 0, 1, 0, 0, 0, -1)

DATA ((KL1 (I), I = 89, 106) = 0, 0, 0, 0, 0, 0, 1, -1, 0, 0, 1, 0,

*-1, 1, 0, 0, 0, 1)

DATA ((KF (I), I = 1, 22) = 0, 0, 2, -2, 2, 0, 2, -2, 2, 0, 2, 2, 2,

*0, 2, 0, 0, 2, 0, 0, 2, 0)

DATA ((KF (I), I = 23, 44) = 2, 0, 0, -2, -2, 0, 0, 2, 2, 0, 2, 2, 0, 2,

*0, 0, 0, 2, 2, 2, 0, 2)

DATA ((KF (I), I = 45, 66) = 2, 2, 2, 0, 0, 2, 0, 2, 2, 2, 0, 2, 0,

*2, 2, 0, 0, 2, 0, -2, 0, 0)

DATA ((KF (I), I = 67, 88) = 2, 2, 2, 0, 2, 2, 2, 2, 0, 0, 0,

*2, 0, 0, 2, 2, 0, 2, 2, 2, 4, 0)

DATA ((KF (I), I = 89, 106) = 2, 2, 0, 1, 2, 2, 2, 0, -2, 2, 0,

*-2, 2, 0, -2, 0, 2, 0)

DATA ((KD (I), I = 1, 22) = 0, 0, 0, 0, 0, -1, -2, 0, -2, 0,

*-2, -2, -2, -2, -2, 0, 0, -2, 0, 2, -2, -2)

DATA ((KD (I), I = 23, 44) = -2, -1,-2, 2, 2, 0, 1, -2, 0, 0, 0, 0,

*-2, 0, 2, 0, 0, 2, 0, 2, 0, -2)

DATA ((KD (I), I = 45, 66) = 0, 0, 0, 2, -2, 2, -2, 0, 0, 2, 2, -2, 2, 2, -2,

*-2, 0, 0, -2, 0, 1, 0)

DATA ((KD (I), I = 67, 88) = 0, 0, 2, 0, 0, 2, 0, -2, 0, 0, 0, 1, 0,

*-4, 2, 4, -4, -2, 2, 4, 0, -2)

DATA ((KD (I), I = 89, 106) = -2, 2, 2, -2, -2, -2, 0, 2, 0, -1, 2,

*-2, 0, -2, 2, 2, 4, 1)

DATA (( (I), I = 1, 22) = 1, 2, 1, 0, 2, 0, 1, 1, 2, 0, 2, 2, 1,

*0, 0, 0, 1, 2, 1, 1, 1, 1)

DATA (( (I), I = 23, 44) = 1, 0, 0, 1, 0, 2, 1, 0, 2, 0, 1, 2, 0, 2, 0,

*1, 1, 2, 1, 2, 0, 2)

DATA ((KOM(I), I = 45, 66) = 2, 0, 1, 1, 1, 1, 0, 2, 2, 2, 0, 2, 1, 1, 1, 1,

*0, 1, 0, 0, 0, 0)

DATA (( (I), I = 67, 88) = 0, 2, 2, 1, 2, 2, 2, 1, 1, 2, 0, 2, 2, 0, 2, 2,

*0, 2, 1, 2, 2, 0)

DATA (( (I), I = 89, 106) = 1, 2, 1, 2, 2, 0, 1, 1, 1, 2, 0, 0, 1, 1, 0,

*0, 2, 0)

DATA ((APCI (I), I = 1, 22) = -171996., 2062., 46., 11., -3.,

*-2., 1., -13187., 1426., -517., 217., 129., 48., -22., 17., -15.,

*-16., -12., -6., -5., 4.)

DATA ((APCI (I), I = 23, 44) = 4., -4., 1., 1., -1., 1., 1., -1.,

*-2274., 712., -386., -301., -158., 123., 63., 63., -58., -59.,

*-51., -38., 29., 29.)

DATA ((APCI (I), I = 45, 66) = -31., 26., 21., 16., -13., -10., -7.,

*7., -7., -8., 6., 6., -6., -7., 6., -5., 5., -5., -4., 4., -4., -3.)

DATA ((APCI (I), I = 67, 88) = 3., -3., -3., -2., -3., -3., 2., -2.,

*2., -2., 2., 2., 1., -1., 1., -2., -1., 1., -1., -1., 1., 1.)

DATA ((APCI (I), I = 89, 106) = 1., -1., -1., 1., 1., -1., 1., 1., -1.,

*-1., -1., -1., -1., -1., -1., 1., -1., 1.)

DATA ((BPCI (I), I = 1, 22) = -174.2, 0.2, 0.0, 0.0, 0.0, 0.0, 0.0,

*0.0, -1.6, -3.4, 1.2, -0.5, 0.1, 0.0, 0.0, -0.1, 0.0, 0.1, 0.0,

*0.0, 0.0, 0.0)

DATA ((BPCI (I), I = 23, 44) = 0., 0., 0., 0., 0., 0., 0., 0., -0.2,

*0.1, -0.4, 0.0, 0.0, 0.0, 0.0, 0.1, -0.1, 0.9, 0.0, 0.0, 0.0, 0.0)

DATA ((BPCI (I), I = 45, 66) = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

DACA ((BPCI (I), I = 67, 88) = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

DATA ((BPCI (I), I = 89, 106) = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

DATA ((AEPC(I), I = 1, 22) = 92025., -895., -24., 0., 1., 0., 1., 0.,

*5736., 54., 224., -95., -70., 1., 0., 0., 9., 7., 6., 3., 3., -2.)

DATA ((AEPC (I), I = 23, 44) = -2., 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*977., -7., 200., 129., -1., -53., -2., -33., 32., 26., 27., 16.,

*-1., -12.)

DATA ((AEPC (I), I = 45, 66) = 13., -1., -10., -8., 7., 5., 0., -3.,

*3., 3., 0., -3., 3., 3., -3., 3., 0., 3., 0.0., 0.0., 0.0, 0.0)

DATA ((AEPC (I), I = 67, 88) = 0., 1., 1., 1., 1., 1., -1., 1., -1., 1.,

*0., -1., -1., 0., -1., 1., 0., -1., 1., 1., 0., 0.)

DATA ((AEPC (I), I = 89, 106) = -1., 0., 0., 0., 0.0, 0.0, 0.0, 0.0, 0.0,

*0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

DATA ((BEPC (I), I = 1, 22) = 8.9, 0.5, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*-3.1, -0.1, -0.6, 0.3, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*0.0)

DATA ((BEPC (I), I = 23, 44) = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*-0.5, 0.0, 0.0, -0.1, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

*0.0)

DATA ((BEPC (I), I = 45, 66) = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

DATA ((BEPC (I), I = 67, 88) = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,

0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

DATA ((BEPC (I), I = 89, 106) = 0.0, 0.0, 0.0, 0.0,

*0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)

TAY = D/36525.

2 = *

3 = 2*

EPCO = 0.4090928042 - 0.2269655*1. - 3* - 0.29*1. - 8*2 +

*0.88*1.E - 8*3

FL1 = 6.24003594 + 928.30195602* - 2.7974*1.E - 6*2 - 5.82*

*1.E - 8*3

F = 1.62790193 + 8433.46615831* - 6.42717*1. - 5*2 + 5.33*

*1. - 8*3

D1 = 5.19846951 + 7771.37714617* - 3.34085*.1. - 5*2 + 9,21*

*1. - 8*3

= 2.182438624 - 33.757045936* + 3.61429*1. - 5*2 + 3.88*

*1. - 8*3

IF (NP1.EQ.1) GO 3

FL = 2,355548393 + 8328.69142288* + 1.517952*1. - 4*2 +

*3.103*1.E - 7*3

HPCI = 0

HEP = 0

DO 2 I = 1,106

APG = KL (I)*FL + KL1 (I)*FL1 + KF (I)*F + KD (I)*D1 + KM (I)*OM

HPCI = HPCI + (APCI (I) + BPCI (I)*TA)*SIN (APG)

HEPC = HEPC + (AEPC (I) + BEPC (I)*TA)*COS (APG)

2 CONTINUE

HPCI = HPCI*4.848136811*1.E - 10

HEPC = HEPC*4.848136811*1.E - 10

GO TO 4

3 SOM = SIN (OM)

COM = COS (CM)

HPCI = -0,83386*1.E - 4*SOM + 1,9994*1.E - 6*COM*SOM - 0,63932*

*1.E - 5*SIN(2*(F - D1 + OM)) +0,6913*1.E - 6*SIN (FL1) - 0,11024*

*1.E - 5*SIN (2*(F + OM))

HEPC = 0,44615*1.E - 4*COM + 0.27809*1.E - 5*COS (2*(F - D1 + OM))

*+ 0.474*1.E - 6*COS (2*(F + OM))

4 CONTINUE

EPC = EPCO + HEPC

RETURN

END

********************************************************

SUBROUTINE (D, NP1, NP2, HM)

DIMENSION HM (3, 3)

CALL (D, NP1, HP, HE, E, E0)

SE = SIN (E)

C = COS (E)

IF (NP2.EQ.1) GO TO 1

SE0 = SIN (0)

CE0 = OS (E0)

SHPCI = SIN (HP)

CH = COS (HP)

HM (1, 1) = CHPCI

HM (1, 2) = -SHPCI*CE0

HM (1, 3) = -SHPCI*SE0

HM (2, 1) = SHPCI*CE

HM (2, 2) = CHPCI*CE*CE0 + SE*SE0

HM (2, 3) = CHPCI*CE*SE0 - SE*CE0

HM (3, 1) = SHPCI*SE

HM (3, 2) = CHPCI*SE*CE0 - CE*SE0

HM (3, 3) = CHPCI*SE*SE0 + CE*CE0

GO TO 2

1 HM (1, 1) = 1

HM (1, 2) = -HP*CE

HM (1, 3) = -HP*SE

HM (2, 1) = HP*CE

HM (2, 2) = 1

HM (2, 3) = -HE

HM (3, 1) = HP*SE

HM (3, 2) = HE

HM (3, 3) = 1

2 CONTINUE

RETURN

END

SUBROUTINE PREMAT (D, PM)

DIMENSION PM (3, 3)

TA = D/36525.

2 = *

3 = 2*

DZETA = 0,111808609*1.E - 1* + 0,146356*1. - 5*2 + 0,872*

*1. - 7*3

Z = 0,0111808609* + 0,53072*1.E - 5*2 + 0,883*1. - 7*3

= 0,97171735*1. - 2* - 0,20685*1. - 5*2 - 0.2028*

*1. - 6*3

SDZ = SIN (DZETA)

CDZ = COS (DZETA)

SZ = SIN (Z)

CZ = COS (Z)

ST = SIN (TETA)

CT = COS (TETA)

PM (1, 1) = CDZ*CZ*CT - SDZ*SZ

PM (1, 2) = -SDZ*CZ*CT - CDZ*SZ

PM (1, 3) = -CZ*ST

PM (2, 1) = CDZ*SZ*CT + SDZ*CZ

PM (2, 2) = -SDZ*SZ*CT + CDZ*CZ

PM (2, 3) = -SZ*ST

PM (3, 1) = CDZ*ST

PM (3, 2) = -SDZ*ST

PM (3, 3) = CT

RETURN

END

*********************************************************

SUBROUTINE RMUMAT (D, NP1, RM)

DIMENSION RM (3, 3)

TA = D/36525.

U = 0,02236172* + 0,67707*1. - 5** + 0,1755*1. - 6*

***

CALL (D, NP1, HP, HE, E, EO)

HALFA = HP*COS (E)

SMUHA = SIN (MU + HALFA)

CMUHA = COS (MU + HALFA)

RM (1, 1) = CMUHA

RM (1, 2) = SMUHA

RM (1, 3) = 0

RM (2, 1) = -SMUHA

RM (2, 2) = CMUHA

RM (2, 3) = 0.

RM (3, 1) = 0.

RM (3, 2) = 0.

RM (3, 3) = 1.

RETURN

END

4

,

1.

NDMG = 6051988 T(UT1) = 0.00 NP1 = 0 NP2 = 0

JD = 2447287 DJD = 5.0000 - 01 D = - 4.257500000 + 03 DM = 0.0000 + 0

SI = 3.910713393 + 00 SC = 3.910706227 + 00 SM = 3.913312726 + 00

9.9999999997 - 01 -7.1659192628 - 06 -3.1070314463 - 06

7.1657857874 - 06 9.9999999905 - 01 -4.2957047299 - 05

3.1073392701 - 06 4.2957022742 - 05 9.9999999907 - 01

9.9999596162 - 01 2.6064761738 - 03 1.1326991650 - 03

-2.6064761739 - 03 9.9999660313 - 01 -1.4761500974 - 06

-1.1326991650 - 03 -1.4762092497 - 06 9.9999935849 - 01

RMU

9.9999999997 - 01 7.1657857874 - 06 0.0000000000 + 00

-7.1657857874 - 06 9.9999999997 - 01 0.0000000000 + 00

0.0000000000 + 00 0.0000000000 + 00 1.0000000000 + 00

= *

9.9999598378 - 01 2.5993102834 - 03 1.1295921461 - 03

-2.5992617571 - 03 9.9999662091 - 01 -4.4426053158 - 05

-1.1297038036 - 03 4.1488766780 - 05 9.9999936102 - 01

2.

NDMG = 6051988 T (UT1) = 153045.15 NP1 = 0 NP2 = 0

JD = 2447288 DJD = 1.463569028 - 01 D = -4.2568536445 + 03

DM = 6.463559029 - 0

SI = 1.699821523 + 00 SC = 1.699813936 + 00 SM = 1.702420040 + 00

9.9999999996 - 01 -7.5873659227 - 06 -3.2897641450 - 06

7.5872245852 - 06 9.9999999905 - 01 -4.2960700284 - 05

3.2900901005 - 06 4.2960672545 - 05 9.9999999907 - 01

9.9999596285 - 01 2.6060804853 - 03 1.1325271999 - 03

-2.6060804853 - 03 9,9999660416 - 01 -1.4757019335 - 06

-1.1325271998 - 03 -1,4757610589 - 06 9.9999935869 - 01

RMU

9,9999999997 - 01 7.5872245852 - 06 0.0000000000 + 00

-7.5872245852 - 06 9.9999999997 - 01 0.0000000000 + 00

0.0000000000 + 00 0.0000000000 + 00 1.0000000000 + 00

= *

9.9999598631 - 01 2.5984931499 - 03 1.1292374490 - 03

-2.5984446347 - 03 9.9999662304 - 01 -4.4427781927 - 05

-1.1293490809 - 03 4.1493339839 - 05 9.9999936142 - 01

3.

NDMG = 23061987 T(UT1) = 0.00 NP1 = 0 NP2 = 0

JD = 2446969 DJD = 5.000000000 - 01 D = -4.575500000 + 03

DM = 0.000000000 + 0

SI = 4.723393770 + 00 SC = 4.723403754 + 00 SM = 4.726204930 + 00

9.9999999994 - 01 9.9845949170 - 06 4.3291889303 - 06

-9.9844175856 - 06 9.9999999911 - 01 -4.0959898342 - 05

-4.3295978945 - 06 4.0959853322 - 05 9,9999999915 - 01

9.9999533584 - 01 2.8011503947 - 03 1.2173044165 - 03

-2,8011503948 - 03 9.9999607676 - 01 -1.7048936386 - 06

-1.2173044164 - 03 -1.7049670601 - 06 9.9999925908 - 01

RMU

9.9999999995 - 01 -9.9844175857 - 06 0.0000000000 + 00

9.9844175857 - 06 9.9999999995 - 01 0.0000000000 + 00

0.0000000000 + 00 0.0000000000 + 00 1.0000000000 + 00

= *

9.9999530254 - 01 2.8111349429 - 03 1.2216335851 - 03

-2.8110849026 - 03 9.9999604797 - 01 -4,2676915707 - 05

-1.2217487277 - 03 3.9242597714 - 05 9.9999925289 - 01

4.

NDMG = 23061987 T (UT1) = 121212.00 NP1 = 0 NP2 = 0

JD = 2446970 DJD = 8.472222222 - 03 D = -4.574991531 + 03

DM = 5.084722222 - 01

SI = 1.643780943 + 00 SC = 1.643790723 + 00 SM = 1.646591587 + 00

9.9999999994 - 01 9.7803786014 - 06 4.2406433872 - 06

-9.7802053348 - 06 9.9999999911 - 01 -4.0856676151 - 05

-4.2410429773 - 06 4.0856632950 - 05 9.9999999915 - 01

9.9999533688 - 01 2.8008391198 - 03 1.2171691360 - 03

-2.8008391198 - 03 9.9999607764 - 01 -1.7045147415 - 06

-1.2171691359 - 03 -1.7045881385 - 06 9.9999925925 - 01

RMU

9.9999999995 - 01 -9.7802053348 - 06 0.0000000000 + 00

9.7802053348 - 06 9.9999999995 - 01 0.0000000000 + 00

0.0000000000 + 00 0.0000000000 + 00 1.0000000000 + 00

= *

9.9999530427 - 01 2.8106194526 - 03 1.2214097595 - 03

-2.8105695476 - 03 9.9999604942 - 01 -4.2573064790 - 05

-1.2215245909 - 03 3.9140006079 - 05 9.9999925317 - 01

 

1.

.. . .. , .. , .. , .. , .. , .. , .. , .. , .. , .. , .. , ..

2. 30.05.89 1375

3. - 1995 .

4.

5. - :

,

, , ,

15855-77

1 ( 1)

 

 

 

 






, .
, , .
150000 - , .
- .

@Mail.ru

Copyright www.gostrf.com, 2008 - 2017