36
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