30
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