РСН 67-87 сторінка 5 з 11
424
IF (HPRI(1).GT.0.0) IPR=IPR+1
425
TPR=TPR+HPRI(IPR)
426
3
KTP=KT-1
427
IBRP=IBR-1
428
IF (KP(1).EO.0) GO TO 143
429
DO 143 LP=1,KW
430
K=KP(LP)
431
DO 145 I=1,NI
432
DO 145 J-1,NJ
433
N1=N(I,J,K)
434
HN1=H(I,J,K)
435
IF ((H1(N1).LE.HN1).AND.(HN1.LE.H2(H1))
436
1
UK(I,J,K)=700.0+100*(HN1-H1(N1))/(H2(N1)-H1(N1))
437
145
CONTINUE
438
PRINT 133,K,SK(K),T,KTP,I8RP,T4
439
DO 143 L=1,NJ,16
440
LK=L*15
441
IF (LK.GT.NJ) LK=NJ
442
PRINT 134,(J,J=L,LK)
443
PRINT 135,(R(J),J=L,LK)
444
DO 143 I=1,NI
445
PRINT 136,(I,Z(I),(UN(I,J,K),J=L,LK))
446
143
CONTINUE
447
IF (JP(1).EO.0) GO TO 165
448
DO 165 LP=1,JN
449
J=JP(LP)
450
DO 148 K=1,NK
451
DO 148 I=1,NI
452
N1=N(I,J,K)
453
HN1=H(I,J,K)
454
IF ((H1(N1).LE.HN1).AND.(HN1.LE.H2(N1)))
455
1
UN(I,J,K)=700.0+100*(HN1-H1(N1))/(H2(N1)-H1(N1))
456
148
CONTINUE
457
PRINT 160,J,R(J),T,KTP,I8RP,T4
458
DO 165 L=1,NK,16
459
LK=L+15
460
IF (LK.GT.NK) LK=NK
461
PRINT 161,(K,K=L,LK)
462
PRINT 135,(SK(I),I=L,LK)
463
DO 166 I=1,N1
464
166
PRINT 136,(I,Z(I),(UN(I,J,K),K=L,LK))
465
165
CONTINUE
466
IF (IF(1).EO.0) GO TO 167
467
DO 167 LP=1,IW
468
I=IP(LP)
469
DO 149 K=1,NK
470
DO 149 J=1,NJ
471
N1=N(I,J,K)
472
HN1=H(I,J,K)
473
IF ((H1(N1).LE.HN1).AND.(HN1.LE.H2(N1)))
474
1
UN(I,J,K)=700.0+100*(HN1-H1(N1))/(H2(N1)-H1(N1))
475
149
CONTINUE
476
PRINT 162,I,Z(I),T,KTP,I8RP,T4
477
DO 167 L=1,NK,16
478
LK=L+15
479
IF (LK.GT.NK) LK=NK
480
PRINT 161,(K,K=L,LK)
481
PRINT 163,(SK(J),J=L,LK)
482
DO 168 J=1,NJ
483
168
PRINT 136,(J,Z(J),(UN(I,J,K) K=L,LK))
484
167
CONTINUE
485
GO TO 10
486
15
CONTINUE
487
RETURN
488
END
*** WNS ***
1
SUBROLTINE WNS(N,A,B,C)
2
DIMENSION A(N),B(N),C(N),T(3,10),V(3,10),Z(10,11),U(11)
3
11
FORMAT (12F6,2)
4
READ 11,((T(I,J),I=1,3),J=1,’)
5
READ 11,((T(I,J),I=1,3),J=1,’)
7
8
FORMAT (22l,7l¢,(11)=,F6.2,3X,ZHUR(T2)=,F5.2,3X,
ZHUR(T3)=,F6.2)
8
PRINT 16
13
U(1)=-30.0
14
DO 6 J=2,6
15
6
U(J)=U(J-1)+5.0
16
DO 5 J=7,11
17
5
U(J)=U(J-1)+1.0
18
DO 10 I=1,4
19
13
W2=W(2,I)-W(3,I)
20
U1=T(1,I)-T(3,I)
21
W1=W(1,I)-W(3,I)
22
U2=T(2,I)-T(3,I)
23
IF (((-0.0001.LT.W2).AND.(W2.LT.0.0001)).OR.((-0.0001.LT.
24
#
W1).AND(W1.LT.0.0001))) GO TO 2
25
GO TO 24
26
2
A(I)=0.0
27
B(I)=0.0001
28
C(I)=W(1,I)
29
GO TO 4
30
74
B(I)=(T(2,I*W2+U1-T(1,I)+W1*U2)/(W2*U1-W1*U2)
31
A(I)=-1*(B(I)-T(1,I))*(B(I)-T(((3,I))/U1
32
C(I)=-(3,I)-A(I)/B(I)-T(3,I)
41
4
DO 3 J=1,11
42
3
Z(I,J)=A(I)/B(I)-U(J))+C(I)
43
10
CONTINUE
44
16
FORMAT (T41,’НЕЗАМЕРЗШАЯ ВОДА ПО СЛОЯМ,)
45*
17
FORMAT (2H0,T25,’ВХОДНЫЕ ДАННЫЕ’,30X,
46
#
‘ВЫЧИСЛЕННЫЕ КОЭФФИЦИНТЫ <PHBO’)
47
18
FORMAT (1H),11X,4HСЛОЙ,I3,6X,4HT1=.F6.2,6X,4HT2=
48
#
,F6.2X,4HT3,2HA(,12,2H)=,F6.4,4X,2HB(,IZ,ZI)=,F8.4)
49
19
FORMAT (22X,7HW(T1)=,F6.2,3X,7HW(T2)=,F62,3X,7HW
50
#
(T3)=,F6.2,3X,2HA(,I2,2H)=,F8.4,4X,2HB(,I2,2H)=,F8.4)
51
PRINT 17
52
DO 20 I=1,4
53
PRINT 18,(I,T(1,I),T(2,I),T(3,I))
54
20
PRINT 19,(W(1,I),W(2,I),W(3,I),I,A(I),I,B(I),I,C(I))
55
21
FORMAT (2H0,T40,’КОНТРОЛЬНАЯ ТАБЛИЦА ЗНАЧЕНИЙ WN ОТ Т’)
56
22
FORMAT (1H0,14X,2HT=,7X,11F7,2)
57
23
FORMAT ((11X,4НСЛОЙ,I3,6X,11F7.2)
58
PRINT 21
59
PRINT 22,U
60
PRINT 23,((I,(Z(I,J),J=1,11)),I=1,N)
61
RETURN
62
END
*** FIND ***
1
SOHROLTINL FIND(HN1,N1,UL1)
11
COMMO\ /TABL/ FT(50,10),FS10(1),F25(15,10),F02(20,
21
#
10),H1(10),H2(10)
31
IF (HN1.GT/IT(30,N1) GO TO 201
41
IF (HN1.GT/IT(40,N1) GO TO 202
51
ML=41
61
ML=50
71
GO TO 203
81
202
ML=31
91
ML=40
101
203
DO 204 IN=ML,MU
111
IF (HN1.GT/IT(IN,N1) GO TO 205
121
204
CONTINUE
131
205
WA=FT(IN,N1)
141
WB=FT(IN-1,N1)
151
XA=IN
161
XB=IN-1
171
GO TO 200
181
201
IF (HN1.GT/IT(10,N1) GO TO 206
191
IF (HN1.GT/IT(20,N1) GO TO 207
201
ML=21
211
ML=30
221
GO TO 203
231
207
ML=11
241
ML=20
251
GO TO 203
261
208
IF (HN1.GT/IT(2,N1) GO TO 208
271
IF (HN1.GT/IT(5,N1) GO TO 209
281
DO 210 IN=2,10
291
IF (HN1.GE.F510(IN,N1)) GO TO 211
301
210
CONTINUE
311
211
WA=F510(IN,N1)
321
WB=F510(IN-1,N1)
331
XA=-5.0-0.5*IN
341
XB=XA+0.5
351
GO TO 200
361
209
DO 212 IN=2,15
371
IF (HN1.GE.F25(IN,N1)) GO TO 213
381
212
CONTINUE
391
213
WA=F25(IN,N1)
401
WB=F25(IN-1,N1)
411
XA=-2.0-0.2*IN
421
XB=XA+0,2
431
GO TO 200
441
208
IF (HN1.GE.F02(1,N1)) GO TO 216
451
GO TO 217
461
216
WA=F02(1,N1)
471
WB=N1(N1)
481
XA=-0.1
491
XB=0.0
501
GO TO 200
511
217
DO 214 IN=2,20
521
IF (HN1.GE.F02(IN,N1)) GO TO 213
531
214
CONTINUE
541
215
WA=F02(IN,N1)
551
WB=F02(IN-1,N1)
561
XA=-0.1+IN
571
XB=XA+0,1
581
200
UL1=(HN1-WA)/(WB-WA)+(XB-XA)+XA
591
RETURN
601
END
*** HF ***
1
FUNCTION HF(N,U1)
11
COMMO\ /RFIS/BCK(10),CCK(10),WC(10),UKR(10)
21
COMMO\ /CUEI/A(10),B(10),C(10)
22
COMMO\ /CONST/RV,CL,CB,TFR,TK
31
D=B(N)
41
D1=TK-J
51
D2=U1-TK
61
HF=BCK(N)*((CCK(N)+RV*CL+WC(A)+C(B)+(CB-R/*CL))*
71
3
D1+\((RV*CL-CB)*ALOG(ABC(U1-D)/C1))*D2/((J1-D)*D1)))
81
RETURN
91
END
Решается задача в 3 измерениях
направление осей К
0----J
/
/
/
/
I
Шаг по времени в часах 121.67
Граница измерений 876000.00
Число блоков по I 7
Номер блока
1
2
3
4
5
6
7
Размер блока
2.00
2.00
2.00
2.00
3.00
3.00
3.00
Коорд. Центра
1.00
3.00
5.00
7.00
9.50
12.50
15.50
Число блоков по J 8
Номер блока
1
2
3
4
5
6
7
8
Размер блока
3.00
3.00
2.00
2.00
2.00
3.00
3.00
3.00
Коорд. центра
1.50
4.50
7.00
9.00
11.00
13.50
16.50
19.50
Число блоков по K 6
Номер блока
1
2
3
4
5
6
Размер блока
3.00
3.00
3.00
2.00
2.00
2.00
Коорд. центра
1.50
4.50
7.50
10.00
12.00
14.00
Условия на границах
K/J
1
2
3
4
5
6
7
8
6
121
121
121
121
111
111
111
111
5
121
121
121
121
131
111
111
111
4
131
131
131
131
131
131
111
111
3
111
111
131
131
131
131
131
131
2
111
111
111
111
111
111
131
131
1
111
111
111
111
111
111
111
111
Шифр XYZW если X=0 - нет источника
Y - номер типа источника
Z - номер зоны краевых условий на верхней границе
W - тоже на нижней
На верхней границе
Номер зоны / месяц
1
2
3
4
5
6
7
8
9
10
11
12
1
темп. возд.
-4.00
-23.20
-19.20
-9.50
3.70
13.10
10.00
12.40
4.60
-4.60
16.40
-22.40
сопротивл.
2.00
2.70
2.80
3.00
1.70
0.0
0.0
0.0
0.0
1.30
2.10
2.30
теплоотд.
Побачили розбіжність з офіційним текстом?
Дивіться також
СанПиН 2.3.2.1293-03
Гигиенические требования по применению пищевых добавок
СНУ Збірник 1