30
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