Effects of process-generated hydrogen on RPV walls

*H)4; 15 - B(OH)3;'

write (4,*) '16 - B(OH)4-; 17 - B2(OH)7-; 18 - B3(OH)10-; 19 - B3(O *H)14--;'

9999 CONTINUE

pH = -log10(Y(4)) O_R = (Y(9)+y(6)+y(3))/(Y(8)+Y(2)+Y(7))

print *, pH,O_R close (3) close (4) close (5)

END

SUBROUTINE RK4SYS(N,X,Y,H,NSTEP) IMPLICIT REAL *8(A-H,O-Z) COMMON /LP/ LP(6)

DIMENSION YY(20),Y(20),F1(20),F2(20),F3(20),F4(20) 94 FORMAT ( 20E12.4)

WRITE (4,94) X,(Y(LP(k1)),K1=1,6) WRITE (3,94) x,(Y(k1),K1=1,19)

step = H tfig = h H2 = 0.5*H

DO 6 K=1,NSTEP

CALL XPSYS(Y,F1) DO 2 I=1,N YY(I) = Y(I) + H2*F1(I) CONTINUE CALL XPSYS(YY,F2) DO 3 I=1,N YY(I) = Y(I) + H2*F2(I) 3 CONTINUE CALL XPSYS(YY,F3) DO 4 I=1,N YY(I) = Y(I) + H*F3(I) 4 CONTINUE 2

CALL XPSYS(YY,F4) DO 5 I=1,N Y(I) = Y(I) + H*(F1(I) + 2.0*(F2(I)+F3(I))+F4(I))/6.0 5 CONTINUE C Y(4) = Y(1)+Y(5)+Y(7)+ Y(7)+ Y(10) - Y(12) C > + Y(16) + Y(17)+ Y(18) +2.*Y(19)

x = y(20) IF (x .ge. tfig) IF (x .ge. tfig)

WRITE (4,94) X,(Y(LP(k1)),K1=1,6)

WRITE (3,94) x,(Y(k1),K1=1,19) IF (x .ge. tfig) WRITE (*,'(1x,a5,E12.4,a3)') 't = ', X, ' s'