.
 
(0)  
 >>>     
   [ + ]  
: .,. .
: 1 1997 .
   30319.0-96* " .    .  "
: (30.12.1996)
: 1 1997 .
   30319.1-96* " .    .     ,      "
: .,. .
: 1 1997 .
   30319.3-96* " .    .      "

30319.2-96



,

1 - , , ( ) ;

2 , ( 9-96 12 1996 .)

:

3 , 30 1996 . 723 30319.2-96 1 1997 .

4

5

30319.2-96

Natural gas. Methods of calculation of physical properties.
Definition of compressibility coefficient

1997-07-01

1

: ( ) .

(, , ), .

, , (. 3.2.1).

30319.0.

2

:

30319.0-96 . .

30319.1-96 . . ,

3

3.1

= z/zc, (1)

z zc - .

, , . pc Tc 30319.0.

3.2

3.2.1

1 . [1-12].

0,1 %.

1 -

 

r, /3

r,

d, %

d, %

di, %

NX19 .

 

32£., /3£40

 

<0,70

<3

0,12

-0,02

+0,07

-0,09

3-7

0,18

-0,01

+0,37

-0,10

>7

0,41

0,17

+0,59

-0,08

0,66£r, /3£1,05

 

0,70 0,75

<3

0,13

0,01

+0,14

-0,13

0£x, .%£15

3-7

0,29

0,12

+0,46

-0,15

0£xy, .%£15

>7

0,42

0,27

+0,66

-0,12

250£, £340

 

>0,75

<3

0,20

0,05

+0,41

-0,13

0,1£, £12,0

3-7

0,57

0,24

+1,06

-0,25

 

>7

1,09

0,34

+1,65

-0,40

 

0,74-1,00 ( H2S)

0,1-11

0,15

-0,02

+0,09

-0,10

GERG-91 .

 

20£., /3£48

 

<0,70

<3

0,11

0,01

+0,13

-0,04

3-7

0,15

0,02

+0,51

-0,06

>7

0,20

0,03

+0,63

-0,06

0,66£r, /3£1,05

 

0,70 0,75

<3

0,12

-0,01

+0,08

-0,17

0£x, .%£15

3-7

0,15

-0,02

+0,11

-0,43

0£xy, .%£15

>7

0,19

0,02

+0,16

-0,34

250£, £340

 

>0,75

<3

0,13

0,01

+0,26

-0,12

0,1£, £12,0

3-7

0,15

-0,01

+0,15

-0,30

 

>7

0,19

0,01

+0,65

-0,31

 

0,74-1,00 ( H2S)

0,1-11

2,10

-0,66

+0,06

-3,10

AGA8-92DC

 

20£., /3£48

 

<0,70

<3

0,10

-0,01

+0,03

-0,06

3-7

0,11

-0,01

+0,15

-0,06

>7

0,12

0,02

+0,19

-0,04

0,66£r, /3£1,05

 

0,70 0,75

<3

0,12

-0,01

+0,08

-0,18

0£x, .%£15

3-7

0,15

-0,03

+0,11

-0,43

0£xy, .%£15

>7

0,19

0,01

+0,16

-0,37

250£, £340

 

>0,75

<3

0,12

0,01

+0,25

-0,11

0,1£, £12,0

3-7

0,15

-0,02

+0,24

-0,24

 

>7

0,17

0,01

+0,31

-0,17

 

0,74-1,00 ( H2S)

0,1-11

1,30

-0,38

+0,06

-1,88

20£., /3£48

 

 

<0,70

<3

0,11

-0,04

+0,01

-0,10

 

3-7

0,12

-0,04

+0,05

-0,11

 

>7

0,12

-0,01

+0,06

-0,14

0,66£r, /3£1,05

 

0,70 0,75

<3

0,12

-0,03

+0,08

-0,17

0£x, .%£15

3-7

0,15

-0,02

+0,11

-0,33

0£xy, .%£15

>7

0,18

0,02

+0,13

-0,27

250£, £340

 

>0,75

<3

0,13

-0,01

+0,25

-0,11

0,1£, £12,0

3-7

0,15

-0,01

+0,18

-0,25

 

>7

0,24

-0,01

+0,28

-0,33

 

0,74-1,00 ( H2S)

0,1-11

0,36

0,10

+0,54

-0,24

:

1 NX19 . GERG-91 . (..) (52) 30319.1.

2 AGA8-92DC (r) (16) 30319.1, (..) 7.2 30319.1 ( . (52) 30319.1).

:

1) NX19 .- ;

2) () GERG-91 . [13, 14] AGA8-92DC [15] - ;

3) . d 1 .

NX19 . GERG-91 . , .

AGA8-92DC , ( . %):

65 - 100 £ 15

£ 3,5 £ 1,5

£ 15 £ 15

£ 30 ( ) £ 0,02 ( AGA8-92DC)

£ 1

(12 - 30) (260 - 340) GERG-91 . AGA8-92DC. : GERG-91 . - 3,0 % [14], AGA8-92DC - 0,5 % [15].

.

1 :

1) d -

, (2)

2) di - i-

, (3)

- ;

3) d - 5168 [16]

, (4)

d - ,

, (5)

d - (0,1 %).

( , . 1).

3.2.2 NX19 .

[17] NX19 .

, (6)

, (7)

, (8)

, (9)

, (10)

, (11)

F D :

0 £ £ 2 0 £ D £ 0,3

, (12)

0 £ < 1,3 -0,25 £ D < 0

, (13)

1,3 £ < 2 -0,21 £ D < 0

, (14)

DTa = Ta - 1,09.

a a :

, (15)

, (16)

- , (48) (49) 30319.1, :

, (17)

. (18)

(17), (18) (ry ra).

(1), (6) (18) , (24) 30319.1

( , . 1).

3.2.3 GERG-91 .

, [12], [18], [13, 14]

, (19)

m m - ;

r - , /3.

:

, (20)

, (21)

-

= 1 - - , (22)

, (23)

, (24)

, (25)

, (26)

, (27)

, (28)

, (29)

, (30)

, (31)

, (32)

. (33)

(23), (27)

, (34)

- ,

, (35)

(35) (x) (22), (z) (24) 30319.1,

, (36)

(19) m m (, ) (, )

, (37)

, (38)

, (39)

, (40)

, (41)

0 = b2Cm, (42)

, (43)

(1),

, (44)

z (36).

( , . 1).

3.2.4 AGA8-92DC

ISO/TC 193 SC1 62 [15]

, (45)

n* - ;

r - , /3.

{bn, cn, kn} (45) A.1.

, (12) 30319.1.

, (46)

(53).

:

, (47)

, (48)

N - .

{an, un, gn, qn, fn} {i, i, Gi, Qi, Fi} (47), (48) .1 .2.

{Eij, Gij} {U, G, Km, Q, F} :

, (49)

(i ¹ j)

, (50)

(i ¹ j)

, (51)

, (52)

, (53)

, (54)

, (55)

{Eij*, Gij*, Uij*, Kij*} - , .3.

, , i=j, .

(45) r (, ) (, ).

r (45) :

1)

, (56)

, (57)

2) k-

. (58)

, (59)

z(k-1) - (45) (k-1), .. r(k-1), 1

,(60)

r = 3r(k-1);

4)

, (61)

(61) , , 2) .

(1) .

( , . 1).

3.2.5

- , , ( )

, (62)

ckl - ;

r = r/r - ;

= / - ;

r - , /3;

r - .

, (63)

{akl, bkl} - , .1.

:

-

, (64)

, (65)

(; )

-

, (66)

, (67)

; (68)

(; )

-

, (69)

, (70)

(64) - (70) N - (, , , -, -, , , ).

{ri, rj, j, j}, {i, j,} {Wi, Wj} .2, {xij, lij} - .3 .4.

, , ( 1 % ), :

- ;

- ;

- - -;

- .

, (12) 30319.1.

(62) r (, ) (, ).

r (62) :

1)

, (75)

, (76)

, (77)

(r), () (W) (64), (66) (69);

2) k-

, (78)

, (79)

z(k-1) (62) (k-1), .. r(k-1), a A1

, (80)

4) .

, (81)

(81) , , 2) .

(1) .

( , . 1).

4

, , (), (T), (rc) (i) . .

5168 [16] , ,

, (82)

d - , ;

dqk - ;

, (83)

, (84)

(82) - (84):

qk - k- (. , rc, i,);

`qk - k- (, , ..);

qk qk - k- ;

Nq - .

(83) qk+ qk- .

` ( ) qk.

:

1) NX 19 . GERG-91 . - Nq = 5 , , , ;

2) AGA8-92DC - Nq = 2 + N (N - ) , , .

, (85)

d - , 3.2.1.

NX19 . GERG-91 . d.

, (86)

d, dp, drc, dxa dxy - , , , , , .

T, , rc, xa x , K, (. (34) - (38) (39) - (43) 30319.1):

- NX19 .

, (87)

, (88)

, (89)

, (90)

, (91)

- GERG-91

, (92)

, (93)

, (94)

, (95)

. (96)

( , . 1).

5

, IBM PC/AT/XT, -77. .

, .

()

AGA8-92DC

.1 - AGA8-92DC

bn

cn

kn

un

gn

qn

fn

1

0,153832600

1

0

0

0,0

0

0

0

2

1,341953000

1

0

0

0,5

0

0

0

3

-2,998583000

1

0

0

1,0

0

0

0

4

-0,048312280

1

0

0

3,5

0

0

0

5

0,375796500

1

0

0

-0,5

1

0

0

6

-1,589575000

1

0

0

4,5

1

0

0

7

-0,053588470

1

0

0

0,5

0

1

0

8

2,29129-9

1

1

3

-6,0

0

0

1

9

0,157672400

1

1

2

2,0

0

0

0

10

-0,436386400

1

1

2

3,0

0

0

0

11

-0,044081590

1

1

2

2,0

0

1

0

12

-0,003433888

1

1

4

2,0

0

0

0

13

0,032059050

1

1

4

11,0

0

0

0

14

0,024873550

2

0

0

-0,5

0

0

0

15

0,073322790

2

0

0

0,5

0

0

0

16

-0,001600573

2

1

2

0,0

0

0

0

17

0,642470600

2

1

2

4,0

0

0

0

18

-0,416260100

2

1

2

6,0

0

0

0

19

-0,066899570

2

1

4

21,0

0

0

0

20

0,279179500

2

1

4

23,0

1

0

0

21

-0,696605100

2

1

4

22,0

0

1

0

22

-0,002860589

2

1

4

-1,0

0

0

1

23

-0,008098836

3

0

0

-0,5

0

1

0

24

3,150547000

3

1

1

7,0

1

0

0

25

0,007224479

3

1

1

-1,0

0

0

1

26

-0,705752900

3

1

2

6,0

0

0

0

27

0,534979200

3

1

2

4,0

1

0

0

28

-0,079314910

3

