Effects of process-generated hydrogen on RPV walls

21

continue pause 'MAXITS exceeded in newt' END

SUBROUTINE funcv(n,y,ff) IMPLICIT REAL *8(A-H,O-Z)

COMMON /BLK01/ TC,TK,R,F,FFF,AV,PH,H2 COMMON /BLK03/ WK,DW,WKM COMMON /BLK23/ CK(50),D0(25),D(25),E(25),V(25) COMMON /Q/ Q(10),TLI,TB DIMENSION Y(n),FF(n)

ff(1) = y(3) + y(4) + y(5) - TLI ff(2) = y(5) + y(6) + y(7) + 2.*y(8) + 3.*y(9) + 4.*y(10) -TB ff(3) = y(1) +y(3) - y(2) - y(7) - y(8) - y(9) - 2.*Y(10)

ff(4) = y(3)*y(2)- y(4)/Q(6) ff(5) = y(3)*y(7) - y(5)/Q(7) ff(6) = y(1)*y(2) - WKM ff(7) = y(6)*y(2) - y(7)/Q(1)

ff(8) = y(6)**2 *y(2) - y(8)/Q(2) ff(9) = y(6)**3 *y(2) - y(9)/Q(3) ff(10) = y(6)**4 *y(2)**2 - y(10)/Q(4) DO 1 k=1,n ff(k) = ff(k)

1 continue return end

SUBROUTINE lnsrch(n,xold,fold,g,p,x,f,stpmax,check,func) IMPLICIT REAL *8(A-H,O-Z) INTEGER n LOGICAL check REAL *8 f,fold,stpmax,g(n),p(n),x(n),xold(n),func,ALF,TOLX PARAMETER (ALF=1.e-4,TOLX=1.e-10) EXTERNAL func CU USES func INTEGER i REAL a,alam,alam2,alamin,b,disc,f2,fold2,rhs1,rhs2,slope,sum *,temp,test,tmplam check= .false. sum=0. do 11 i=1,n sum=sum+p(i)*p(i) 11 continue sum=sqrt(sum) if (sum .gt. stpmax) then do 12 i=1,n p(i)=p(i)*stpmax/sum 12 continue endif slope=0. do 13 i=1,n slope=slope+g(i)*p(i)