Пример обработки полученных результатов приведен на рис. 3 (вклейка).

 

 

Рис. 3. Результаты расчета температурного режима грунта по программе PROGNOZ

 

Исходные данные для расчета контрольного примера приложения 1 по программе PROGNOZ-3S

 

Перфорировал _______ ___                 Проверил_____________

                                                                перфорацию

 

Перфорировал _______ ___                 Проверил_____________

                                                                перфорацию

 

Перфорировал _______ ___                 Проверил_____________

                                                                перфорацию

 

Выдача на печать для контроля ввода начальных данных подготовленных для расчета по программе PROGNOZ-3S контрольного примера приложения 1

 

-0.3

-1.0

-10.0

-0.3

-1.0

-10.0

0.0

0.0

0.0

 

 

 

 

0.14

0.12

0.08

0.12

0.08

0.05

0.0

0.0

0.0

 

 

 

 

1.30

1.55

2.35

1.15

1.45

2.15

2.0

2.0

2.0

2.0

3.0

3.0

3.0

3.0

3.0

2.0

2.0

2.0

3.0

3.0

3.0

3.0

3.0

3.0

2.0

2.0

2.0

121.67

 

 

 

 

 

 

 

 

 

 

 

876000.0

 

 

 

 

 

 

 

 

 

 

 

-24.6

-23.2

-19.2

-9.5

3.7

13.1

18.0

12.4

4.6

-4.6

-16.4

-22.4

 

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

 

-24.6

-23.2

-19.2

-9.5

3.7

13.1

18.0

12.4

4.6

-4.6

-16.4

-22.4

 

2.6

2.7

2.8

3.0

1.7

0.0

0.0

0.0

0.0

1.3

2.1

2.3

 

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

1.8

 

4.4

4.4

4.7

5.0

2.82

0.0

0.0

0.0

0.0

2.14

3.57

3.59

 

13.8

12.3

12.3

12.0

12.1

12.6

13.4

13.0

14.3

13.4

13.5

12.7

 

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

20.0

 

13.8

12.3

12.3

12.0

12.1

12.6

13.4

13.0

14.3

13.4

13.5

12.7

 

-1.4

0.0

0.0

0.0

0.0

 

 

 

 

 

 

 

 

3

3

3

1

2

2

2

2

1

4

6

8

1

4

6

0

 

 

 

 

 

 