1

3

1,0

1

0

0

29

-1,418465000

3

1

3

9,0

1

0

0

30

-5,99905-17

3

1

4

-13,0

0

0

1

31

0,105840200

3

1

4

21,0

0

0

0

32

0,034317290

3

1

4

8,0

0

1

0

33

-0,007022847

4

0

0

-0,5

0

0

0

34

0,024955870

4

0

0

0,0

0

0

0

35

0,042968180

4

1

2

2,0

0

0

0

36

0,746545300

4

1

2

7,0

0

0

0

37

-0,291961300

4

1

2

9,0

0

1

0

38

7,294616000

4

1

4

22,0

0

0

0

39

-9,936757000

4

1

4

23,0

0

0

0

40

-0,005399808

5

0

0

1,0

0

0

0

41

-0,243256700

5

1

2

9,0

0

0

0

42

0,049870160

5

1

2

3,0

0

1

0

43

0,003733797

5

1

4

8,0

0

0

0

44

1,874951000

5

1

4

23,0

0

1

0

45

0,002168144

6

0

0

1,5

0

0

0

46

-0,658716400

6

1

2

5,0

1

0

0

47

0,000205518

7

0

0

-0,5

0

1

0

48

0,009776195

7

1

2

4,0

0

0

0

49

-0,020487080

8

1

1

7,0

1

0

0

50

0,015573220

8

1

2

3,0

0

0

0

51

0,006862415

8

1

2

0,0

1

0

0

52

-0,001226752

9

1

2

1,0

0

0

0

53

0,002850906

9

1

2

0,0

0

1

0

.2 -

,

, 3/

G

Q

F

16,0430

151,3183

0,4619255

0,0

0,0

0,0

30,0700

244,1667

0,5279209

0,079300

0,0

0,0

44,0970

298,1183

0,5837490

0,141239

0,0

0,0

-

58,1230

337,6389

0,6341423

0,281835

0,0

0,0

-

58,1230

324,0689

0,6406937

0,256692

0,0

0,0

28,0135

99,73778

0,4479153

0,027815

0,0

0,0

44,0100

241,9606

0,4557489

0,189065

0,69

0,0

34,0820

296,3550

0,4618263

0,088500

0,0

0,0

-

72,1500

370,6823

0,6798307

0,366911

0,0

0,0

-

72,1500

365,5999

0,6738577

0,332267

0,0

0,0

-

86,1770

402,8429

0,7139987

0,432254

0,0

0,0

-

100,2040

427,5391

0,7503628

0,512507

0,0

0,0

-

114,2310

450,6472

0,7851933

0,576242

0,0

0,0

4,0026

2,610111

0,3589888

0,0

0,0

0,0

28,0100

105,5348

0,4533894

0,038953

0,0

0,0

31,9988

122,7667

0,4186954

0,021000

0,0

0,0

39,9480

119,6299

0,4216551

0,0

0,0

0,0

18,0153

514,0156

0,3825868

0,332500

0,0

0,0

.3 -

i

J

Eij*

Uij

Kij

Gij*

0,971640

0,886106

1,003630

 

0,960644

0,963827

0,995933

0,807653

0,996050

1,023960

 

 

0,990126

 

 

 

-

1,019530

 

 

 

-

0,995474

1,021280

 

 

-

1,002350

 

 

 

-

1,003050

 

 

 

-

1,012930

 

 

 

-

0,999758

 

 

 

-

0,988563

 

 

 

1,022740

0,835058

0,982361

0,982746

0,970120

0,816431

1,007960

 

0,945939

0,915502

 

 

1,005710

 

 

 

-

0,946914

 

 

 

-

0,973384

0,993556

 

 

-

0,959340

 

 

 

-

0,945520

 

 

 

-

0,937880

 

 

 

-

0,935977

 

 

 

-

0,933269

 

 

 

0,925053

0,969870

1,008510

0,370296

0,960237

 

 

 

1,500000

0,900000

 

 

-

0,906849

 

 

 

-

0,897362

 

 

 

-

0,726255

 

 

 

-

0,859764

 

 

 

-

0,766923

 

 

 

-

0,782718

 

 

 

-

0,805823

 

 

 

1,035020

1,080500

1,000460

 

-

 

1,250000

 

 

-

1,013060

1,250000

 

 

-

 

1,250000

 

 

-

1,005320

1,250000

 

 

-

1,004900

 

 

 

()

.1 -

k

l

akl

bkl

k

l

akl

bkl

1

0

6,087766 × 10-1

-7,187864 × 10-1

8

2

4,015072 ×10-1

-9,576900 × 100

2

0

-4,596885 ×10-1

1,067179 × 101

9

2

-1,016264 × 10-1

2,419650 × 100

3

0

1,149340 × 100

-2,576870 × 101

10

2

-9,129047 × 10-3

2,275036 × 10-1

4

0

-6,075010 × 10-1

1,713395 × 101

1

3

-2,837908 × 100

1,571955 × 101

5

0

-8,940940 × 10-1

1,617303 × 101

2

3

1,534274 × 101

-3,020599 × 102

6

0

1,144404 × 100

-2,438953 × 101

3

3

-2,771885 × 101

6,845968 × 102

7

0

-3,457900 × 10-1

7,156029 × 100

4

3

3,511413 × 101

-8,281484 × 102

8

0

-1,235682 × 10-1

3,350294 × 100

5

3

-2,348500 × 101

5,600892 × 102

9

0

1,098875 × 10-1

-2,806204 × 100

6

3

7,767802 × 100

-1,859581 × 102

10

0

-2,193060 × 10-2

5,728541 × 10-1

7

3

-1,677977 × 100

3,991057 × 101

1

1

-1,832916 × 100

6,057018 × 100

8

3

3,157961 × 10-1

-7,567516 × 100

2

1

4,175759 × 100

-7,947685 × 101

9

3

4,008579 × 10-3

-1,062596 × 10-1

3

1

-9,404549 × 100

2,167887 × 102

1

4

2,606878 × 100

-1,375957 × 101

4

1

1,062713 × 101

-2,447320 × 102

2

4

-1,106722 × 101

2,055410 × 102

5

1

-3,080591 × 100

7,804753 × 101

3

4

1,279987 × 101

-3,252751 × 102

6

1

-2,122525 × 100

4,870601 × 101

4

4

-1,211554 × 101

2,846518 × 102

7

1

1,781466 × 100

-4,192715 × 101

5

4

7,580666 × 100

-1,808168 × 102

8

1

-4,303578 × 10-1

1,000706 × 101

6

4

-1,894086 × 100

4,605637 × 101

9

1

-4,963321 × 10-2

1,237872 × 100

1

5

-1,155750 × 100

6,466081 × 100

10

1

3,474960 × 10-2

-8,610273 × 10-1

2

5

3,601316 × 100

-5,739220 × 101

1

2

1,317145 × 100

-1,295347 × 101

3

5

-7,326041 × 10-1

3,694793 × 101

2

2

-1,073657 × 101

2,208390 × 102

4

5

-1,151685 × 100

2,077675 × 101

3

2

2,395808 × 101

-5,864596 × 102

5

5

5,403439 × 10-1

-1,256783 × 101

4

2

-3,147929 × 101

7,444021 × 102

1

6

9,060572 × 10-2

-9,775244 × 10-1

5

2

1,842846 × 101

-4,470704 × 102

2

6

-5,151915 × 10-1

2,612338 × 100

6

2

-4,092685 × 100

9,965370 × 101

3

6

7,622076 × 10-2

-4,059629 × 10-1

7

2

-1,906595 × 10-1

5,136013 × 100

1

7

4,507142 × 10-2

-2,298833 × 10-1

.2 - ,

i

rci, /3

Wi

pki,

rki, /3

Tki, K

zki,

4

16,043

4,5988

163,03

190,67

0,2862

0,6682

0,0006467

26

30,070

4,88

205,53

305,57

0,2822

1,2601

0,1103

38

44,097

4,25

218,54

369,96

0,2787

1,8641

0,1764

-

-410

58,123

3,784

226,69

425,40

0,2761

2,4956

0,2213

-

-410

58,123

3,648

225,64

407,96

0,2769

2,488

0,2162

N2

28,0135

3,390

315,36

125,65

0,2850

1,16490

0,04185

2

44,010

7,386

466,74

304,11

0,2744

1,8393

0,2203

H2S

34,082

8,940

349,37

373,18

0,2810

1,4311

0,042686

1 (rki), (Tki) (Wi) .

2 rci - i-

.3 - xij

j

i

4

C2H6

38

-C4H10

-410

N2

CO2

H2S

4

0,0

0,036

0,076

0,121

0,129

0,060

0,074

0,089

C2H6

-

0,0

0,0

0,0

0,0

0,106

0,093

0,079

38

-

-

0,0

0,0

0,0

0,0

0,0

0,0

-C410

-

-

-

0,0

0,0

0,0

0,0

0,0

-410

-

-

-

-

0,0

0,0

0,0

0,0

N2

-

-

-

-

-

0,0

0,022

0,211

CO2

-

-

-

-

-

-

0,0

0,089

H2S

-

-

-

-

-

-

-

0,0

.4 - lij

j

i

4

26

38

-410

-C4H10

N2

2

H2S

4

0,0

-0,074

-0,146

-0,258

-0,222

-0,023

-0,086

0,0

26

-

0,0

0,0

0,0

0,0

0,0

0,0

0,0

38

-

-

0,0

0,0

0,0

0,0

0,0

0,0

-C4H10

-

-

-

0,0

0,0

0,0

0,0

0,0

-410

-

-

-

-

0,0

0,0

0,0

0,0

N2

-

-

-

-

-

0,0

-0,064

0,0

2

-

-

-

-

-

-

0,0

-0,062

H2S

-

-

-

-

-

-

-

0,0

()

C **********************************************************

C * *

* *

* ( ) *

* *

C **********************************************************

IMPLICIT REAL*8(A-H,O-Z)

CHARACTER*26 AR(25)

DIMENSION PI(100),TI(100),ZP(100,100)

COMMON/P/P/T/T/RON/RON/YI/YC(25)/Z/Z/NPR/NPR

DATA AR/ (4), (26), (38),

* - (-410), - (-410), (N2),

* (O2), (H2S),

* (22), (24), (36),

* - (-512), - (-C5H12),

* - (-512), - (-614),

* (66), - (-716), (78),

* - (-818), - (-920),

* - (-1022), (), (2),

* (), (2)/

200 WRITE(*,100)

CALL VAR(NVAR)

IF(NVAR.EQ.5) GO TO 134

WRITE(*,l00)

100 FORMAT(25(/))

WRITE(*,1)

1 FORMAT( ./)

IF(NVAR.LE.2) THEN

WRITE(*,(A\))

* 293.15 101.325 , /.

READ(*,*)RON

WRITE(*,53)

53 FORMAT( 0, ,

* /

* 1, ,

* \)

