С. 24 ГОСТ 25645.115—84
ТЕКСТ ПРОГРАММЫ
SUBROUTINE DBNSIT (ХК, ТН, SZ, НКМ, SUN,
* АР, FI, FSH, Г, PR, RO, RU)
IMPLICIT REAL*8(A—Н, O -Z )
REAL*® Kl, К2, КЭ, К4, КО, КР
DIMENSION ХК(3), SUN (2)
COMMON/UKOEF/U(8, 4), GR, WZ
COMMON/QKOEF/Q (27,3)/YEAR/Р (38)
H=HKM
RAD= DSQRT(XK(1)*XK(1)+XK(2)*XK(2)+XK(3)*XK(3))
X=XK(1)/RAD
Y= XK(2)/RAD
Z= XK(3)/RAD
IF (HLT.120.D0) GOTO 5
J=1
IF(H.GE.I8OD0) J = 2’
IF (H GE 600 DO) J= 3
HD= H*H
RO= DEXP(Q(2, J)—Q(3, J)*DSQRT(H—Q(4, J)))
RO=RO*GR
K1 = Q(8, J)+Q(9, J)*H+Q(10, J)*HD+Q(1T,
J)*HD*H
SD= DSIN(SUN(2))
CD= DCOS(SUN(2))
G= SUN(1)+Q(14, J)—SZ—WZ*(TH—10800.00)
SG= DSIN(G)
CG=DCOS(G)
CF=Z*SD+[CD*(X*CG+Y*SG)
CP=((1D0+CF)*&D0)**((Q(1>2, J)+Q (13, J)*H)*.5DO)
Kl’= 1.D0+K1*CP
RO= Kl*RO
I= INT (D/ lOiDO)+1)
SI = P(1) + (P(I+1)—P(I))*(D—(I—1) *10 DO)/10.DO
K2=Q(16, J)+Q(16, J)*H+Q(17, J)*HD
K2=1.D0+K2*SI
RO= K2*RO
K3=Q(5, J)+Q(6, J)*H+Q(7, J)*HD
A = (FI—FSH)/FSH
K3=l.DO+K3*A
RO= K3*RO
A=FSH—Q(l, J)
KO=Q(25, J)+Q(26, J)*H+Q(27, J)+HD
KO=l.DO+KO*A
RO= KO*RO
IF (PR GT 0.D0) GOTO 3
CALL INTERP(AP, KP)