Effects of process-generated hydrogen on RPV walls

DO 21 K=1,N

print *, k, V(K),VV(k)

21 continue

3333

continue

c

y= 0

y(1) = 0. y(2) = 0. y(3) = 0. y(4) = V(1) y(5) = V(2) y(6) = O2 y(7) = 0. y(8) = H2 y(9) = H2O2 y(10) = 0. y(11) = 0.

DO 44 K=1,8 y(11+k) = V(2+k)

44 continue 444 CONTINUE

DX= Tmax H= dt NX = DX/H

X=0.

Y(20) = X

CALL RK4SYS(20,X,Y,H,NX)

C Y(4) = Y(1)+Y(5)+Y(7)+ Y(10) - Y(12) C > + Y(16) + Y(17)+ Y(18) +2.*Y(19)

elect = 0. DO 66 i=1, Jmax elect = elect +Z(i)*Y(i)

66 continue PRINT *, ELECT

PRINT *, '[H]=', Y(2)

c c

STOP '5555'

write (5,94) Y(2)

WRITE (*,*) 'THE RESULTS ARE WRITTEN IN "PWR_CHEM_RK.out"' write (4,*) write (4,*) '1 - e(aq); 2 - H; 3 - OH; 4 - H+; 5 - OH-; 6 - O2; * 7 - O2-; 8 - H2; 9 - H2O2;'

write (4,*) '10 - HO2-; 11 - HO2; 12 - Li+; 13 - LiOH; 14 - LiB(O