READ(*,*)NPR

IF(NPR.EQ.0) WRITE(*,3)

3 FORMAT ( , . %)

IF(NPR.EQ.l) WRITE(*,33)

33 FORMAT( , . %)

WRITE(*,(A\)) (N2)

READ(*,*)YA

YA = YA/100.

WRITE(*,(A\)) (02)

READ(*,*)YY

YY = YY/100.

ELSE

WRITE(*,35)

35 FORMAT( 0, /

* 1, \)

READ(*,*)NPR

IF(NPR.EQ.0) WRITE(*,3)

IF(NPR.EQ.l) WRITE(*,33)

DO 5 I=1,25

WRITE(*,(A\)) AR(I)

READ(*,*)YC(I)

5 YC(I) = YC(I)/100.

ENDIF

WRITE(*,(A\))

* :

READ(*,*)NP

WRITE(*,(A\))

* :

READ(*,*)NT

WRITE(*,(A\))

* :

READ(*,*)(PI(I),I=1,NP)

WRITE(*,(A\))

* :

READ(*,*)(TI(I),I=1,NT)

WRITE(*,(A\))

* .

P=.101325D0

T=293.15D0

ICALC=1

GO TO (10,20,30,40) NVAR

10 CALL NX19(YA,YY)

ZN=Z

GO TO 50

20 CALL GERG2(ICALC,YA,YY)

ZN=Z

GO TO 50

30 CALL AGA8DC(ICALC)

ZN=Z

GO TO 50

40 CALL VNIC(ICALC)

ZN=Z

50 CONTINUE

IF(Z.EQ.0D0) THEN

CALL RANGE(NRANGE)

IF(NRANGE) 134,134,200

ENDIF

ICALC=2

NTS=0

DO 7 I=1,NP

P=PI(I)

D07 J=1,NT

T=TI(J)

IF(NVAR.EQ.l) CALL NX19(YA,YY)

IF(NVAR.EQ.2) CALL GERG2(ICALC,YA,YY)

IF(NVAR.EQ.3) CALL AGA8DC(ICALC)

IF(NVAR.EQ.4) CALL VNIC(ICALC)

IF(Z.NE.0D0) NTS=NTS+1

ZP(I,J)=Z/ZN

7 CONTINUE

IF(NTS.EQ.0) THEN

CALL RANGE(NRANGE)

IF (NRANGE) 134,134,200

ELSE

I=1

9 I=0

DO 11 J=1,NT

IF(ZP(I,J).EQ.0D0)

IC=IC+1

11 CONTINUE

IF(IC.EQ.NT) THEN

IF(I.NE.NP) THEN

DO 13 J=I,NP-1

PI(J)=PI(J+1)

DO 13 K=1,NT

13 ZP(J,K)=ZP(J+1,K)

ENDIF

NP=NP-1

ELSE

I=I+1

ENDIF

IF(I.LE.NP) GO TO 9

J=l

15 JS=0

DO 17 I=1,NP

IF(ZP(I,J).EQ.0D0) JS=JS+1

17 CONTINUE

IF(JS.EQ.NP) THEN

IF(J.NE.NT) THEN

DO 19 I=J,NT-1

I(I)=I(I+1)

DO 19 K=1,NP

19 ZP(K,I)=ZP(K,I+1)

ENDIF

NT=NT-1

ELSE

J=J+1

ENDIF

IF(J.LE.NT) GO TO 15

CALL TABL(YA,YY,PI,TI,ZP,NP,NT,NVAR,AR)

ENDIF

GO TO 200

134 STOP

END

SUBROUTINE VAR(NVAR)

WRITE(*,1)

1 FORMAT(//

*10X, //

*10, ---------------- ----------------- /

*10, /

*10, 1. NX 19 /

*10, /

*10, 2. GERG-91 /

*10, /

*10, 3. AGA8-92DC /

*10, /

*10, 4. /

*10, /

*10,---------------------------------------------------/)

WRITE(*,5)

5 FORMAT(/,3X,

* 5 ђ,

*\)

READ(*,*)NVAR

RETURN

END

SUBROUTINE RANGE(NRANGE)

IMPLICIT REAL*8(A-H,-Z)

COMMON/Z/Z

WRITE(*,1)

1 FORMAT(//

* /

* ? 0 - , 1 - \)

READ(*,*)NRANGE

RETURN

END

SUBROUTINE TABL(YA,YY,PI,TI,ZP,NP,NT,NVAR,AR)

IMPLICIT REAL*8(A-H,-Z)

CHARACTER*26 AR(25), FNAME

CHARACTER METH(4)*31,A*6,LIN1(5)*9,LIN2(5)*9,LIN3(6)*9,LIN4*9,

*AT(06)*28

CHARACTER*70 F,FZ(11,2)

DIMENSION PI(100),TI(100),ZP(100,100),ZPP(6)

COMMON/RON/RON/YI/YC(25)/NPR/NPR

DATA METH/

*( NX19),

*( GERG-91),

*( AGA8-92DC),

*( )/

DATA LIN1/5*------/,LIN2/5*------/,LIN3/6*------/,

*LIN4/------/,A/ - /

DATA AT/

* T, K, T, K, T, K, T,K,

* T, K, T, K/

DATA FZ/

*(3X,F5.2,2X,6(3X,F6.4)),(3X,F5.2,5X,A6,5(3X,F6.4)),

*(3X,F5.2,2X,2(3X,A6),4(3X,F6.4)),(3X,F5.2,2X,3(3X,A6),

*3(3X,F6.4)),

*(3X,F5.2,2X,4(3X,A6),2(3X,F6.4)),(3X,F5.2,2X,5(3X,A6),

*3X,F6.4),

*(3X,F5.2,2X,5(3X,F6.4),3X,A6),(3X,F5.2,2X,4(3X,F6.4),

*2(3X,A6)),

*(3X,F5.2,2X,3(3X,F6.4),3(3X,A6)),(3X,F5.2,2X,2(3X,F6.4),

*4(3X,A6)),

*(3X,F5.2,5X,F6.4,5(3X,A6)),(3X,F9.6,1X,F6.4,5(3X,F6.4)),

*(X,F9.6,lX,A6,5(3X,F6.4)),(3X,F9.,lX,A6,3X,A6,4(3X,F6.4)),

*(3X,F9.6,1X,A6,2(3X,A6),3(3X,F6.4)),(3X,F9.6,1X,A6,3(3X,A6),

*2(3X,F6.4)),

*(3X,F9.6,1X,A6,4(3X,A6),3X,F6.4),(3X,F9.6,1X,F6.4,4(3X,F6.4),

*3X,A6),

*(3X,F9.6,1X,F6.4,3(3X,F6.4),2(3X,A6)),(3X,F9.6,1X,F6.4),

*2(3X,F6.4),3(3X,A6)),

*(3X,F9.6,1X,F6.4,3X,F6.4,4(3X,A6)),(3X,F9.6,1X,F6.4,5(3X,A6))/

22 WRITE(*,44)

