Effects of process-generated hydrogen on RPV walls
Z(7) = -1. Z(8) = 0. Z(9) = 0 Z(10) = -1. Z(11) = 0. Z(12) = 1. Z(13) = 0. Z(14) = 0. Z(15) = 0. Z(16) = -1. Z(17) = -1. Z(18) = -1. Z(19) = -2.
C c c
V = 0 VV = 0 N=10
DO 20 K=1,N V(K)= 1.e-0 VV(k) = 0.
20 continue
call newt(V,N,check) call funcv(N,V,y) if (check) then
write (*,*) 'Convergence problems.' endif
c
go to 3333
QQ(1)=10.**(1573.21/TK+28.6059+0.012078*TK-13.2258*LOG10(TK)) QQ(2)=10.**(2756.1/TK-18.966+5.835*LOG10(TK)) QQ(3)=10.**(3339.5/TK-8.084+1.497*LOG10(TK)) QQ(4)=10.**(12820/TK-134.56+42.105*LOG10(TK)) QQ(5)=10.**(14099/TK-118.115+36.237*LOG10(TK)) QQ(6) = 1./(10**(-0.7532 - 0.0048*Tc + 6.746e-6*Tc**2)) c Q(6)=1.99 QQ(7)=2.12
VV(3) = TLi VV(6) = TB VV(2) = TLi/(1. + QQ(1)*TB)
VV(1) = WKM/VV(2) VV(7) = TLi - VV(2)
VV(4) = VV(2)*VV(3)*QQ(6) VV(5) = VV(3)*VV(7)*QQ(7) VV(8) = VV(6)**2*VV(2)*QQ(2) VV(9) =VV(6)**3*VV(2)*QQ(3) VV(10) = VV(6)**4*VV(2)**2*QQ(4)
49 continue
Made with FlippingBook Publishing Software