*** 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 (22??,7,(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 |
