С**************************************************
С PROGRAM SP29
С**************************************************
DIMFNSION EN(41), R(41), FR(41), F(41)
COMMON/PART/A(29), D(29), ALF(29), GAM(29), T0, AM
COMMON/SUNSPW(456)
С ВЫЧИСЛЕНИЕ СПЕКТРОВ ЯДЕР ЖЕЛЕЗА НА 15 ЯНВАРЯ 1987 ГОДА
DATA IZ, Т/26, 1987 04/
С ЗАДАЕТСЯ МАССИВ ЗНАЧЕНИИ ЭНЕРГИИ, РАВНОМЕРНЫЙ НА ЛОГАРИФМИЧЕСКОЙ
С ШКАЛЕ
DO 1 I=1,41
U=FLOAT(I)
1 EN (I)=10 **((U-1.)*.* 10
C ВЫЧИСЛЕНИЕ ПОТОКОВ ЧАСТИЦ
DO 2 I=1.41
2 CALL FLXGCR (IZ, EN(I), T, FR(I), F(I))
C РАСПЕЧАТКА РЕЗУЛЬТАТОВ ВЫЧИСЛЕНИЙ
WRITE (1,11) IZ A(IZ),AM,T
11 FORMAT (18X, ЛУЧИ' КОСМЕТИЧЕСКИЕ ГАЛАКТИЧЕСКИЕ'
*/5X, 10НЧАСТИЦЫ Z=, 12,3Н А=. F4 1, 10Н МАССЫ АМ =,Е10.3,
*11 Н ГЭВ/(НУКЛ), 4Н В ,F7.2, 5H ГОДУ)
WRITE (1. 12)
12 FORMAT (13X, 50Н ЭНЕРГИЯ ЖЕСТКОСТЬ СПЕКТР Ф(R) СПЕКТР F(E),/
*13Х. 58Н[МЭВ/НУКЛ] [ГВ] [С*М**2*СР*ГВ] [С*М**2*СР*МЭВ/НУКЛ])
WRITE (1. 13) I, EN(I), R(I), FR(I), F(I), 1 = 1,41)
13 FORMAT(5I (5X, 13, 4(1PE13.1)/))
END
SUBROUTINE FLXGCR(IZ, EN, T, R, FR, F)
COMMON/SUNSP/W (456)
COMMON/PART/A(29), D(29), ALF(29), GAM(29), TO, AM
IF(IZ. LE. 28) AM=0.939
IF(IZ. EO 29) AM=5. IE-04
IF(IZ. LE. 28)Q=1.
IF(IZ. EQ.29)Q = - 1.
IF(IZ. LF 28)T1 = 1982.5
IF(IZ. EQ.29)T1 = 1980.5
IF(IZ IE 28)Z=FLOAT(IZ)
IF(IZ. EQ.29)Z=1.
С ОПРЕДЕЛЕНИЕ ПОРЯДКОВОГО НОМЕРА МЕСЯЦА NT В ТАБЛИЦЕ "W"
DO 1 I = l, 480
U=FLOAT(I)
AT= 1949.96+0,083*U
IF (Т. LE AT) GOTO 2
1. CONTINUE
2 NT=I
С ВЫЧИСЛЕНИЕ ЖЕСТКОСТИ И СКОРОСТИ ЧАСТИЦЫ
R = A(IZ)/Z*SQRT((EN/1000. +2 *AM)*EN/1000 )
B = R/SQRT(R**2 + (АМ*А (IZ) *Z) **2)
С ВЫЧИСЛЕНИЕ ВЕЛИЧИНЫ ЗАПАЗДЫВАНИЯ ПОТОКОВ - NM [МЕСЯЦЕВ]
X=SIN (.285* (T-T0))
ST=(ABS (X))**.333
SN=SIGN (ST.X)
DT=10.+5. *SN
NM=INT (DT+5)
C ОПРЕДЕЛЕНИЕ ЭФФЕКТИВНОГО ЧИСЛА СОЛНЕЧНЫХ ПЯТЕН
AN=(FLOAT(NM)+1)*2.
L1=NT-NM*2
L2=NT
S=0
SS=0.
V=0.
DO 4 J=L1, L2
V=V+l.
U=V
IF(V.GT.NM+l) U=AN-V
SS = SS + U
4 S = S + U*W(J)
WQ=S/SS
R0=. 375+3.0E-04*WQ**1. 455
X = SIN(.285*(T-T1))
ST=(ABS(XH))** 333
SN=SIGN(ST, X)
Y=B*R/R0
IF (Y.GT. 15.)Y=15.
IF(IZ.EQ.29.)WW = B*R/. 106
IF(IZ.LE.28.)WW = B*R/.012*Z/A(IZ)
IF (WW.GT.15.)WW=15.
IF(1Z.EQ.29.)VR=5.5*(1.-.9/EXP(WW))
IF(IZ.LE.28.)VR = 5.5*(l.-1.2/EXP(WW))
GAPA = VR+1.13*Q*Y/EXP(Y)*SN
YY=R
IF(YY.GT.15.)YY=15.
IF(IZ. EQ.29)GAM(IZ) = 3.0-1.4/EXP(YY)
F0=D(IZ)*B**ALF(IZ)/R**GAM(IZ)
FR=F0*(R/(R + R0))**GAPA
F=A(IZ)/Z*FR/B/1000.
RETURN
END
BLOCK DATA
DIMENSION W50(120), W60(120), W70(120), W80(96)
COMMON/SUNSP/W (456)
CCMMON/PART/A(29), D(29), ALF(29), GAM(29), T0, AM
EQUIVALENCE(W(1), W50(l)), (W(121), W60(l)), (W(241), W70(l))
EQUIVALENCE(W(361), W80(l))
DATA W50/
*101.6, 94,8, 109.7, 113.4, 106.2, 83.6, 91.0, 85.2, 51.3, 61.4, 54.8, 54.1,
*59.9, 59.9, 55.9, 92.9, 108.5, 100.6, 61.5, 61.0,83.1,51.6,52.4,45.8,
*40.7, 22.7, 22.0,29.1, 23.4, 36.4,39.3,54.9,28.2,23.8,22.1,34.3,
26.5, 3.9, 10.0, 27.8, 12.5, 21.8, 8.6, 23.5, 19.3, 8.2, 1.6, 2.5,
* 0.2, 0.5,10.9, 1.8, 0.8, 0.2, 4.8, 8.4, 1.5, 7.0, 9.2, 7.6,
* 23.1, 20.8, 4.9, 11.3, 28.9, 31.7, 26.7, 40.7, 42.7, 58.4, 89.2,76.9.
* 73.9, 124.0, 118.4, 110.7, 136.6, 116.6, 129.1, 169.6, 173.2, 155.3, 201.3,
* 192.1
* 165.0, 170.2, 157.4, 175.2, 164.6, 200.7, 187.2, 158.0,235.8,253.8,210 9,
* 239. 4,
* 202.5, 164.9, 190.7, 196.0, 175.3, 171.5, 191.4, 200.2, 201.2, 181.5, 152.3,
* 187.6,
* 217.4, 143.1, 185.7, 163.3, 172.0, 168.7, 149.6, 199.6, 145.2, 111.4, 124.0,
* 125.0/
DATA W60/146.3,
* 106.0, 102,2, 122 0, 119.6, 110.2, 121.7, 134.1, 127.2,82.8,89.6,85.6,
* 57.9, 46.1, 53.0, 61.4, 51.0, 77.4, 70.2, 55.9, 63.6, 37.7, 32.6, 39 9,
* 38.7, 50.3, 45.6, 46.4, 43.7, 42.0, 21.8, 21.8, 51.3, 39.5, 26.9, 23.2,
*19.8, 24.4, 17.1, 29.3, 43.0, 35.9, 19.6, 33.2, 38.8, 35.3, 23.4, 14.9,
* 15.3, 17.7, 16.5, 8.6, 9.5, 9.1, 3.1, 9.3, 4.7, 6.1, 7.4, 15.1,
* 17.5, 14.2, 11.7, 6.8, 24.1, 15.9, 11.9, 8.9, 16.8, 20.1, 15.8, 17.0,
* 28.2, 24.4, 25.3, 48.7, 45.3, 47.7, 56.7, 51.2, 50.2, 57.2, 57.2, 70.4,
* 110.9, 93.6, 111.8, 69.5, 86.5, 67.3, 91.5, 107.2, 76.8, 88.2, 94.3, 126.4,
* 121.8, 111.9, 92.2, 81.2, 127.2, 110.3, 96.1, 109.3, 117.2, 107.7,86.0,
* 109.8,
* 104.4, 120.5, 135.8, 106.8, 120. , 106. , 96.8, 98. , 91.3, 95.7, 93.5, 97.9
DATA W70111.5,
* 127.8, 102.9, 109.5, 127.5, 106.8, 112.5, 93.0, 99.5, 86.6, 95.2, 83.5,
* 91.3, 79.0, 60.7, 71.8, 57.5, 49.8, 81.0, 61.4, 50.2, 51.7, 63.2, 82.2.
* 61.5, 88.4, 80.1, 63.2, 80.5, 88.0, 76.5, 76.8, 64.0, 61.3, 41.6,45.3,
* 43.4, 42.9, 46.0, 57.7, 42.4, 39.5, 23.1, 25.6, 59.3, 30.7, 23.9, 23.3,
*27.6, 26.0, 21.3, 40.3, 39.5, 36.0, 55.8, 33.6, 40.2, 47.1, 25.0,20.5,
* 18.9, 11.5, 11.5, 5.1, 9.0, 11.4, 28.2, 39.7, 13.9, 9.1, 19.4, 477.8,
*8.1, 4.3, 21.9, 18.8, 12.4, 12.2, 1.9, 16.4, 13.5, 20.6, 5.2, 15.3,
* 16.4, 23.1, 8.7, 12.9, 18.6, 38.5, 21.4, 30.1, 44.0, 43.8, 29.1, 43.2,
* 51.9, 93.6, 76.5, 99.7, 82.7, 95.1, 70.4, 58.1, 138.2, 125.1, 97.9, 122.7,
* 166.6, 137.5, 138.0, 101.5, 134.4, 149.5, 159.4, 142.4, 188.4, 186.2, 183.3,
* 176.3/
DATA W80/159.6, 155.0, 126.2,
* 164.1, 179.9, 157.3, 136.3, 135.4, 155.0, 164.7, 147.9, 174.4, 114.0,
* 141.3, 135.5, 156.4, 127.5, 90.9, 143.8, 158.7, 167.3, 162.4, 137.5, 150.1,
* 111.2,
* 163.6, 153.8, 122.0, 82.2, 110.4, 106.1, 107.6, 118.8, 94.7, 98.1, 127.0,
* 84.3, 51.0, 66.5, 80.7, 99.2, 91.1, 82.2, 71.8, 50.3, 55.8, 33.3, 33.4,
* 57.0, 85.4, 83.5, 69.7, 76.4, 46.1, 37.4, 25.5, 15.7, 12.6, 22.4, 18.2,
* 16.5, 15.9, 17.2, 18.1, 27.4, 24.2, 30.7, 11.1, 3.9, 18.5, 16.6, 17.3,
* 2.5, 23.5, 15.1, 18.5, 13.7, 1.1, 18.1,7.4, 3.8, 35.4, 15.2, 6.8,
* 10.4, 2.4, 14.7, 39.6, 33.0, 17.4, 33.0, 38.6, 33.5, 61.1, 40.9, 55.0/
DATA TO 1978.5
DATA A/1. , 4. , 6.9, 9 . 10.8, 12. , 14. , 16. , 19. ,20.2,
*23. , 24.3, 27. , 28.1, 31. , 32.1, 35.4, 39.9, 39.1, 40.1,44.9, 47.9,
* 50.9, 52.9,54.9, 55.8, 58.9, 58.7.1./
DATA D/2.0EO4, 3. 5E03, 1.7E + 01, l.6+01, 5.1E+01, 9.6E + 01, 3 5E + 01,
* 8.4E+01, 36Е00,1.5Е+01,4.2Е00, 1.8Е+01, 39Е00, 1.2Е+01, 1.0Е00,
* 2.7Е00, 1.2Е00, 2.3Е00, 1.8Е00, 2.6Е00, 6.9Е-01, 2.5Е00, 1.13Е00,
* 2.1Е00,1.40Е00, 9.2Е00, 8.7Е-02, 4.5Е-01,1.7Е02/
DATA ALF/3. , 3,. , 3.4, 4.5, 3.9, 3.1, 3.6, 3.0, 3.8, 3.1, 3.4, 3.0, 3.2,
* 3.0, 4.0, 3.4, 4.5, 4.5, 4.2, 3.2, 3.6, 3.6, 3.3, 3.3, 3.0, 3.1, 4.0, 3.2, 3./
DATA GAM/2.75, 2.75, 2.7, 2.9, 3.0, 2.75, 2.9, 2.7, 3.0, 2.75, 2.9, 2.7,
* 2.8, 2.65, 2.95,2.7, 3.0, 2.9, 3.0,2.75, 2.9, 2.95, 2.9, 2.85, 2.70, 2.6,
* 2.75, 2.6,3./
END
1. Рассчитывают жесткостной и энергетический спектры ядер железа ГКЛ на 15 января 1987 г.
2. Жесткостной и энергетический спектры частиц ГКЛ вычисляют при помощи программы SP29, в которой вызывается расчетная подпрограмма FLXGCR. Чтобы вычислить жесткостной и энергетический спектры ядер железа на 15 января 1987 г., нужно в программе SP29 в операторе DATA IZ, Т задать значения порядкового номера элемента железа (26) и время, соответствующее 15 января 1987 г. - 1987.04: DATA IZ, T/26, 1987.04/.
3. Результаты расчета жесткостного и энергетического спектров ядер железа ГКЛ на 15 января 1987 г. приведены в табл.2. Здесь EN(41), R(41), FR(41) и F(41) - массивы значений энергии, жесткости, жесткостного и энергетического спектров соответственно, размерность которых определяется требуемым количеством точек энергетического спектра (в примере - 41).
Комментарии (0)
Чтобы оставить комментарий вам необходимо авторизоваться