q
С. б ГОСТ 25645.150-90
ПРОГРАММА РАСЧЕТА ПОТОКА ЧАСТИЦ ГКЛ
с
PROGRAM SP29
j
DIMENSION EN(4I). КГ41|. FR(4I). F<4I)
COMMOX/PART/A(29). D(29>. Af.P{29). GAM(291. TO. AM
COMMON ’SUNSPAV <4551
СВЫЧИСЛЕНИЕ СПЕКТРОВ ЯДЕР ЖЕЛЕЗА НА Г, ЯНВАРЯ 1967 ГОДА
DATA IZ. T/2S. 1987 04,’
СЗАДАЕТСЯМАССИВЗНАЧЕНИИЭНЕРГИИ.РАВНОМЕРНЫЙ НА ЛОГАРИФМИЧЕСКОИ
СШКАЛЕ
ГЮ 1 1-1.41
U-FLOAT(I)
1 ENG)-1 0 **’<ГЪ* 1.)М )»10
СВЫЧИСЛЕНИЕ ПОТОКОВ ЧАСТИЦ
DO 2 1-1.41
2 CALI. FLXGCRfIZ. EN(I). T. R(I). FR(I). F(l>)
СРАСПЕЧАТКА РЕЗУЛЬТАТОВ ВЫЧИСЛИНИИ
W R lT EG .tniZ . A(IZ). AM. T
И FORMAT!IBX. ’ЛУЧИ КОСМИЧЕСКИЕ ГАЛАКТИЧЕСКИЕ’
*#X. 10НЧАСТИЦЫ Z -. 12. ЗН A -. F4 1. ЮН МАССЫ AM *. ЕЮ.З.
•1III ГЭВДНУКЛ). 4Н В . F7.2. 5И ГОЛУ)
WRITE П. 12)
12 FORMAT!IЗХ. 50Н ЭНЕРГИЯЖЕСТКОСТЬСПЕКТР Ф(й)СПЕКТР I <Е>./
•I3X. 58НГМЭВ/ПУКЛ1ГГВ1(С*М**2*СР’ГВ] [С*.И’*2»СР*МЭВ/НУКЛ0
WRITE (I. 131(1. FN
(I).
Rfl). FR!I>, F(I). 1-1.41)
13 FORMATS!(5X. 13. 4ЦРЕ13.1)/))
END
SUBROUTINE FLXGCRflZ. EN. T. R. FR. F)
COMMON ‘SUNSPAV (456)
COMMON/PART/A(29). D(29l. ALF(29). GAM(29). TO. AM
IFfIZ. LE. 28)AM- 0.939
IFfIZ. EQ.29)AM~5. IE—04
JFfIZ.LE. 28)0=1.
IFfIZ. EQ.29)0——I.
IFfIZ. !.F. 28)T1 = 1982. 5
IF(IZ. EQ. 29)T I- 1980.5
IFfIZ I F 28)Z~ FLOAT(IZ)
IFfIZ. EQ.29)Z—I.
СОПРЕДЕЛЕНИЕ ПОРЯДКОВОГО НОМЕРА МЕСЯЦА NT В ТАБЛИЦЕ <W.
DO 1 I—1. 480
U—FLOATfl)
AT—1949.96-r 0.083* LI
IF(T L0 AT) GOTO 2
1 CONTINUE
2 NT-1
СВЫЧИСЛЕНИЕ ЖЕСТКОСТИ И СКОРОСТИ ЧАСТИЦЫ
R—А(IZ)/7*SQRT( (FN/1000. Ь2. *АМ)*Е\/101» >
В—R/SQRT !Р**2+ (AM’AflZ) fZ)**2)
С ВЫЧИСЛЕНИЕ ВЕЛИЧИНЫ ЗАПАЗДЫВАНИЯ ПОТОКОВ - NM [.МЕСЯЦЕВ]
X=SINf. 285* (Т—ТО))
s t
—(
a b s
(X) 1**. ззз
ЯМ—SICN(ST. X)
D T-I0. +5. -SN
NM -lN TfD T+ 5)
С ОПРЕДЕЛЕНИЕ ЭФФЕКТИВНОГО ЧИСЛА СОЛНЕЧНЫХ ПЯГЕН
AN= (FLOATfNM) -И . )#2.
LI —NT—NM*2
1.2* NT
S—о
ss-o.
v-o.
DO
4
.1*1.1,12
V -V + l.
U=V
IFfV.GT. NM+ 1) U -A N —V
SS-SS+U
4S -S + U*W(J)
WQ-S/SS
R0~. 375+ 3.0E—4)4•WQ**1. 455
X-S1N E 285*(T—TI))
ST* (ABSfX))**. 333
S N - SIGN (ST. X)
Y -B’R /R0
IF(Y. GT. 15. )Y -I5.