В указанной последовательности эти параметры помещаются в форм'ате 313 языка Фортран.
Вслед за первой записью следует KF групп записей, каждая из которых содержит данные о действии того или иного фунгицида, включая и контрольную - не содержащую фунгицид.
Первая запись в каждой группе содержит в формате 13 языка Фортран количество концентраций (7V) данного фунгицида, включая данные о контрольных пробах, не содержащих фунгицид (N < 20). Далее следует N подгрупп, каждая из которых содержит лянные о каждой концентрации
Первая запись в подгруппе - заголовок подгруппы, содержит сведения о концентрации фунгицида ‘ (С), лаг-фазе (LE), количестве отборов (Р) и количестве проб в каждом отборе (FN). Все перечисленные данные в указанной последовательности содержатся в записи в формате 416 языка Фортран (Р < 20) .
Вслед за заголовком подгруппы располагаются Р пар записей, содержащих данные о каждом из отборов. Первая запись в упомянутой паре записей содержит данные о времени (TR), прошедшем с начала эксперимента до взятия отбора и количества проб (FIR) в данном отборе1. Значения TR и FIR представляют в записи в формате 216 языка Фортран. Вторая запись содержит данные о биомассах не более чем в 16 пробах каждого из отборов. Значения биомасс в каждом отборе представляются в записи в формате I6F 5.2.
Выходные данные
Выходом программ является распечатка поступающих на обработку данных, а также сообщения выдаваемые программой в процессе счета. Все выдаваемые программой сообщения и промежуточные результаты счета не нуждаются в дополнительных комментариях.
Текст вызывающей (первой) подпрограммы
С ***********************************************************************
INTEGER Р(20), С(20), LE(20), Т(20, 20), N, KF, N0, Pft TR, FIR, FI, FN, F, Pl C C C С ПРОГРАММА КОНТРОЛЯ
С С ВВОДА И ФОРМИРОВАНИЯ
С С МАССИВОВДАННЫХДЛЯ
С С РАСЧЕТАПАРАМЕТРОВ
С С ЗАВИСИМОСТИ РОСТА
С С БИОМАССЫ ВО ВРЕМЕНИ
С
REAL М(20, 20), МВ(20), R(20) С ВВОД: KF - КОЛИЧЕСТВО ФУНГИЦИДОВ; С Р0 - ПРЕДЕЛЬНО ДОПУСТИМОЕ ЧИСЛО ОТБОРОВ; (МИНИМАЛЬНОЕ) С N0-ПРЕДЕЛЬНО ДОПУСТИМОЕ ЧИСЛО КОНЦЕНТРАЦИЙ; С (МИНИМАЛЬНОЕ)
С с ***********************************************************************
31 FORMAT (IX, 100(’*’)/’ * ДАННЫЕ ДЛЯ ’,13,’ - ОЙ КОНЦЕНТРАЦИИ’, *’ ОБРАБАТЫВАТЬСЯ НЕ БУДУТ *’,/1Х, 100(’*’))
С************************************************************************ с
READ (5, 1) KF, Р0, N0
1 FORMAT (313, F7, 5) PRINT 101 101 FORMAT (’KF P0 N0') PRINT 1, KF, P0.N0 C Q ***|И:П**М******************************
C * ЦИКЛ ПО ФУНГИЦИДАМ *
c *****************************************
c
DO 999 LL= 1, KF C
С ВВОД: N - ЧИСЛО КОНЦЕНТРАЦИЙ ДЛЯ ДАННОГО ФУНГИЦИДА (N<= 20).
С
READ (5, 2) N
2 FORMAT (13) IF (N.GE.N0) GO TO 3
9 PRINT 71, LL, N, N0 GO TO 999
710 FORMAT (’ КОЛИЧЕСТВО КОНЦЕНТРАЦИЙ N=’, 12)
71 FORMAT (’ **** ДЛЯ ФУНГИЦИДА ’, 12,’ РАСЧЕТ HE ВОЗМОЖЕН’ *, ’(N=’,12,’ МЕНЬШЕ N0 ’,12,’)’)
3 PRINT 710, N
С
С ЦИКЛ ПО ЧИСЛУ КОНЦЕНТРАЦИЙ
С ДЛЯ ДАННОГО ФУНГИЦИДА
С D0 99 I = l,N
с
С ВВОД: С - ЗНАЧЕНИЕ КОНЦЕНТРАЦИИ;
С LE - ЛАГ-ФАЗА ЭКСПЕРИМЕНТАЛЬНАЯ;
С Р - ЧИСЛО ОТБОРОВ ДАННОЙ КОНЦЕНТРАЦИИ;
С FN - ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ (НЕ БОЛЕЕ 16);
С
10 READ (5, 4) C(I), LE(I), P(I), FN
4 FORMAT (416) Р1=Р(1)
41 FORMAT (’*********** КОНЦЕНТРАЦИЯ’, 14, 10Х,
*’***********’)
PRINT 41,1
PRINT 42
PRINT 4, C(I), LE(I), P(I), FN
C C C
42 FORMAT (4X, ’C LE P FN’)
IF (P(I).GE.P® GO TO 5
KI =1
21 PRINT 72, C(I), P(I), P0
72 FORMAT (’ ******** ДЛЯ ДАННОЙ КОНЦЕНТРАЦИИ С =’, 16,
*’ ОБРАБОТКА НЕ ПРОИЗВОДИТСЯ’,
*/ 15Х,’ (Р =’,16,’ МЕНЬШЕ Р0 =’, 16,’) ’)
С
С ПРОГОН ДАННЫХ ЭТОЙ КОНЦЕНТРАЦИИ
DO 6 К = К1,Р1
READ (5, 7) TR, FIR
7 FORMAT (216)
PRINT 7, TR, FIR
READ (5, 8) (R(J), J = 1, 16)
8 FORMAT (16F5.2)
PRINT 8, (R(J), J = 1, 16)
6 CONTINUE
8
С УМЕНЬШИТЬ ЧИСЛО КОНЦЕНТРАЦИЙ ДАННОГО ФУНГИЦИДА
С
И N = N-1
IF (N.LT.N0) GO ТО 9
GO ТО 10
С 7
С ЦИКЛ ПО ЧИСЛУ ОТБОРОВ ДАННОЙ КОНЦЕНТРАЦИИ
С
5 DO 12 К=1,Р1<J<_> О и ио 19 13 14 15 18 23 20 16 22 17 24 25 12 |
READ (5,7)T(I,K), FI IF (FI.EQ.0) GO ТО 13 F = FI GO TO 14 IF (FN, EQ.(J) GO TO 15 F = FN IF (F.GT. 16) GO TO 16 GOTO 17 PRINT 18, T (I, K) FORMAT (’**N3** ДАННЫЕ ДЛЯ T =’, 16,’ HE ОБРАБАТЫВАЮТСЯ Т.К.' * ; ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ РАВНО НУЛЮ’) READ (5,8) (R(J), J = 1, 16) PRINT 8, (R(J), J = 1,16) P(I) = P(I) - 1 IF (P(I), GE.P0) GO TO 19 PRINT 20 FORMAT (IX, 7j0(’-’)) K1 = K+1 GO TO 21 PRINT 22, T (I, K) FORMAT (’ **N4** ДАННЫЕ ДЛЯ T=’,16,’ HE ОБРАБАТЫВАЮТСЯ T.K.’ *; ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ БОЛЬШЕ 16’) GO ТО 23 PRINT 7, Т (I, К), F READ (5,8) (R(J), J = 1, F) PRINT 24, (R(J), J = 1, F) FORMAT (’ **’, 16 (IX, F5.2)) M(I, K) = 0 DO 25 KK=1,F M(I, K) = M(I, K) + R(KK) CONTINUE M(I,K) = M(I, K)/F * CONTINUE |
28 27 26 99 999 30 |
ЗАПОЛНЕНИЕ МАССИВА MB К = 2 IF ((M(I, K) - M(I, K-1)).LT.0) GO TO 26 IF (K.EQ.P(I) GO TO 27 K = K+ 1 GO TO 28 MB(I) = M(I, K) P(I) = P(I)-1 GO TO 99 MB(I) = M(I, K-l) P(I) = K-2 • IF (P(I).LE.P(J) GO TO 30 CONTINUE CALL DATVLA (C, LE, MB, M, T, N, P, P0, N0) CONTINUE STOP PRINT 31,1 GO TO 99 END |
ВВОД: Т - ВРЕМЯ С НАЧАЛА РОСТА БИОМАССЫ FI - ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ
1
С
С С - КОНЦЕНТРАЦИЯ С LE - ЛАГ-ФАЗА ЭКСПЕРИМЕНТАЛЬНАЯ ***** **
С Т - ВРЕМЯ НАБЛЮДЕНИЯ *** *
С Р - ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ *** **
с *** ***
с *** ***
Q ******
Q ****
с ♦***
с ******
с *** ♦♦*
с ***
Q ♦♦♦ *♦
Q ***** **
с
SUBROUTINE DATVLA (С, LE, МВ, М, Т, N, Р, Р0, N0) INTEGER МКО(20), С(20), LE(20), Т (20, 20),Р(20), ММ(20), TT(2Qf) INTEGER N.PI.N0P0
REAL F(20), M(20, 20), MB(20), KL, КС, А(2® ,В(20), LR(20), SAM(20)
PRINT 1013
1012 FORMAT (2(/lX, 100(’*’)))
1013 FORMAT (3(/lX, 100(’*’)))
102 FORMAT (6X, ’Г, 5X, ’2’, 5X, ’З’, 5X, ’4’, 5X, ’5’
5X, ’б’, 5X, ’Г, 5X, ’8’, 5X, ’9’, 4X, ’10’, 4X
; ’П’, 4X, ’12’, 4X, ’ІЗ’, 4X, ’14’, 4X, ’15’, 4X,
’16’, 4X, ’17’, 4X, ’18’, 4X, ’19’, 4X, ’20’)
103 FORMAT (’С’, 14, 19 (IX, 15))
104 FORMAT (’LE’, 14,19 (IX, 15))
105 FORMAT (’MB’, F4.2, 19 (IX, F5.2))
Yc c *
С С ЦИКЛ ПО КОНЦЕНТРАЦИЯМ *
с с *
PRINT 103, (C(J), J= 1,N)
PRINT 104, (LE(J),J =1, N)
PRINT 105, (MB(J), J = 1, N)
C C *
C C ♦
1020 FORMAT (2(/lX, 100(’*’)), /50Х, ’PAC4ETKL’)
C
C ,
NKO = 0 NK = N DO 1001 = 1,NK PI = P(I) 1000 PRINT 1012 110 FORMAT (/IX, 20(’=’), ’ ДАННЫЕ ДЛЯ КОНЦЕНТРАЦИИ ’, 12, 20 (’=’)) 111 FORMAT (’М’, F4.2, 19(1X,F5.2)) 112 FORMAT (’Т’, 14, 19(1X, 15)) 113 FORMAT (’М-ЭКСП.МАССА, MR-PAC4.MACCA, Т-ВРЕМЯ, СКО-ОТКЛОНЕНИЕ’) с *
YС ♦
С ПЕЧАТЬ ИСХОДНЫХ ДАННЫХ, УЧАСТВУЮЩИХ В РАСЧЕТАХ А И В »
С *
С ' *
151 PRINT 110,1
PRINT 102
PRINT 112, (T(I, J), J = 1, P I)
PRINT 111, (M(I, J), J = 1,PI)
YC *
c *
С РАСЧЕТ КОЭФФИЦИЕНТОВ А И В *
С *
с : *
SE = 0 SC = 0 SR = 0 SD = 0 SL = LE® DO = 120 J= 1,PI ST = T(I, J)
F(J) - ALOG ((MB® - M(I, J))/M(I, J)) H = ST - SL SE = SE + F(J) SC = SC + H SD = SD + F(J) * H SR = SR + H*H
120 CONTINUE
PT = PL
B(I) = (SC»SE - PT«SD)/(PT*SR—SC* SC)
A® = (SE»SR - SC*SD)/(PT*SR - SC*SC)
106 FORMAT (’+’, 50X,’ LN(A)=’, F12.7,’ B=’,F12.7,’ CKO=’, E16.10) GO TO 3000
1 FN распространяются на всю подгруппу, FIR - на конкретный отбор. Если указаны оба, то предпочтение отдается FIR. Нельзя, чтобы FN и FIR отсутствовали. FN и FIR не должны превышать 16.