/*

 

 

 

 

 

 

 

 

 

 

 

 

//

 

 

 

 

 

 

 

 

 

 

 

 

Приложение 2

Рекомендуемое

Пример текста подпрограммы ISTO

 

С                ПОДПРОГРАММА ISTO

С       ОПРЕДЕЛЯЕМ КОЛИЧЕСТВО ВЫДЕЛЕННОЙ ИЛИ

С       ПОГЛОЩЕННОЙ F В КАЛЛ (ИЛИ КДЖ) В ЕДИНИЦУ

С       ВРЕМЕНИ (ЧАС) В ФИКСИРОВАННОМ ЭЛЕМЕНТЕ С

С       КООРДИНАТАМИ (I,J,K)

С       Т - ТЕКУЩЕЕ ВРЕМЯ ОТ НАЧАЛА РАСЧЕТА

SUBROUTINE ISTO (I,J,K,T,F)

IF (T.GT.152.0).AND.(T.(T.154.0)) GO TO 10

GO TO 11

10 IF ((((K.EQ.3).AND.((J.EQ.3).OR.(J/EQ.5))).OR((J.EQ.4).AND.

# (K.EQ.5))).AND.(I.LE.7)) GO TO 12

GO TO 11

12 IF (I.LE.6) F=16666.6

IF (I.EQ.7) F=25000.0

11 OCNTINUE

RETURN

END

Приложение 3

Обязательное

Текст программы PROGNOZ-3S с результатами расчета контрольного примера приложения 1

 

Раздел MAIN-3S содержит подпрограммы MAIN и BLOCK DATA составленные для решения контрольного примера приложения 1.

Вызов программы на счет из библиотеки исходных модулей имеет вид??

// JOB

// CLASS=C, MSGLEVEL=1, TIME=15

// EXEC FORTHCLG, PARM, FORT=(OPT=02¢),

// REGION, FORT=300K, REGION, LKED=100 K, REGION,

// GO=100 K

// FORT. SYSIN DD DSN=BIBL 4,S (MAIN3S),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL 4,S (GRUN3S),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL4,S (WNS),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL4,S (FIND),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

//                DD DSN=BIBL4,S (HF),

//                          VOL=SER=DTN033, DISP=SHR, UNIT=5061

Здесь  BIBL4,S - имя библиотеки исходных модулей

BIBL4,S - имя тома.

Результаты расчета контрольного примера представлены на 1-й и 5-й расчетный год.

*** MAIN3S ***

1

 

 

DIBENSION H(7,8,6),UR(7,8,6),HRDV(7,8,6),J(7,8,6),

2

 

#

A1(3),B1(3),C1(3)

3

 

 

IAIEGER NN(7,8,6)/336*3/XL(3,6)/14*111,2*131,2*111,12

4

 

#

*131,2*114*121/131,3*111,4*121,4*111/

5

 

 

REAL HPRI(2/-1.0,8760.0/

6

 

 

CALLWNS(3,A1,B1,C1)

7

 

 

TKUR=17601.0

8

 

 

DO 1 K=1,8

9

 

 

DO 1 J=1,8

10

 

1

N(1,J,K)=1

11

 

 

DO 2 K=1,3

12

 

 

DO 2 J=3,8

13

 

 

DO 2 I=2,5

14

 

2

N(I,J,K)=2

15

 

 

DO 3 K=1,5

16

 

 

DO 3 J=1,8

17

 

 

U(1,J,K)=-0,1

18

 

 

U(2,J,K)=-0,4

19

 

 

U(3,J,K)=-1,0

20

 

 

U(4,J,K)=-1,2

21

 

 

U(5,J,K)=-1,4

22

 

 

U(6,J,K)=-1,4

23

 

3

U(7,J,K)=-1,4

24

 

 

CALL GRAND(1,7,5,6,7,12,3,1,TKDH,A1,B1,C1,H,KL,U,

25

 

#

1,4,3,2,HPRI,DV,HRDV

26

 

 

STOP

27

 

 

END

 

C

 

 

 

C

 

 

 

C

 

 

 

C

 

 

 

C

 

 

28

 

 

BLOCK DATA

29

 

 

COMMON/RBIS/BCK(10),CCK(10),WC(10),UKR(10)

30

 

 

DATA BCK/1390,0.1520,0.1500.0,7*0/. CCK/.22,0.24,0.23

31

 

#

7*0/,WC/0.25,0.27,7*0.0/,UKR/10*0.0/

32

 

 

END

*** GRAN3S ***

1

 

 

SLBROLTINE GRUND(SL,NI,NJ,NK,NT,LI,M1,IK,AX,KL

2

 

#

,IW,JW,KW,HPRI,H,HNDV,SI)

3

 

 

DIMENSION H(NI,NJ,NK),Z(100),R(100),SK(100),

4

 

1

CIAL(10),CI(10), TLN(10),TLT(10)HI(100),HJ(100)HK(100)

5

 

2

TR(10),U(NI,NJ,NK)TB(36,10)RB(36,10)AB(36,10)

6

 

3

KL(NJ,NK),NJ(10),NK(10),JP(100),IP(100),KP(100)

7

 

4

AX(A,SL),BX(B,SL),CX(C,SL),UN(NI,NJ,NK)

8

 

 

COMMO /BCK(10),CCK(10),WC(10),HKP(10)

9

 

 

COMMO /A(10),B(10),C(10)

10

 

 

COMMO /RV,CL,CB,TFP

11

 

 

COMMO /FT(50,10),F510(10,10),F25(15,10),F02(20,15),

12

 

#

H1(10),H2(10)

13

 

 

RV=1

14

 

 

CL=2.0234

15

 

 

CB=4.0368

16

 

 

TFP=337.944

 

 

 

PP=3.6

17

 

 

IF (RSI.ET.1) GO TO 1

18

 

 

CL=0.5

19

 

 

CB=1.0

20

 

 

PP=1.0

21

1

 

TK=273.0

22

 

 

TPER=3760.0

23

 

 

DO 180 I=1,HSL

24

 

 

A(I)=AX(I)

25

 

 

B(I)=BX(I)

26

180

 

C(I)=CX(I)

27

 

 

DO 216 K=1,HSL

28

 

 

DO 217 IN=1,SJ

29

 

 

UL1=-1

30

217

 

FI(IN,K)=NF(K,UL1)

31

 

 

DO 218 IN=1,10

32

 

 

UL1=-5.0-0.5*IN

33

218

 

F510(I1,K)=IF(K,UL1)

34

 

 

DO 219 IN=1,15

35

 

 

UL1=-2.0-0.2*IN

36

219

 

F25(IN,K)=HI(K,UL1)

37

 

 

DO 220 IN=1,20

38

 

 

UL1=-0.1*IN

39

220

 

FCZ(IN,K)=IN(K,UL1)

40

216

 

CONTINUE

41

181

 

FORMAT (13F5.2)

42

182

 

FORMAT (12F6.2)

43

183

 

FORMAT (6F12.2)

44

184

 

FORMAT (27IN)

45

 

 

READ 131,(TLH(I),I=1,NSL),(TLT(I),I=1,NSL),(HI(I),I=1,NI),

46

 

#

(HJ(I),II=1,NJ),(HK(I),I=1,NK),(SLT(I),I,NI)

47

 

 

READ 133,(TR(I),I=1,NI)

48

 

 

READ 132, ((IB(I,J),I=1,LI),J=1,N),((RF(I,J),I=1,LI),J=1,N),

49

 

#

((AB(I,J),I=1,LI),J=1,N),(TBIK(I),I=1,N1),TRJ,TBJK,TBK,TB

50)

 

 

READ 134,(HR(I),I=1,N),(NK(I),I=1,5),(JP(I),J=1,JP),(KP(I),

51

 

#

I=1,KP),(IP(I),I=1,IP)

52

 

 

T=0.0

53

 

 

IT=1

54

 

 

KT=1

55

 

 

T2=TPER/LI

56

 

 

HI=SHI(IT)

57

133

 

FORMAT (5X,2IK=,I3,2X,F7.2,1H,2X,I=,F2.2,12X,1H=,

58

 

1

2X,I3,7YПЕРИОДАБ,I3,9H’ИНТЕГРАЛ’F7.2,4HЧАСА)

59

134

 

FORMAT (2X,12HБЛОКИ,15,16,17)

60

135

 

FORMAT (2X,12HМЕТРЫ,13,14,15)

61

136

 

FORMAT (2X,15,17F7.2)

62

160

 

FORMAT (5A,HJ=,13,2X,17F.2

63

 

 

IH=,2X,I3,7HПЕРИОДА,I3,9H’ИНТЕГРАЛА’F7.2,4НЧАСА)

64

161

1

FORMAT (2X,H2’,БЛОК’,15,16I7)

65

162

 

FORMAT (5X,2HПЕРИОДА,I3,9HИНТЕРВАЛА,

66

 

 

F7.2,4HЧАСА)

67

163

 

FORMAT (2X,H2МЕТРЫ,13F3.2)

68

 

 

DO 5 I=1,NSL

69

 

 

H1(I)=HF(2,UKR(I))-15.0

70

 

 

H2(I)=HF(1)+TFR*(CCK(I)*BCK(I)*(WC(I)-A(I)/(B(I)-UKP

 

 

 

(I))-C(I))+15.0

71

 

 

CIAL(I)=BCK(I)*(CCK(I)+WC(I)+CB)

72

5

 

CW(I)=CCK(I)+CL*HC(I)+C(I)*(CB-RV*CL)

73

 

 

DO 5 K=1,K

74

 

 

DO 6 J=1,J

75

 

 

DO 5 I=1,I

76

 

 

DV=NI(I)*NJ(J)*NK(K)

77

 

 

N1=N(I,J,K)

78

 

 

UL=(I,J,K)

79

 

 

IF (UL.LT.(NKR(N1)-0.0001)) GO TO 7

80

 

 

IF (UL.LT.(NKR(N1)+0.0001)) GO TO 9

81

 

 

GO TO 3

82

 

 

H(I,J,K)=HF(N1,UL)

83

 

 

GO TO 3

84

 

 

H(I,J,K)=(O.5*(H1(N1)+H2(N1))

85

 

 

GO TO 3

86

 

 

H(I,J,K)=(H(N1)+CTIL(N1)*(UL-UKZ(N1)))

87

 

 

CONTINUE

88

 

 

IPR=Z

89

 

 

IPR=HPRI(Z)

90

 

 

SK(1)=HK(J)

91

 

 

IF (NK.EO.1) GO TO 175

92

 

 

DO 132 K=2,NK

93

132

 

SN(K)=SK(K)

94

175

 

DO 138 K=1,NK

95

178

 

SN(K)=SK(K)-0.5*HK(K)

96

 

 

R(1)=KJ(1)

97

 

 

IF (NJ.EO.1,GO TO 176

98

 

 

DO 139 J=2,HJ

99

130

 

R(J)=R(J-1)*HJ(J)

100

176

 

DO 140 J=1,HJ

101

140

 

R(J)=R(J)-0.5*HJ(J)

102

 

 

Z(1)=HI(1)

103

 

 

IF (NI.EO.1) GO TO 177

104

 

 

DO 141 I=2,HI

105

141

 

Z(I)=Z(I-1)*HI(I)

106

172

 

DO 142 I=1,HI

107

142

 

Z(I)=Z(I)-0.5*HI(I)

108

11

 

FORMAT (20X,17HРЕШАЕТСЯ ЗАДАЧА 313??11Н

 

 

 

ИЗММЕРЕНИЯХ)

109

12

 

FORMAT (15X,16HНАПРАВЛЕНИЕ ОСЕЙ6Х??1НК)

110

13

 

FORMAT (140,30Х6Н0----J)

111

14

 

FORMAT (31X,1H/)

112

22

 

FORMAT (30X,1HI)

113

23

 

FORMAT (10X,22HШАГ ПО ВРЕМЕНИ В ЧАСАХ??

 

 

 

(8F11.2))

114

24

 

FORMAT (12X,17HГРАНИЦА ИЗМЕРЕНИЯ5Х, (8F11.2))

115

25

 

FORMAT (140,9X,17HЧИСЛО БЛОКОВ ПО II4)

116

26

 

FORMAT (2X,12HНОМЕР БЛОКАI4,15I7)

117

27

 

FORMAT (2X,12HРАЗМЕР БЛОКА16F7.2)

118

28

 

FORMAT (2X,12HКООРД. ЦЕНТРА,16F7.2)

119

29

 

FORMAT (1409X,17HЧИСЛО БЛОКОВ ПО J,I4)

120

30

 

FORMAT (140,9X,17HЧИСЛО БЛОКОВ ПО KI4)

121

33

 

FORMAT (140,14X,21HУСЛОВИЯ НА ГРАНИЦАХ)

122

34

 

FORMAT (3X,3HK/J,I4,23I5)

123

35

 

FORMAT (2X,I424I5)

124

36

 

FORMAT (5X,32HШИФР XYZW ЕСЛИ X=0 - НЕТ

 

 

 

ИСТОЧНИКА)

125

37

 

FORMAT (20X,25HY - НОМЕР ТИПА ИСТОЧНИКА)

126

38

 

FORMAT (20X,49HZ - НОМЕР ЗОНЫ КРАЕВЫХ

 

 

 

УСЛОВИЙ НА ВЕРХНЕЙ ГРАНИЦЕ)

127

39

 

FORMAT (20X,18H - ТОЖЕ НА НИЖНЕЙ)

128

40

 

FORMAT (140,21X,18H - НА ВЕРХНЕЙ ГРАНИЦЕ)

129

41

 

FORMAT (2X,19HНОМЕР ЗОНЫ/МЕСЯЦ,I4,12I7)

130

43

 

FORMAT (2X,16,4X,5HТЕМП ВОЗД.,15F7.2)

131

46

 

FORMAT (12X,9HСОПРОТИВЛ.,15F7.2)

132

47

 

FORMAT (12X,9HТЕПЛООТД.,15F7.2)

133

48

 

FORMAT (2X,24HПРОДОЛЖИТЕЛЬНОСТЬ МЕСЯЦА,

 

 

 

F9.2,5H ЧАСА)

134

49

 

FORMAT (2X,24HПЕРИОДИЧНОСТЬ ИЗМЕРЕНИЯ,

 

 

 

F9.2,5H ЧАСА)

135

121

 

FORMAT (140,14X,29HТЕМПЕРАТУРА НА ВЕРХНЕЙ

 

 

 

ГРАНИЦЕ)

136

122

 

FORMAT (5X,10HНОМЕР ЗОНЫ,1415I7)

137

123

 

FORMAT (5X,10HТЕМПЕРАТУРА,15F7.2)

138

124

 

FORMAT (15X,10H)

139

125

 

FORMAT (5X,10HTBJ=,F7.2,2X,5HTBJ<=,F7.2)

140

 

1

1F7.2

141

126

 

FORMAT (10X,11HНОМЕР СЛОЯ2X)

142

 

*

25HТЕПЛОПРОВ??3X)

143

127

 

FORMAT (5X,15F12.2,5F7.23F7.3)

144

128

 

FORMAT (14015X,28HНАЧАЛЬНАЯ ТЕМПЕРАТУРА

 

 

 

ГРУНТА)

145

130

 

FORMAT (15X,16I7)

146

131

 

FORMAT (5X,25HШАГ ГЕНЕРАТОРА РЕЗУЛЬТАТА

 

 

 

(10F9.2))

147

132

 

FORMAT (5X,2HK=,13,2X,F7.2)

148

164

 

FORMAT (5X,2HJ=,13,2X,F7.2)

149

110

 

FORMAT (5X,2HI=,13,2X,F7.2)

150

 

 

L=2

151

 

 

IF ((NK.ST.1,ND.(NJ.ST.1)) GO TO 150

152

 

 

IF ((NK.EO.1,ND.(NJ.EO.1)) L=1

153

 

 

DO TO 151

154

150

 

L=3

155

151

 

PRINT 11,L

156

 

 

PRINT 12

157

 

 

PRINT 13

158

 

 

DO 152 I=I,4

159

152

 

PRINT 14

160

 

 

PRINT 22

161

 

 

PRINT 23,(SHT(I),I=1,NI)

162

 

 

PRINT 24,(TR(I),I=1,NT)

163

 

 

PRINT 25I

164

 

 

DO 153 L=1,NI,15

165

 

 

LK=L+14

166

 

 

IF (LK.ST.(I)LK=NI)

167

 

 

PRINT 26,(I,I=L,LK)

168

 

 

PRINT 27,(I,I=L,LK)

169

153

 

PRINT 28,(Z(I),I=L,LK)

170

 

 

PRINT 29,J

171

 

 

DO 154 L=1,HJ,15

172

 

 

LK=L+14

173

 

 

IF (LK.ST.J)LK=NJ)

174

 

 

PRINT 26,(J,J=L,LK)

175

 

 

PRINT 27,(J,J=L,LK)

176

154

 

PRINT 28,(R(J),J=L,LK)

177

 

 

PRINT 30,K

178

 

 

DO 155 L=1,NK,15

179

 

 

LK=L+14

180

 

 

IF (LK.ST.NK) LK=NK

181

 

 

PRINT 26,(K,K=L,LK)

182

 

 

PRINT 27,(NK(K),K=L,LK)

183

155

 

PRINT 28,(SK(K),K=L,LK)

184

 

 

PRINT 33

185

 

 

DO 157 L=1,NJ,24

186

 

 

LK=L+23

187

 

 

IF (LK.GT.NJ) LK=NJ

188

 

 

PRINT 34,(J,J=L,LK)

189

 

 

DO 157 K=1,NK

190

 

 

I=NK=K+1

191

157

 

PRINT 35,(I,(KL(J,I),J=L,LK)

192

 

 

PRINT 36

193

 

 

PRINT 37

194

 

 

PRINT 38

195

 

 

PRINT 39

196

 

 

PRINT 40

197

 

 

DO 154 K=1,NK

198

 

 

PRINT 41,(I,I=1,LI)

199

 

 

PRINT 43,(I,TB(I,K),I=1,LI)

200

 

 

PRINT 46,(I,RB(I,K),I=1,LI)

201

158

 

PRINT 47,(I,AB(I,K),I=1,LI)

202

 

 

PRINT 48,TZ

203

 

 

PRINT 49,TPER

204

 

 

PRINT 121

205

 

 

PRINT 122,(I,I=1,MI)

206

 

 

PRINT 123,(TBIK(I),I=1,MI)

207

 

 

PRINT 124

208

 

 

PRINT 125TBJ,TBJK,TBK,TBKK

209

 

 

PRINT 126

210

 

 

DO 159 K=1,NSL

211

 

 

PRINT 127,K,BCK(K),CCK(K),WC(K),UKR(K),TLM(K),

212

 

2

TLT(K),(K),B(K),C(K)

213

159

 

CONTINUE

214

 

 

PRINT 129

215

 

 

IF (KP(1).EO,J) GO TO 160

216

 

 

DO 170 LP=1,KP

217

 

 

K=KP(JP)

218

 

 

PRINT 132,K,SK(K)

219

 

 

DO 170 L=1,HJ,16

220

 

 

LK=L+15

221

 

 

IF (LK.GT.NJ) LK=NJ

222

 

 

PRINT 134,(J,J=L,LK)

223

 

 

PRINT 135,(R(J),J=L,LK)

224

 

 

DO 170 I=1,NI

225

 

 

PRINT 136,(I,Z(I),(U(I,J,K),J=L,LK))

226

170

 

PRINT 130,(H(I,J,K),J=L,LK)

227

160

 

IF (JP(1).EO.0) GO TO 171

228

 

 

DO 172 LP=1,JN

229

 

 

J=JN(LP)

230

 

 

PRINT 164,J,R(J)

231

 

 

DO 172 L=1,NK,16

232

 

 

LK=L+15

233

 

 

IF (LK.GT.NK) LK=NK

234

 

 

PRINT 161,(K,K=L,LK)

235

 

 

PRINT 135,(SK(I),I=L,LK)

236

 

 

D0 172 I=1,NI

237

 

 

PRINT 136,(I,Z(I),(U(I,J,K),K=L,LK))

238

172

 

PRINT 130,(H(I,J,K),K=L,LK)

239

171

 

IF (IP(1).EO.O) GO TO 173

240

 

 

D0 174 LP=1,IP

241

 

 

I=IP(LP)

242

 

 

PRINT 110,I,Z(I)

243

 

 

D0 174 L=1,NK,16

244

 

 

LK=L+15

245

 

 

IF (LK.GT.NK) LK=NK

246

 

 

PRINT 161,(K,K=L,LK)

247

 

 

PRINT 163,(SK(J),J=L,LK)

248

 

 

DO 174 J=1,NJ

249

 

 

PRINT 136,(J,R(J),(U(I,J,K),K=L,LK))

250

174

 

PRINT 130,(H(I,J,K),K=L,LK)

251

173

 

PRINT 131,(PRHI(I),I=2,INT)

252

10

 

T=T+NT

253

201

 

CONTINUE

254

 

 

IF (T.GT.TKON) GOTO 15

255

 

 

IF ((1-KT*TPC).GT.0.0) KT=KT+1

256

 

 

IF (T-TR(IT) 16,17,17

257

17

 

IT=IT+1

258

 

 

HI=SHI(IT)

259

16

 

T3=T-(KT-1).TPER

260

 

 

IBR=1

261

18

 

IF (T3-IBR*T2) 19,19,20

262

20

 

IBR=IBR+1

263

 

 

GO TO 13

264

19

 

T4=T3-(IBR-1)*T2

265

21

 

D0 31 K=1,NK

266

 

 

D0 32 J=1,NJ

267

 

 

IK=KL(J,K)

268

 

 

MC=IK/10-IK/100*10

269

 

 

MCK=IK/10*10

270

 

 

DO 50 I=1,NI

271

 

 

H1=,(I,J,K)

272

 

 

HN=H(I,J,K)

273

 

 

UL=U(I,J,K)

274

 

 

IF (HN-H1(NI) 51,52,52

275

51

 

TL=TLX(N1)

276

 

 

GO TO 53

277

52

 

IF (HN-H2(NI) 45,45,44

278

45

 

TL=0.5*(TLX(N1)+TLT(N1)

279

 

 

GO TO 53

280

44

 

TL=TLT(N1)

281

53

 

IF (I.GE.0) GO TO 54

282

55

 

GO TO 56

283

 

 

TE10=IB(IBR-1,MC)

284

 

 

GO TO 57

285

56

 

TE10=IB(LI,MC)

286

57

 

U1=TB10+(TB(IBR,MC)-TB10)/T2*T4

287

 

 

IJ=BR(10)

288

 

 

GO TO (58,59,60),IJ

289

58

 

R1=TL/TI(1)*(J1-CL)*2

290

 

 

GO TO 66

291

59

 

R1=U1

292

 

 

GO TO 66

293

60

 

GO TO (61),IBR

294

 

 

RB10=RB(IBR-1,MC)

295

 

 

AB10=AB(IBR-1,MC)

296

 

 

GO TO 62

297

61

 

RB10=RB(LI,MD)

298

 

 

AB10=AB(LI,MD)

299

62

 

RB11=RB10+(RB(IBR,MD)-RB10)/T2*T4

300

 

 

AB11=AB10+(AB(IBR,MD)-AB10)/T2*T4

301

 

 

R1=2/(TI(1)/TL+2*(1/AB11+RB11))*(H1-UL)

302

 

 

GO TO 66

303

54

 

U1=U(I-1,J,K)

304

 

 

N2=H(I-1,J,K)

305

 

 

HN2=H(I-1,J,K)

306

 

 

IF (HN2-H1(N2)) 63,64,64

307

63

 

TL1=TLT(N2)

308

 

 

GO TO 65

309

64

 

TL1=TLT(N2)

310

 

 

IF (HN2.LS.H2(N2)) TL1=0.5*(TLN(N2)+TL1)

311

65

 

UX=HI(I-1)+HI(I)

312

 

 

R1=2*(HI(I-1)*TL1+HI(I)*TL)/)UX+LX)*(H1-UL)

313

66

 

IF (I-I1) 67,68,67

314

67

 

U2=U(I+1,J,K)

315

 

 

H2=H(I+1,J,K)

316

 

 

HN2=H(I+1,J,K)

317

 

 

IF (HN3-H1(N2)) 69,70,70

318

69

 

TLZ=TLH(N2)

319

 

 

GO TO 71

320

70

 

TLZ=TLT(N2)

321

 

 

IF (HNR.LE.12(N2)) TLZ=0.5*(TLN(N2)+TL2

322

71

 

UX=HI(I+1)+HI(I)

323

 

 

R2=2*(HI(I+1)*TLZ+HI(I)*TL)/(UX*LX)*(U2-UL)

324

 

 

GO TO 72

325

68

 

U2=UB.K(M,X)

326

 

 

IF (MX(1)-1) 73,74,73

327

74

 

R2=TL/HI(I)*(U2-UL)*2

328

 

 

GO TO 72

329

73

 

R2=U2

330

72

 

IF (J.UE.1) GO TO 75

331

 

 

U1=TBJ

332

 

 

IF (MK(2)-1) 76,77,76

333

77

 

R3=TL/HJ(1)-(U3-UL)*2

334

 

 

GO TO 70

335

76

 

R3=TBJ

336

 

 

GO TO 70

337

75

 

UJ=U(I,J-1,K)

338

 

 

N2=N(I,J-1,K)

339

 

 

HN2=H(I,J-1,K)

340

 

 

IF (HN2-H1(N2)) 80,81,81

341

80

 

TL3=TLT(N2)

342

 

 

GO TO 82

343

81

 

TL3=TLT(N2)

344

 

 

IF (HN2.LE.12(N2)) TL3=0.5*(TLM(N2)+TL3)

345

82

 

UX=HJ(J-1)+HJ(J)

346

 

 

R3=2*(HJ(J-1)*TL3+HJ(J)*TL)/(UX-LX)*(U3-UL)

347

79

 

IF (J.LE.HJ) GO TO 83

348

 

 

U4=TBJK

349

 

 

IF (MX(3)-1) 84,85,84

350

85

 

R4=TL/HJ(NJ)*(U4-UL)*2

351

 

 

GO TO 90

352

84

 

R4=TBJK

353

 

 

GO TO 90

354

83

 

U4=U(I,J+1,K)

355

 

 

N2=N(I,J+1,K)

356

 

 

HN2=H(I,J+1,K)

357

 

 

IF (H/2-H1(N2)) 87,88,88

358

87

 

TL4=TLT(N2)

359

 

 

GO TO 89

360

88

 

TL4=TLT(N2)

361

 

 

IF (HN2.LE.H2(N2)) TL4=0.5*(TLM(N2)+TL4)

362

89

 

UX=HJ(J+1)+HJ(J)

363

 

 

R4=2*(HJ(J+1)*TL4+HJ(J)*TL)/(UX*UX)*(U4-UL)

364

90

 

IFF (K.NE.1) GO TO 91

365

 

 

U5=TBK

366

 

 

IF (MX(4)-1) 92,93,92

367

93

 

R5=TL/HK(1)*(U5-UL)*2

368

 

 

GO TO 98

369

92

 

R5=TBK

370

 

 

GO TO 98

371

91

 

U5=U(I,J,K-1)

372

 

 

N2=N(I,J,K-1)

373

 

 

HN2=H(I,J,K-1)

374

 

 

IF (HN2-H1(N2)) 95,96,96

375

95

 

TL5=TLT(N2)

376

 

 

GO TO 97

377

96

 

TL5=TLT(N2)

378

 

 

IF (HN2.LE.H2(N2)) TL5=0.5*(TLM(N2)+TL5)

379

97

 

UX=HK(K-1)+HK(K)

380

 

 

R5=2*(HK(K-1)*TL5+HK(K)*TL)/(UX*UX)*(U5-UL)

381

98

 

IF (K.NE.NK) GO TO 99

382

 

 

U6=TBKK

383

 

 

IF (MK(5)-1) 100,101,100

384

101

 

R6=TL/MK(NK)*(U6-UL)*2

385

 

 

GO TO 103

386

100

 

R6=TBKK

387

 

 

GO TO 103

388

99

 

U6=U(I,J,K+1)

389

 

 

N2=N(I,J,K-1)

390

 

 

HK2=H(I,J,K+1)

391

 

 

IF (HK2-H(N2)) 104,105,105

392

104

 

TL6=TLT(N2)

393

 

 

GO TO 106

394

105

 

TL6=TLT(N2)

395

 

 

IF (HN2.LE.H2(N2)) TL6=0.5*(TLM(N2)+TL6)

396

106

 

UX=HK(K+1)*HK(K)

397

 

 

R6=2*(HK(K+1)*TL6+HK(K)*TL)/(UX*UX)*(U6-UL)

398

103

 

CONTINUE

399

 

 

F=0.0

400Р

108

 

HNDV(I,J,K)=HJ+HT*((R1+R2)/HI(I)+(R3+R4)/HJ(J)+

 

 

 

(R5+R6)/HK(K)+1)*NP

402

 

 

HK1=HNDV(I,J,K)

403

 

 

IF (HK1-H1(H1)) 111,112,112

404

112

 

IF (HK1-H2(H2)) 113,113,114

405

113

 

UL1=UKR(N1)

406

 

 

GO TO 115

407

114

 

UL1=(HN1-H2(N1))/CTAL(N1)+UKR(N1)

408

 

 

GO TO 115

409

111

 

IF ((A(N1).EO.0.0).AND.(C(N1).EO.0.0)) GO TO 4

410

 

 

CALL FIND(HN1,N1,UL1)

411

 

 

GO TO 115

412

4

 

UL1=-HN1/H1(N1)*TK+TK

413

115

 

UN(I,J,K)=UL1

414

50

 

CONTINUE

415

32

 

CONTINUE

416

31

 

CONTINUE

417

 

 

DO 42 K=1,NK

418

 

 

DO 42 J=1,NJ

419

 

 

DO 42 I=1,NI

420

 

 

U(I,J,K)=UN(I,J,K)

421

 

 

H(I,J,K)=HNDV(I,J,K)

422

42

 

CONTINUE

423

 

 

IF (T.LT.(TPR-0.01)) GO TO 10

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

<