Инженерные изыскания для строительства. Составление прогноза изменений температурного режима вечномерзлых грунтов численными методами

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),l(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