44 FORMAT(// ?,)

WRITE(*,(A\))

* 0 - , 1 - , 2 -

READ(*,*)NYST

IF(NYST.EQ.0) OPEN(1,FILE=CON)

IF(NYST.EQ.l) OPEN(1,FILE=PRN)

IF(NYST.EQ.2) WRITE(*,(A\))

IF(NYST.EQ.2) READ(*,(A))FNAME

IF(NYST.EQ.2) OPEN(1,FILE=FNAME)

IF(NYST.EQ.0) WRITE(*,100)

100 FORMAT(25(/))

IF(NYST.EQ.l) PAUSE

* , <>

WRITE(1,88)METH(NVAR)

88 FORMAT(

*13X, ./

*18,31/)

NW=3

IF(NVAR.LE.2) THEN

WRITE(1,1)RON

1 FORMAT( 293.15 101.325 ,F6.4, /.)

NW=NW+1

IF(YA.NE.0D0.OR.YY.NE.0D0) THEN

IF(NPR.EQ.0) WRITE(1,3)

3 FORMAT( . %)

IF(NPR.EQ.l) WRITE(1,33)

33 FORMAT( .%)

NW=NW+1

IF(YA.NE.0D0) THEN

WRITE(1,5)AR(6),YA* 100.

5 FORMAT(2(A26,F7.4))

NW=NW+1

ENDIF

IF(YY.NE.0D0) THEN

WRITE(1,5)AR(7),YY*100.

NW=NW+1

ENDIF

ENDIF

ELSE

IF(NPR.EQ.0) WRITE(1,3)

IF(NPR.EQ.l) WRITE(1,33)

NW=NW+1

I=1

9 J=I+1

13 CONTINUE

IF(YC(J).NE.0D0) THEN

WRITE(1,5)AR(I),YC(I)*100.,AR(J),YC(J)*100.

NW=NW+1

DO 11 I=J+1,25

IF(YC(I).NE.0D0.AND.I.NE.25) GO TO 9

IF(YC(I).NE.0D0.AND.I.EQ.25) THEN

WRITE(1,5)AR(I),YC(I)*100.

nw=nw+1

GO TO 99

ENDIF

11 CONTINUE

ELSE

J=J+1

IF(J.LE.25) THEN

GO TO 13

ELSE

WRITE(1,5)AR(I),YC(I)*100.

NW=NW+1

ENDIF

ENDIF

ENDIF

99 CONTINUE

IF(NW.GT.12.AND.NYST.EQ.0) THEN

WRITE(*,7)

7 FORMAT(/)

PAUSE <>

WRITE(*,100)

NW=0

ENDIF

DO 15 I=1,NT,6

IF(NW.GT.12.AND.NYST.EQ.0) THEN

WRITE(*,7)

PAUSE <>

WRITE(*,100)

NW=0

ENDIF

IF(NW.GT.46.AND.NYST.NE.O) THEN

WRITE(1,7)

WRITE(*,7)

IF(NYST.EQ.l)

PAUSE

* <>

NW=0

ENDIF

IF(I+5.LE.NT) THEN

NL=6

ELSE

NL=NT-I+1

ENDIF

WRITE(1,7)

IF(NL.GT.1) WRITE(1,17)LIN2(1),(LIN1(K),K=1,NL-1)

IF(NL.EQ.l) WRITE(1,17)LIN2(1)

17 FORMAT( ------,6A9)

WRITE(1,19)AT(NL)

19 FORMAT( ------,A28)

IF(NL.GT.1)WRITE(1,21)LIN4,(LIN2(K),K=1,NL-1)

IF(NL.EQ.l) WRITE(1,21)LIN4

21 FORMAT( p, ,69)

WRITE(1,23)(TI(K),K=I,I+NL-1)

23 FORMAT(10X,6(:,|,F6.2))

WRITE(1,17)(LIN3(K),K=1,NL)

NW=NW+6

DO 25 J=1,NP

JP=1

IF(PI(J).EQ.0.101325D0) JP=2

NL1=0

NLN=0

DO 27K=I,I+NL-1

NL1=NL1+1

IF(ZP(J,K).EQ.0D0) THEN

ZPP(NL1)=A

NLN=NLN+1

ELSE

ZPP(NL1)=ZP(J,K)

ENDIF

27 CONTINUE

IF(NLN.EQ.NL) GO TO 133

IF(NLN.EQ.0) THEN

F=FZ(1,JP)

ELSE

IF(ZP(J,I).EQ.0D0) F=FZ(NLN+1,JP)

IF(ZP(J,I+NL-1).EQ.0D0) F=FZ(NLN+12-NL,JP)

ENDIF

IF(NLI.EQ.1)WRITE(1,F)PI(J),ZPP(1)

IF(NL1.EQ.2)WRITE(1,F)PI(J),ZPP(1),ZPP(2)

IF(NL1.EQ.3)WRITE(1,F)PI(J),ZPP(1),ZPP(2),ZPP(3)

IF(NL1.EQ.4)WRITE(1,F)PI(J),ZPP(1),ZPP(2),ZPP(3),ZPP(4)

IF(NL1.EQ.5)

*WRITE(1,F)PI(J),ZPP(1),ZPP(2),ZPP(3),ZPP(4),ZPP(5)

IF(NL1.EQ.6)

*WRITE(1,F)PI(J),ZPP(1),ZPP(2),ZPP(3),ZPP(4),ZPP(5),ZPP(6)

NW=NW+1

133 CONTINUE

IF(NW.EQ.20.AND.NYST.EQ.0) THEN

IF(J.EQ.NP.AND.I+NL-1.EQ.NT) GO TO 29

WRITE(*,7)

PAUSE <>

WRITE(*,100)

NW=0

WRITE(1,7)

IF(NL.GT.1)WRITE(1,17)LIN2(1),(LIN1(K),K=1,NL-1)

IF(NL.EQ.l) WRITE(1,17)LIN2(1)

WRITE(1,19)AT(NL)

IF(NL.GT.1)WRITE(1,21)LIN4,(LIN2(K),K=1,NL-1)

IF(NL.EQ.l) WRITE(1,21)LIN4

WRITE(1,23)(TI(K),K=I,I+NL-1)

WRITE(1,17)(LIN3(K),K=1,NL)

NW=NW+6

ENDIF

IF(NW.EQ.54.AND.NYST.NE.0) THEN

IF(J.EQ.NP.AND.I+NL-1.EQ.NT) GO TO 29

WRITE(1,7)

WRITE(*,7)

IF(NYST.EQ.l) PAUSE

* <>

NW=0

IF(NL.GT.1) WRITE(1,17)LIN2(1),(LIN1(K),K=1,NL-1)

IF(NL.EQ.l) WRITE(1,17)LIN2(1)

WRITE(1,19)AT(NL)

IF(NL.GT.1) WRITE(1,21)LIN4,(LIN2(K),K=1,NL-1)

IF(NL.EQ.l) WRITE(1,21)LIN4

WRITE(1,23)(TI(K),K=I,I+NL-1)

WRITE(1,17)(LIN3(K),K=1,NL)

NW=NW+6

ENDIF

25 CONTINUE

15 CONTINUE

29 CLOSE(l)

WRITE(*,7)

PAUSE , <>

WRITE(*,66)

66 FORMAT(/ ?,

*, 0 - , 1 - \)

READ(*,*)NBOLB

IF(NBOLB.EQ.l) GO TO 22

RETURN

END

C **********************************************************

* *

* *

* NX19. *

C **********************************************************

SUBROUTINE NX19(YA,YY)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/NCONT/NCONT/YA/Y(2)/RON/RON

Y(1)=YA

Y(2)=YY

CALL PTCONT

IF(NCONT.EQ.l) GO TO 134

CALL EA

CALL PHASEA

134 RETURN

END

SUBROUTINE PTCONT

IMPLICIT REAL*8(A-H,O-Z)

COMMON/NCONT/NCONT/Z/Z/P/P/T/T/YA/Y(2)/RON/RON

NCONT=0

IF(RON.LT.0.66D0.OR.RON.GT.1D0) NCONT=1

IF(Y(1).GT.0.2D0.OR.Y(2).GT.0.15D0) NCONT=l

IF(P.LE.0.D0.OR.T.LE.0.D0) NCONT=1

IF(T.LT.250.D0.OR.T.GT.340.D0) NCONT=1

IF(P.GT.12.D0) NCONT=1

IF(NCONT.EQ.1) Z=0D0

RETURN

END

SUBROUTINE EA

IMPLICIT REAL*8(A-H,O-Z)

COMMON/T/T/YA/Y(2)/RON/RON/P/P/PT/PA,TA/BI/B1,B2/T0/T0

PCM=2.9585*(1.608D0-0.05994*RON+Y(2)-.392*Y(1))

TCM=88.25*(0.9915D0+1.759*RON-Y(2)-1.681*Y(1))

PA=0.6714*P/PCM+0.0147

TA=0.71892*T/TCM+0.0007

DTA=TA-1.09D0

F=0D0

IF(PA.GE.0D0.AND.PA.LT.2D0.AND.DTA.GE.0D0.AND.DTA.LT.0.3D0)

F=75D-5*PA**2.3/DEXP(20.*DTA)+

*11D-4*DTA**0.5*(PA*(2.17D0-PA+1.4*DTA**0.5))**2

IF(PA.GE.0D0.AND.PA.LT.1.3D0.AND.DTA.GE.-0.25D0.AND.DTA.LT.0D0)

*F=75D-5*PA**2.3*(2D0-DEXP(20.*DTA))+

*1.317*PA*(1.69D0-PA**2)*DTA**4

IF(PA.GE.1.3D0.AND.PA.LT.2D0.AND.DTA.GE.-0.21D0.AND.DTA.LT.0D0)

*F=75D-5*PA**2.3*(2D0-DEXP(20.*DTA))+

*0.455*(1.3D0-PA)*(1.69*2.D0**1.25-PA**2)*(DTA*(0.03249D0+

*18.028*DTA**2)+DTA**2*(2.0167D0+DTA**2*(42.844D0+200.*DTA**2)))

T1=:TA**5/(TA**2*(6.60756*TA-4.42646D0)+3.22706D0)

T0=(TA**2*(1.77218D0-0.8879*TA)+0.305131D0)*T1/TA**4

B1=2.*T1/3.-TO**2

B0=T0*(T1-T0**2)+0.1*T1*PA*(F-1D0)

B2=(B0+(B0**2+B1**3)**0.5)**(1D0/3D0)

RETURN

END

SUBROUTINE PHASEA

IMPLICIT REAL*8(A-H,O-Z)

COMMON/Z/Z/PT/PA,TA/BI/B1,B2/T0/T0

Z=(1D0+0.00132/TA**3.25)**2*0.1*PA/(B1/B2-B2+T0)

RETURN

END

C *************************************************************

* *

* *

* GERG-91. *

* *

C *************************************************************

$NOTRUNCATE

SUBROUTINE GERG2(ICALC,YA,YY)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/T/T1/P/PRESS/RON/RON/Z/Z

COMMON/XBLOK/X1,X2,X3,X11,X12,X13,X22,X23,X33

COMMON/MBLOK/GM2,GM3,FA,FB,TO,R

DATABMO/.0838137D0/,BM1/-.00851644D0/,WD0/134.2153D0/,

*WD1/1067.943D0/

Z=-1D0

IF(ICALC.EQ.2) GO TO 3

X2=YA

X3=YY

IF(RON.LT.0.66D0.OR.RON.GT.1D0)Z=0D0

IF(X2.LT.0D0.OR.X2.GT.0.2D0)Z=0D0

IF(X3.LT.0D0.OR.X3.GT.0.15D0) Z=0D0

IF(Z.EQ.0D0) GO TO 133

X1=1D0-X2-X3

11=1*1

X12=X1*X2

X13=X1*X3

X22=X2*X2

X23=X2*X3

X33=X3*X3

Z=1D0-(.0741*RON-.006D0-.063*YA-.0575*YY)**2

BMNG=24.05525*Z*RON

Y1=1D0-YA-YY

BMY=(BMNG-28.0135*YA-44.01*YY)/Y1

()

H=47.479*BMY+128.64D0

RETURN

3 =1

=1-0

P=PRESS

IF(PRESS.LE.0D0.OR.PRESS.GT.12D0)Z=0D0

IF(T1.LT.250D0.OR.T1.GT.340D0)Z=0D0

IF(Z.EQ.0D0) GO TO 133

CALL B11BER(T,H,B11)

CALL BBER(T,B11,B,Z)

IF(Z.EQ.0D0) GO TO 133

CALL CBER(T,H,C,Z)

IF(Z.EQ.0D0) GO TO 133

CALL ITER2(P,T,B,C,Z)

133 RETURN

END

SUBROUTINE B11BER(T,H,B11)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/BBLOK/BR11H0(3),BR11H1(3),BR11H2(3),BR22(3),BR23(3),BR33(3)

T2=T*T

B11=BR11H0(1)+BR11H0(2)*T+BR11H0(3)*T2+

*(BR11H1(1)+BR11H1(2)*T+BR11H1(3)*T2)*H+

*(BR11H2(1)+BR11H2(2)*T+BR11H2(3)*T2)*H*H

END

SUBROUTINE BBER(T,B11,BEFF,Z)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/BBLOK/BR11H0(3),BR11H1(3),BR11H2(3),BR22(3),BR23(3),BR33(3)

COMMON/ZETA/Z12,Z13,Y12,Y13,Y123

COMMON/XBLOK/X1,X2,X3,X11,X12,X13,X22,X23,X33

T2=T*T

B22=BR22(1)+BR22(2)*T+BR22(3)*T2

B23=BR23(1)+BR23(2)*T+BR23(3)*T2

B33=BR33(1)+BR33(2)*T+BR33(3)*T2

BA13=B11*B33

IF(BA13.LT.0D0) THEN

Z=0D0

RETURN

ENDIF

ZZZ=Z12+(320D0-T)**2*1.875D-5

BEFF=X11*B11+X12*ZZZ*(B11+B22)+2.*X13*Z13*DSQRT(BA13)+

*X22*B22+2.*X23*B23+X33*B33

END

SUBROUTINE CBER(T,H,CEFF,Z)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/CBLOK/CR111H0(3),CR111H1(3),CR111H2(3),CR222(3),CR223(3),

*CR233(3),CR333(3)

COMMON/ZETA/Z12,Z13,Y12,Y13,Y123

COMMON/XBLOK/X1,X2,X3,X11,X12,X13,X22,X23,X33

T2=T*T

C111=CR111 H0(1)+CR111H0(2)*T+CR111H0(3)*T2+

*(CR111H1(1)+CR111H1(2)*T+CR111H1(3)*T2)*H+

*(CR111H2(1)+CR111H2(2)*T+CR111H2(3)*T2)H*H

C222=CR222(1)+CR222(2)*T+CR222(3)*T2

C223=CR223(1)+CR223(2)*T+CR223(3)*T2

C233=CR233(1)+CR233(2)*T+CR233(3)*T2

C333=CR333(1)+CR333(2)*T+CR333(3)*T2

CA112=C111*C111*C222

CA113=C111*C111*C333

CA122=C111*C222*C222

CA123=C111*C222*C333

CA133=C111<C333*C333

IF(CA112.LT.0D0.OR.CA113.LT.0D0.OR.CA122.LT.0DO0.OR.

*CA123.LT.0D0.OR.CA133.LT.0D0)THEN

Z=0D0

RETURN

ENDIF

D3REP=1D0/3D0

CEFF=X1*X11*C111+3D0*X11*X2*(CA112)**D3REP*(Y12+(T-270D0)*.0013D0)

*+3.*X11*X3*(CA113)**D3REP*Y13+

*3.*X1*X22*(CA122)**D3REP*(Y12+(T-270D0)*.0013D0)+

*6.*X1*X2*X3*(CA123)**D3REP*Y123+3.*X1*X33*(CA133)**D3REP*Y13+

*X22*X2*C222+3.*X22*X3*C223+3.*X2*X33*C233+X3*X33*C333

END

,

SUBROUTINE ITER2(P,T,Bm,Cm,Z)

IMPLICIT REAL*8(A-H,O-Z)

B1=1D3*P/2.7715/T

B0=Bl*Bm

C0=Bl**2*Cm

A1=1D0+B0

A0=1D0+1.5*(B0+C0)

A01=A0**2-A1**3

IF(A01.LE.0D0) THEN

Z=0D0

RETURN

ENDIF

A=A0-A01**0.5

A2=DABS(A)**(1D0/3D0)

IF(A-LT.0D0) A2=-A2

Z=(1D0+A2+A1/A2)/3.

END

BLOCK DATA BDGRG2

IMPLICIT REAL*8(A-H,O-Z)

COMMON/BBLOK/BR11H0(3),BR11H1(3),BR11H2(3),BR22(3),BR23(3),

*BR33(3)/CBLOK/CR111H0(3),CR111H1(3),CR111H2(3),CR222(3),

*CR223(3),CR233(3),CR333(3)

COMMON/ZETA/Z12,Z13,Y12,Y13,Y123

COMMON/MBLOK/GM2,GM3,FA,FB,TO,R

DATA BR11H0/-.425468D0,.2865D-2,-.462073D-5/,

* BR11H1/.877118D-3,-.556281D-5,.881514D-8/,

* BR11H2/-.824747D-6,.431436D-8,-.608319D-11/,

* BR22/-.1446D0,.74091D-3,-.91195D-6/,

* BR23/-.339693D0,.161176D-2,-.204429D-5/,

* BR33/-.86834D0,.40376D-2,-.51657D-5/

DATA CR111H0/-.302488D0,.195861D-2,-.316302D-5/,

* CR111 H1/.646422D-3,-.422876D-5,.688157D-8/,

* CR111H2/-.332805D-6,.22316D-8,-.367713D-11/,

* CR222/.78498D-2,-.39895D-4,.61187D-7/,

* CR223/.552066D-2,-.168609D-4,.157169D-7/,

* CR233/.358783D-2,.806674D-5,-.325798D-7/,

* CR333/.20513D-2,.34888D-4,-.83703D-7/

DATA Z12/.72D0/,Z13/-.865D0/,Y12/.92D0/,Y13/.92D0/,Y123/1.1D0/

DATA GM2/28.0135D0/,GM3/44.01D0/,

* FA/22.414097D0/,FB/22.710811D0/,

* TO/273.15D0/,R/.0831451D0/

END 46

C **********************************************************

* *

* *

* AGA8-92DC. *

C * *

C **********************************************************

SUBROUTINE AGA8DC(ICALC)

IMPLICIT REAL*8(A-H,O-Z)

REAL*8 KI,KIJ,KD

COMMON/RM/RM/Y1/Y(19)/NC1/NC/NI1/NI(19)/EFI/EI(19),KI(19),

*GI(19),QI(19),FI(19)

*/INTER1/EIJ(19,19),UIJ(19,19),KIJ(19,19),GIJ(19,19)

*/EFD/ED(19),KD(19),GD(19),QD(19),FD(19)/Z/Z

RM=8.31448D0

IF(ICALC.NE.l) GO TO 3

CALL COMPO1

IF(Z.EQ.0D0) GO TO 133

CALL PARIN1

DO 75 I=1,NC

EI(I)=ED(NI(I))

KI(I)=KD(NI(I))

GI(I)=GD(NI(I))

QI(I)=QD(NI(I))

FI(I)=FD(NI(I))

DO 123 J=1,NC

IF(I.GE.J) GO TO 123

EIJ(I,J)=EIJ(NI(I),NI(J))

UIJ(I,J)=UIJ(NI(I),NI(J))

KIJ(I,J)=KIJ(NI(I),NI(J))

GIJ(I,J)=GIJ(NI(I),NI(J))

123 CONTINUE

75 CONTINUE

CALL PARMI1

3 CALL PHASE1

133 RETURN

END

SUBROUTINE COMPO1

IMPLICIT REAL*8(A-h,O-Z)

DIMENSION ZNI(25),YI(25)

COMMON/YI/Y(19)/YI/YC(25)/NC1/NC/NT1/NI(19)/NPR/NPR

DATA ZNI/.9981D0,.992D0,.9834D0,.9682D0,.971D0,.9997D0,.9947D0,

*.99D0,.993D0,.994D0,985D0,.945D0,.953D0,1D0,.919D0,

*.936D0,.876D0,.892D0,3*1D0,1.0005D0,1.0006D0,.9996D0,.9993D0/

DO 100 I=1,25

100 YI(I)=YC(I)

YI(13)=YI(13)+YI(14)

YI(14)=0D0

IF(NPR.EQ.0D0) GO TO 5

YI(17)=YI(17)+YI(19)+YI(20)+YI(21)

YI(19)=0D0

YI(20)=0D0

YI(21)=0D0

SUM=0D0

DO 7 I=1,25

7 SUM=SUM+YI(I)/ZNI(I)

DO 9 I=1,25

9 YI(I)=YI(I)/ZNI(I)/SUM

5 YI(2)=YI(2)+YI(9)+YI(10)

YI(9)=0D0

YI(10)=O0D0

YI(3)=YI(3)+YI(11)

YI(11)=0D0

YI(15)=YI(15)+YI(16)

YI(16)=0D0

YI(17)=YI(17)+YI(18)

YI(18)=0D0

NC=0

=0

YSUM=0D0

DO 11 1=1,25

IF((I.GE.9.AND.I.LE.11).OR.I.EQ.14.0R.I.EQ.16.0R.I.EQ.18)

*=+1

IF(YI(I).EQ.0D0) GO TO 11

NC=NC+1

NI(NC)=I-C

Y(NC)=YI(I)

YSUM=YSUM+Y(NC)

11 CONTINUE

CALL MOLDOl(YI)

DO 13 I=1,NC

13 Y(I)=Y(I)/YSUM

RETURN

END

SUBROUTINE MOLDO1(YI)

IMPLICIT REAL*8(A-H,O-Z)

DIMENSION YI(25)

COMMON/Z/Z

Z=-1D0

YS=0D0

DO 1 I=9,25

1 YS=YS+YI(I)

1F(YI(1).LT.0.65D0.OR.YI(2).GT.0.15D0.OR.YI(3).GT.0.035D0.OR.

*YI(4).GT.0.015D0.OR.YI(5).GT.0.015D0.OR.YS.GT.0.01D0) Z=0D0

IF(YI(6).GT.0.2D0.OR.YI(7).GT.0.15D0.OR.Y1(8).GT.5D-5) Z=O0D0

RETURN

END

SUBROUTINE PARIN1

IMPLICIT REAL*8(A-H,O-Z)

REAL*8 KIJ

COMMON/INTER1/EIJ(19,19),UIJ(19,19),KIJ(19,19),GIJ(19,19)

DO 1 I=1,19

DO 1 J=1,19

EIJ(I,J)=1D0

UIJ(I,J)=1D0

KIJ(I,J)=1D0

1 GIJ(I,J)=1D0

EIJ(1,6)=0.97164D0

UIJ(1,6)=0.886106D0

KIJ(1,6)=1.00363D0

EIJ(1,7)=0.960644D0

UIJ(1,7)=0.963827D0

KIJ(1,7)=0.995933D0

GIJ(1,7)=0.807653D0

EIJ(1,3)=0.99605D0

UIJ(1,3)=1.02396D0

EU(1,17)=1.17052D0

UIJ(1,17)=1.15639D0

KIJ(1,17)=1.02326D0

GIJ(1,17)=1.95731D0

EIJ(1,18)=0.990126D0

EIJ(1,5)=1.01953D0

EIJ(1,4)=0.995474D0

UIJ(1,4)=1.02128D0

EIJ(1,10)=1.00235D0

EIJ(1,9)=1.00305D0

EIJ(1,11)=1.01293D0

EIJ(1,12)=0.999758D0

EIJ(1,13)=0.988563D0

EIJ(6,7)=1.02274D0

UIJ(6,7)=0.835058D0

KIJ(6,7)=0.982361D0

GIJ(6,7)=0.982746D0

EIJ(2,6)=0.97012D0

UIJ(2,6)=0.816431D0

KIJ(2,6)=1.00796D0

EIJ(3,6)=0.945939D0

UIJ(3,6)=0.915502D0

EIJ(6,17)=1.08632D0

UIJ(6,17)=0.408838D0

KIJ(6,17)=1.03227D0

EIJ(6,18)=1.00571D0

EIJ(5,6)=0.946914D0

EIJ(4,6)=0.973384D0

UIJ(4,6)=0.993556D0

EIJ(6,10)=0.95934D0

EIJ(6,9)=0.94552D0

EIJ(6,11)=0.93788D0

EIJ(6,12)=0.935977D0

EIJ(6,13)=0.933269D0

EIJ(2,7)=0.925053D0

UIJ(2,7)=0.96987D0

KIJ(2,7)=1.00851D0

GIJ(2,7)=0.370296D0

EIJ(3,7)=0.960237D0

EIJ(7,17)=1.28179D0

EIJ(7,18)=1.5D0

UIJ(7,18)=0.9D0

EIJ(5,7)=0.906849D0

EIJ(4,7)=0.897362D0

EIJ(7,10)=0.726255D0

EIJ(7,9)=0.859764D0

EIJ(7,11)=0.766923D0

EIJ(7,12)=0.782718D0

EIJ(7,13)=0.805823D0

EIJ(2,3)=1.03502D0

UIJ(2,3)=1.0805D0

KIJ(2,3)=1.00046D0

EIJ(2,17)=1.16446D0

UIJ(2,17)=1.61666D0

KIJ(2,17)=1.02034D0

UIJ(2,5)=1.25D0

EIJ(2,4)=1.01306D0

UIJ(2,4)=1.25D0

UIJ(2,10)=1.25D0

EIJ(2,9)=1.00532D0

UIJ(2,9)=1.25D0

EIJ(3,17)=1.034787D0

EIJ(3,4)=1.0049D0

EIJ(17,18)=1.1D0

EIJ(5,17)=1.3D0

EIJ(4,17)=1.3D0

RETURN

END

SUBROUTINE PARMI1

IMPLICIT REAL*8(A-H,O-Z)

REAL*8 KI,KIJ,KM

INTEGER GN,QN,FN

DIMENSION EIJM(19,19),GIJM(19,19)

COMMON/Y1/Y(19)/NC1/NC/EFI/EI(19),KI(19),GI(19),QI(19),FI(19)

*/INTER1/EIJ(19,19),UIJ(19,19),KIJ(19,19),GIJ(19,19)

*/KM/KM/COEF1/B1(13),C1(53)/AN/AN(53)

*/GQFN/GN(53),QN(53),FN(53)/UN/UN(53)

DO 1 I=1,NC

EIJM(I,I)=EI(I)

GIJM(I,I)=GI(I)

DO 1 J=1,NC

IF(I.GE.J) GO TO 1

EIJM(I,J)=EIJ(I,J)*(EI(I)*EI(J))**.5

GIJM(I,J)=GIJ(I,J)*(GI(I)+GI(J))/2.

1 CONTINUE

KM=0D0

UM=0D0

KM=0D0

UM=0D0

GM=0D0

QM=0D0

FM=0D0

DO 3 I=1,NC

KM=KM+Y(I)*KI(I)**2.5

UM=UM+Y(I)*EI(I)**2.5

GM=GM+Y(I)*GI(I)

QM=QM+Y(I)*QI(I)

3 FM=FM+Y(I)**2*FI(I)

KM=KM*KM

UM=UM*UM

DO 5 I=1,NC-1

DO 5 J=I+1,NC

UM=UM+2.*Y(I)*Y(J)*(UIJ(I,J)**5-1D0)*(EI(I)*EI(J))**2.5

GM=GM+2.*Y(I)*Y(J)*(GIJ(I,J)-1D0)*(GI(I)+GI(J))

5 KM=KM+2.*Y(I)*Y(J)*(KIJ(I,J)**5-1D0)*(KI(I)*KI(J))**2.5

KM=KM**.6

UM=UM**.2

DO 7 N=1,13

B1(N)=0D0

DO 9 I=1,NC

9 1(N)=B1(N)+Y(I)*Y(I)(GIJM(I,I)+ 1D0-GN(N))**GN(N)*

*(QI(I)*QI(I)+1D0-QN(N))**QN(N)*(FI(I)+1D0-FN(N))*FN(N)*

*EIJM(I,I)"UN(N)*KI(I)*KI(I)*KI(I)

DO 11 I=1,NC-1

DO 11 J=I+1,NC

11 1(N)=B1(N)+2.*Y(I)*Y(J)(GIJ(I,J)+1D0-GN(N))**GN(N)*

*(QI(I)*QI(J)+1D0-QN(N))**QN(N)*((FI(I)*FI(J))**.5+

1D0-FN(N))**FN(N)*EIJM(I,J)**UN(N)*(KI(I)*KI(J)**1.5

7 CONTINUE

DO 13 N=8,53

13 C1(N)=AN(N)*(GM+1D0-GN(N))**GN(N)*(QM**2+1D0-QN(N))**

*QN(N)*(FM+1D0-FN(N))**FN(N)*UM**UN(N)

RETURN

END

SUBROUTINE PHASE1

IMPLICIT REAL*8(A-H,O-Z)

COMMON/Z/Z/RM/RM/T/T/P/P/AI1/AO,A1/AN/AN(53)

*/COEF1/B1(13),C1(53)/COEF2/B,C(53)/UN/UN(53)

CALL PCONT1(P,T)

IF(Z.EQ.0D0) GO TO 134

B=0D0

DO 1 N=1,13

1 B=B+AN(N)/T**UN(N)*B1(N)

DO 3 N=8,53

3 C(N)=C1(N)/T**UN(N)

PR=P/5.

RO=9D3*P/(RM*T*(1.1*PR+0.7D0))

CALL FUN1(RO)

Z=1D0+AO

134 RETURN

END

,

( )

SUBROUTINE FUN1(X)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/P/P/RM/RM/T/T/AI1/AO,A1

ITER=1

1 CONTINUE

CALL COMPL1(X)

Z=1.D0+AO

FX= 1 .D6*(P-(1.D-3*RM*T*Z*X))

F= 1 .D3*RM*(1.D0+A1)

DR=FX/F

X=X+DR

IF(ITER.GT.10) GO TO 4

ITER=ITER+1

IF(DABS(DR/X).GT.1.D-6) GO TO 1

4 CALL COMPL1(X)

RETURN

END

SUBROUTINE PCONT1(P,T)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/Z/Z

Z=-1D0

IF(T.LT.250D0.OR.T.GT.340D0)Z=0D0

IF(P.LE.0D0.OR.P.GT.12D0) Z=0D0

RETURN

END

SUBROUTINE COMPL1(RO)

IMPLICIT REAL*8(A-H,O-Z)

REAL*8 KM

INTEGER BN,CN

COMMON/KM/KM/COEF2/B,C(53)/BCKN/BN(53),CN(53),KN(53)/AI1/AO,A1

ROR=KM*RO

S1=0D0

S2=0D0

S3=0D0

DO 1 N=8,53

EXP=DEXP(-CN(N)*ROR**KN(N))

IF(N.LE.13) S1=S1+C(N)

S2=S2+C(N)*(BN(N)-CN(N)*KN(N)*ROR**KN(N))*ROR**BN(N)*EXP 1

S3=S3+C(N)*(-CN(N)*KN(N)**2*KM*ROR**(KN(N)-1)*ROR**BN(N)*

*EXP+(BN(N)-CN(N)*KN(N)*ROR**KN(N))*BN(N)*KM*ROR**(BN(N)-1)*

*EXP-(BN(N)-CN(N)*KN(N)*ROR**KN(N))*ROR**BN(N)*EXP*CN(N)*KN(N)*

*KM*ROR**(KN(N)-1))AO1=B*RO-ROR*S1

AO=AO1+S2

A1=AO+AO1+RO*S3

RETURN

END

BLOCK DATA DCAGA8

IMPLICIT REAL*8(A-H,O-Z)

REAL*8 KD

INTEGER BN,CN,GN,QN,FN

COMMON/EFD/ED(19),KD(19),GD(19),QD(19),FD(19)

*/BCKN/BN(53),CN(53),KN(53)/UN/UN(53)

*/AN/AN(53)/GQFN/GN(53),QN(53),FN(53)

DATA ED/151.3183D0,244.1667D0,298.1183D0,337.6389D0,324.0689D0,

*99.73778D0,241.9606D0,296.355D0,370.6823D0,365.5999D0,

*402.8429D0,427.5391D0,450.6472D0,472.1194D0:488.7633D0,

*2.610111D0,26.95794D0,105.5348D0,122.7667D0/

DATA KD/.4619255D0,.5279209D0,.583749D0,.6341423D0,.6406937D0,

*.4479153D0,.4557489D0,.4618263D0,.6798307D0,.6738577D0,

*.7139987D0,.7503628D0,.7851933D0,.8157596D0,.8389542D0,

*.3589888D0,.3514916D0,.4533894D0,.4186954D0/

DATA GD/0D0,.0793D0,.141239D0,.281835D0,.256692D0,.027815D0,

*.189065D0,.0885D0,.366911D0,.332267D0,.432254D0,.512507D0,

*.576242D0,.648601D0,.716574D0,0D0,.034369D0,.038953D0,.021D0/

DATA QD/6*0D0,.69D0,12*0D0/,FD/16*0D0,1D0,2*0D0/

DATA AN/.1538326D0,1.341953D0,-2.998583D0,-.04831228D0,

*.3757965D0,-1.589575D0,-.05358847D0,2.29129D-9,1576724D0,

*-.4363864D0,-.04408159D0,-.003433888D0,.03205905D0,.02487355D0, -

*.07332279D0,-.001600573D0,.6424706D0,-.4162601D0,-.06689957D0,

*.2791795D0,-.6966051D0,-.002860589D0,-.008098836D0,3.150547D0,

*.007224479D0,-.7057529D0,.5349792D0,-.07931491D0-1.418465D0,

*-5.99905D-17,.1058402D0,.03431729D0,-.007022847D0,.02495587D0,

*.04296818D0,.7465453D0,-.2919613D0,7.294616D0,-9.936757D0,

*-.005399808D0,-.2432567D0,.04987016D0,.003733797D0,1.874951D0,

*.002168144D0,-.6587164D0,.000205518D0,.009776195D0,-.02048708D0,

*.01557322D0,.006862415D0,-.001226752D0,.002850906D0/

DATA BN/13*1,9*2,10*3,7*4,5*5,2*6,2*7,3*8,2*9/

DATA CN/7*0,6*1,2*0,7*1,0,9*1,2*0,5*1,0,4*1,0,1,0,6*1/

DATA KN/7*0,3,3*2,2*4,2*0,3*2,4*4,0,2*1,2*2,2*3,3*4,2*0,3*2,

*2*4,0,2*2,2*4,0,2,0,2,1,4*2/

DATA UN/0D0,.5D0,1D0,3.5D0,-.5D0,4.5D0,.5D0,-6D0,2D0,3D0,2*2D0,

*11D0,-.5D0,.5D0,0D0,4D0,6D0,21D0,23D0,22D0,-1D0,-.5D0,7D0,-1D0,

*6D0,4D0,1D0,9D0,-13D0,21D0,8D0,-.5D0,0D0,2D0,7D0,9D0,22D0,23D0,

*1D0,9D0,3D0,8D0,23D0,1.5D0,5D0,-.5D0,4D0,7D0,3D0,0D0,1D0,0D0/

DATA GN/4*0,2*1,13*0,1,3*0,1,2*0,3*1,16*0,1,2*0,1,0,1,2*0/

DATA QN/6*0,1,3*0,1,9*0,1,0,1,8*0,1,4*0,1,4*0,1,0,1,2*0,1,5*0,1/

DATA FN/7*0,1,13*0,1,2*0,1,4*0,1,23*0/

END

C ***********************************************************

* *

* *

* . *

* *

C ***********************************************************

SUBROUTINE VNIC(ICALC)

IMPLICIT REAL*8(A-H,O-Z)

REAL*8 LIJ(8,8)

DIMENSION VC(8),TC(8),PII(8),DIJ(8,8)

COMMON/PARCD/VCD(8),TCD(8),PIID(8)/ABIJ/AIJ(10,8),BIJ(10,8)

*/B/B(10,8)/RM/RM/Y/Y(8)/BM/BM(8)/NI/NI(8)/NC/NC/Z/Z

RM=8.31451DO

IF(ICALC.NE.1) GO TO 1

CALL COMPON

IF(Z.EQ.0D0) GO TO 133

CALL DDIJ(DIJ,LIJ)

DO 75 I=1,NC

TC(I)=TCD(NI(I))

VC(I)=BM(I)/VCD(NI(I))

PII(I)=PIID(NI(I))

DO 123 J=1,NC

IF(I.GE.J) GO TO 123

DIJ(I,J)=DIJ(NI(I),NI(J))

LIJ(I,J)=LIJ(NI(I),NI(J))

123 CONTINUE

75 CONTINUE

CALL PARMIX(DIJ,LIJ,TC,VC,PII,PIM)

DO 27 I=1,10

DO 27 J=l,8

27 B(I,J)=AIJ(I,J)+BIJ(I,J)*PIM

1 CALL PHASE

133 RETURN

END

SUBROUTINE COMPON

IMPLICIT REAL*8(A-H,O-Z)

DIMENSION BMI(25),ROI(8),GI(8),YI(25)

COMMON/Y/Y(8)/BMM/BMM/BM/BM(8)/YI/YC(25)/NI/NI(8)/NC/NC/NPR/NPR

DATA BMI/16.043D0,30.07D0,44.097D0,2*58.123D0,28.0135D0,

*44.01D0,34.082D0,26.038D0,28.054D0,42.081D0,3*72.15D0,

*86.177D0,78.114D0,100.204D0,92.141D0,114.231D0,128.259D0,

*142.286D0,4.0026D0,2.0159D0,28.01D0,31.9988D0/

DATAR0I/0.6682D0,1.2601D0,1.8641D0,2.4956D0,2.488D0,

*1.1649D0,1.8393D0,1.4311D0/

DO 100 I=1,25

100 YI(I)=YC(I)

IF(NPR.EQ.1) GO TO 333

BMM=0D0

DO 3333 I=1,25

3333 =+YI(I)*I(I)

333 YS=0D0

DO 55 I=9,25

YS=YS+YI(I)

55 CONTINUE

YS1=0D0

DO 67 I=12,21

67 YS1=YS1+YI(I)

YS2=0D0

DO 69 1=22,25

69 YS2=YS2+YI(I)

YI(2)=YI(2)+YI(9)+YI(10)

YI(3)=YI(3)+YI(11)

YI(4)=YI(4)+YS1

YS3=YI(4)+YI(5)

IF(NPR.EQ.1.AND.YI(5).LT.0.01D0.AND.YS3.LT.0.03D0) YI(4)=YS3

IF(NPR.EQ.1.AND.YI(5).LT.0.01D0.AND.YS3.LT.0.03D0) YI(5)=0D0

IF(NPR.EQ.0.AND.Y1(5).LT.0.01D0.AND.YS3.LE.0.03D0) YI(4)=YS3

IF(NPR.EQ.0.AND.YI(5).LT.0.01D0.AND.YS3.LE.0.03D0)YI(5)=0D0

YI(6)=YI(6)+YS2

IF(NPR.EQ.0) GO TO 555

ROM=0D0

DO 7 I=1,8

7 ROM=ROM+YI(I)*ROI(I)

DO 9 I=1,8

9 GI(I)=YI(I)*ROI(I)/ROM

SUM=0D0

DO 11 1=1,8

11 SUM=SUM+GI(I)/BMI(I)

SUM=1./SUM

DO 13 I=1,8

13 YI(I)=GI(I)*SUM/BMI(I)

555 NC=0

YSUM=0D0

DO 155 I=1,8

IF(YI(I).EQ.0D0) GO TO 155

NC=NC+1

NI(NC)=I

Y(NC)=YI(I)

YSUM=YSUM+Y(NC)

BM(NC)=BMI(I)

155 CONTINUE

CALL MOLDOL(YI,YS)

DO 551 I=1,NC

551 Y(I)=Y(I)/YSUM

RETURN

END

SUBROUTINE MOLDOL(YI,YS)

IMPLICIT REAL*8(A-H,O-Z)

DIMENSION YI(25)

COMMON/Z/Z

Z=-1D0

IF(YI(1).LT.0.65D0.OR.YI(2).GT.0.15D0.OR.YI(3).GT.0.035D0.OR.

*YI(4).GT.0.015D0.OR.YI(5).GT.0.015D0.OR.YS.GT.0.01D0)Z=0D0

IF(YI(6).GT.0.2D0.OR.YI(7).GT.0.15D0.OR.YI(8).GT.0.3D0)Z=0D0

RETURN

END

SUBROUTINE DDIJ(DIJ,LIJ)

IMPLICIT REAL-8(A-H,O-Z)

REAL*8 LIJ(8,8)

DIMENSION DIJ(8,8)

DO 1 I=1,8

DO 1 J=l,8

LIJ(I,J)=0.D0

1 DIJ(I,J)=0.D0

DIJ(1,2)=0.036D0

DIJ(1,3)=0.076D0

DIJ(1,4)=0.121D0

DIJ(1,5)=0.129D0

DIJ(1,6)=0.06D0

DIJ(1,7)=0.074D0

DIJ(2,6)=0.106D0

DIJ(2,7)=0.093D0

DIJ(6,7)=0.022D0

DIJ(1,8)=0.089D0

DIJ(2,8)=0.079D0

DU(6,8)=0.211D0

DIJ(7,8)=0.089D0

LIJ(1,2)=-0.074D0

LIJ(1,3)=-0.146D0

LIJ(1,4)=-0.258D0

LIJ(1,5)=-0.222D0

LIJ(1,6)=-0.023D0

LIJ(1,7)=-0.086D0

LIJ(6,7)=-0.064D0

LIJ(7,8)=-0.062D0

RETURN

END

SUBROUTINE PARMIX(DIJ,LIJ,TC,VC,PII,PIM)

IMPLICIT REAL*8(A-H,O-Z)

REAL*8 LIJ(8,8)

DIMENSION Y(8),DIJ(8,8),VCIJ(8,8),TCIJ(8,8),V13(8),TC(8),VC(8),

*PII(8),PIIJ(8,8)

COMMON/PARCM/TCM,VCM/Y/Y/NC/NC/PCM/PCM

DO 1 I=1,NC

1 V13(I)=VC(I)**(1.DO/3.DO)

DO 3 I=1,NC

VCIJ(I,I)=VC(I)

PIIJ(I,I)=PII(I)

TCIJ(I,I)=TC(I)

DO 3 J=1,NC

IF(I.GE.J) GO TO 3

VCIJ(I,J)=(1.DO-LIJ(I,J))*((V13(I)+V13(J))/2.)**3

PIIJ(I,J)=(VC(I)*PII(I)+VC(J)*PII(J))/(VC(I)+VC(J))

TCU(I,J)=(1.D0-DIJ(I,J))*(TC(I)*TC(J))**0.5

VCIJ(J,I)=VCIJ(I,J)

PIIJ(J,I)=PIIJ(I,J)

TCIJ(J,I)=TCIJ(I,J)

3 CONTINUE

VCM=0.D0

PIM=0.D0

TCM=0.D0

DO 5 I=1,NC

DO 5 J=1,NC

VCM=VCM+Y(I)*Y(J)*VCIJ(I,J)

PIM=PIM+Y(I)*Y(J)*VCIJ(I,J)*PIIJ(I,J)

5 TCM=TCM+Y(I)*Y(J)*VCIJ(I,J)*TCIJ(I,J)**2

PIM=PIM/VCM

TCM=(TCM/VCM)**0.5

PCM=8.31451D-3*(0.28707D0-0.05559*PIM)*TCM/VCM

RETURN

END

SUBROUTINE PHASE

IMPLICIT REAL*8(A-H,O-Z)

COMMON/Z/Z/RM/RM/T/T/P/P/PCM/PCM/AI/AO,A1

IF(T.LT.250D0.OR.T.GT.340D0.OR.P.LE.0D0.OR.P.GT.12D0) THEN

Z=0D0

GO TO 134

ENDIF

PR=P/PCM

RO=9D3*P/(RM*T*(1.1*PR+0.7D0))

CALL FUN(RO)

CALL OMTAU(RO,T)

IF(Z.EQ.0D0) GO TO 134

Z=1.D0+AO

134 RETURN

END

,

( )

SUBROUTINE FUN(X)

IMPLICIT REAL*8(A-H,-Z)

COMMON/P/P/RM/RM/T/T/AI/AO,A1

ITER=1

1 CONTINUE

NPRIZ=0

IF(ITER.NE.l) NPRIZ=1

CALL COMPL(X,T,NPRIZ)

Z=1.D0+AO

FX=1.D6*(P-(1.D-3*RM*T*Z*X))

F=1.D3*RM*T*(1.D0+A1)

DR=FX/F

X=X+DR

IF(ITER.GT.10) GO TO 4

ITER=ITER+1

IF(DABS(DR/X).GT.1.D-6) GO TO 1

4 CALL COMPL(X,T,NPRIZ)

RETURN

END

SUBROUTINE OMTAU(RO,T)

IMPLICIT REAL*8(A-H,O-Z)

COMMON/PARCM/TCM,VCM/Z/Z

Z=-1D0

TR=T/TCM

ROR=RO*VCM

IF(TR.LT.1.05D0) Z=0D0

IF(ROR.LT.0.D0.OR.ROR.GT.3.D0) Z=0D0

RETURN

END

SUBROUTINE COMPL(RO,T,NPRIZ)

IMPLICIT REAL*8(A-H,O-Z)

DIMENSION B(10,8),BK(10)

COMMON/PARCM/TCM,VCM/B/B/AI/AO,A1

IF(NPRIZ.NE.0) GO TO 7

TR=T/TCM

DO 1 I=1,10

BK(I)=0

DO 1 J=1,8

1 BK(I)=BK(I)+B(I,J)/TR**(J-1)

7 ROR=RO*VCM

AO=0.D0

A1=0.D0

DO 33 I=1,10

D=BK(I)*ROR**I

AO=AO+D

33 A1=A1+(I+1)*D

RETURN

END

BLOCK DATA BDVNIC

IMPLICIT REAL*8(A-H,O-Z)

COMMON/PARCD/VCD(8),TCD(8),PIID(8)/ABIJ/AIJ(10,8),BIJ(10,8)

DATA TCD/190.67D0,305.57D0,369.96D0,425.4D0,407.96D0,

*125.65D0,304.11D0,373.18D0/

DATA VCD/163.03D0,205.53D0,218.54D0,226.69D0,225.64D0,

*315.36D0,466.74D0,349.37D0/

DATA PIID/0.0006467D0,0.1103D0,0.1764D0,0.2213D0,0.2162D0,

*0.04185D0,0.2203D0,0.042686D0/

DATA AIJ/.6087766D0,-.4596885D0,1.14934D0,-.607501D0,

*-.894094D0,1.144404D0,-.34579D0,-.1235682D0,.1098875D0,

*-.219306D-1,-1.832916D0,4.175759D0,-9.404549D0,10.62713D0,

*-3.080591D0,-2.122525D0,1.781466D0,-.4303578D0,-.4963321D-1,

*.347496D-1,1.317145D0,-10.73657D0,23.95808 D0,-31.47929D0,

* 18.42846D0,-4.092685D0,-. 1906595D0,.4015072D0,-.1016264D0,

*-.9129047D-2,-2.837908D0,15.34274D0,-27.71885D0,35.11413D0,

*-23.485D0,7.767802D0,-1.677977D0,.3157961D0,.4008579D-2,0.D0,

*2.606878D0,-11.06722D0,12.79987D0,-12.11554D0,7.580666D0,

*-1.894086D0,4*0.D0,

*-1.15575D0,3.601316D0,-.7326041D0,-1.151685D0,.5403439D0,

*5*0.D0,.9060572D-1,-.5151915D0,.7622076D-1,7*0.D0,

*.4507142D-1,9*0.D0/

DATA BIJ/-.7187864D0,10.67179D0,-25.7687D0,17.13395D0,

*16.17303D0,-24.38953D0,7.156029D0,3.350294D0,-2.806204D0,

*.5728541D0,6.057018D0,-79.47685D0,216.7887D0,-244.732D0,

*78.04753D0,48.70601D0,-41.92715D0,10.00706D0,1.237872D0,

*-.8610273D0,-12.95347D0,220.839D0,-586.4596D0,744.4021D0,

*-447.0704D0,99.6537D0,5.136013D0,-9.5769D0,2.41965D0,

*.2275036D0,15.71955D0,-302.0599D0,684.5968D0,-828.1484D0,

*560.0892D0,-185.9581D0,39.91057D0,-7.567516D0,-.1062596D0,

*0.D0,-13.75957D0,205.541D0,-325.2751D0,284.6518D0,

*-180.8168D0,46.05637D0,4*0.D0,

*6.466081D0,-57.3922D0,36.94793D0,20.77675D0,-12.56783D0,

*5*0.D0,-.9775244D0,2.612338D0,-.4059629D0,7*0.D0,

*-.2298833D0,9*0.D0/

END

()

.1 NX19

0,101325 293,15 : 0,6799 /3

:

........................................................................................................ 0,8858 . %

................................................................................. 0,0668 . %

................................................................................................. 2,001

........................................................................................... 270,00

................................................................. 0,9520

................................................................................................. 2,494

........................................................................................... 280,00

................................................................. 0,9473

................................................................................................. 0,900

........................................................................................... 290,00

................................................................. 0,9844

.2 GERG-91

0,101325 293,15 : 0,6799 /3

:

........................................................................................................ 0,8858 . %

................................................................................. 0,0668 . %

................................................................................................. 2,001

........................................................................................... 270,00

................................................................. 0,9521

................................................................................................. 3,997

........................................................................................... 290,00

................................................................. 0,9262

................................................................................................. 7,503

........................................................................................... 330,00

................................................................. 0,9244

.3 AGA8-92DC

:

....................................................................................................... 98,2722

......................................................................................................... 0,5159

..................................................................................................... 0,1607

- .................................................................................................... 0,0592

.......................................................................................................... 0,8858

................................................................................... 0,0668

- .................................................................................................. 0,0157

- .................................................................................................. 0,0055

- .................................................................................................. 0,0016

- .................................................................................................... 0,0009

....................................................................................................... 0,0157

0,101325 293,15 : 0,6799 /3

................................................................................................. 2,001

........................................................................................... 270,00

................................................................. 0,9520

................................................................................................. 3,997

........................................................................................... 290,00

................................................................. 0,9262

................................................................................................. 7,503

........................................................................................... 330,00

................................................................. 0,9246

.4

:

....................................................................................................... 89,2700

......................................................................................................... 2,2600

..................................................................................................... 1,0600

- .................................................................................................... 0,0100

.......................................................................................................... 0,0400

................................................................................... 4,3000

............................................................................................ 3,0500

................................................................................................ 0,0100

0,101325 293,15 : 0,7675 /3

................................................................................................. 1,081

........................................................................................... 323,15

................................................................. 0,9853

................................................................................................. 4,869

........................................................................................... 323,15

................................................................. 0,9302

................................................................................................. 9,950

........................................................................................... 323,15

................................................................. 0,8709

()

( )

.1 NX19

( )

, %

,

1,991

2,011

1,00

,

269,50

270,50

0,35

, /3 (0,101325 , 293,15 )

0,6790

0,6808

0,25

, . %:

 

 

 

(N2)

0,8769

0,8947

2,00

(2)

0,0661

0,0675

2,00

( ) - 0,9520

: (82) - 0,09 %; (86) - 0,07 %.

.2 GERG-91

( )

, %

,

1,991

2,011

1,00

,

269,50

270,50

0,35

, /3 (0,101325 , 293,15 )

0,6790

0,6808

0,25

, . %:

 

 

 

(N2)

0,8769

0,8947

2,00

(2)

0,0661

0,0675

2,00

( ) - 0,9521

: (82) - 0,09 %; (86) - 0,09 %.

.3 AGA8-92DC

( )

, %

,

1,991

2,011

1,00

,

269,50

270,50

0,35

, . %:

 

 

 

(4)

97,2722

99,2722

2,00

(26)

0,5030

0,5288

5,00

(38)

0,1607

0,1607

-

- (-410)

0,0592

0,0592

-

(N2)

0,8769

0,8947

2,00

(2)

0,0661

0,0675

2,00

- (-512)

0,0157

0,0157

-

- (-614)

0,0055

0,0055

-

- (-716)

0,0016

0,0016

-

- (-C8H18)

0,0009

0,0009

-

()

0,0157

0,0157

-

( ) - 0,9520

- 0,08 %

.4

( )

, %

,

1,076

1,086

1,00

,

322,65

323,65

0,31

, . %:

 

 

 

(4)

88,3700

90,1700

2,00

(26)

2,2030

2,3170

5,00

(C3H8)

1,0600

1,0600

-

- (-410)

0,0100

0,0100

-

(N2)

0,0396

0,0404

2,00

(2)

4,2570

4,3430

2,00

(H2S)

3,0500

3,0500

-

(36)

0,0100

0,0100

-

( ) - 0,9853

- 0,03 %

()

[1] .. . . - ., - , 1979, 348

[2] Kleinrahm R., Duschek W., Wagner W. Measurement and correlation of the (pressure, density, temperature) relation of methane in the temperature range from 273.15 to 323.15 at pressures up to 8 MPa. - J. Chem. Thermodynamics, 1988, v.20, p.621-631

[3] Robinson R.L., Jacoby R.H. Better compressibility factors. - Hydrocarbon Processing, 1965,v.44,No.4,p.141-145

[4] Achtermann H.-J., Klobasa F.,Rogener H. Realgasfaktoren von Erdgasen. Teil I: Bestimmung von Realgasfaktoren aus Brechungsindex-Messungen. - Brennstoff-Warme-Kraft, 1982, Bd.34, No.5, s.266-271

[5] Achtermann H.-J., Klobasa F.,Rogener H. Realgasfaktoren von Erdgasen. Teil II: Bestimmung von Realgasfaktoren mit eener Burnett-Apparatur. - Brennstoff-Warme-Kraft, 1982, Bd.34, No.6, s.311-314

[6] Eubank Ph.T., Scheloske J., Hall K.R., Holste J.C. Densities and mixture virial coefficients for wet natural gas mixtures. - Journal of Chemical and Engineering Data, 1987, v.32, No.2, p.230-233

[7] Jaeschke ., Julicher H.P. Realgasfaktoren von Erdgasen. Bestimmung von Realgasfaktoren nach der Expansionsmethode. - Brennstoff-Warme-Kraft, 1984, Bd.36, No.11, s.445-451

[8] Jaeschke . Realgasverhalten Einheitliche Berechnungsmoglichkeiten von Erdgas L und H. - Gas und Wasserfach. Gas/Erdgas, 1988, v.129, No.l, s.30-37

[9] Blanke W., Weiss R. pvT-Eigenschaften und Adsorptions- verhalten von Erdgas bei Temperaturen zwischen 260 und 330 mit Drucken bis 3 MPa. - Erdol-Erdgas-Kohle, 1988, Bd.104, H.10, s.412-417

[10] Samirendra N.B. et al Compressibility Isotherms of Simulated Natural Gases. - J. Chem. Eng. Data, 1990, v.35, No.l, p.35-38

[11] Fitzgerald M.P., Sutton C.M. Measurements of Kapuni and Maui natural gas compressibility factors and comparison with calculated values. - New Zealand Journal of Technology, 1987, v.3, No.4, p.215-218

[12] Jaeschke ., Humphreys A.E. The GERG Databank of High Accuracy Compressibility Factor Measurements. GERG TM4 1990. - GERG Technical Monograph, 1990, 477 p

[13] Jaeschke ., Humphreys A.E. Standard GERG Virial Equation for Field Use. Simplification of the Input Data Requirements for the GERG Virial Equation - an Alternative Means of Compressibility Factor Calculation for Natural Gases and Similar Mixtures. GERG TM5 1991. - GERG Technical Monograph, 1991, 173 p

[14] ISO/TC 193 SC1 63. Natural gas - calculation of compression factor. Part 3: Calculation using measured physical properties

[15] ISO/TC 193 SC1 62. Natural gas - calculation of compression factor. Part 2: Calculation using a molar composition analysis

[16] ISO 5168:1978 International Standard. Measurement of fluid flow - Estimation of uncertainty of a flow-rate measurement

[17] VDI/VDE 2040, part 2, 1987. Calculation principles for measurement of fluid flow using orifice plates, nozzles and venturi tubes. Equations and formulas

[18] Jaeschke . et al. High Accuracy Compressibility Factor Calculation for Natural Gases and Similar Mixtures by Use of a Truncated Virial Equation. GERG TM2 1988. - GERG Technical Monograph, 1988, 163 p

: , , , , , , , , , , , , ,

 >>>     
   [ − ]  
: (18.04.1986)
: 18 1986 .
  "    -         "
: (21.03.1996)
: 1 1997 .
   127.2-93 " .  "
: (18.05.1982)
: 1 1983 .
   25.504-82* "    .     "

    :  140 2 1

    9 :  910      153.5 4 1 910

    -11 :    913   201 2 2 913

    -1 :   910 139 3 2 910

    :   114 2 1


 @Mail.ru   adm@stroyplan.ru.
- cc, , .
© 2006-2016 ""
     .