Effects of process-generated hydrogen on RPV walls

Galpha(1) = 0.06 Galpha(2) = 0.21 Galpha(3) = 0.24 Galpha(4) = 0.06 Galpha(5) = 0. Galpha(6) = 0. Galpha(7) = 0. Galpha(8) = 1.3 Galpha(9) = 0.985 Galpha(10) = 0.0 Galpha(11) = 0.22

c GAMN = 1.e20 c GAMR = 3.e21

S1 = 0. S2 = 0. DO I=1,11

E(I)=(GAMn*GN(I)/100./AV+GAMg*GRA(I)/100./AV)*1000. > + (GAMa*Galpha(I)/100./AV)*1000.

E(i) = E(I)*DW

ENDDO

ccc

stop 'property'

RETURN

END SUBROUTINE newt(x,n,check)

IMPLICIT REAL *8(A-H,O-Z) INTEGER n,nn,NP,MAXITS LOGICAL check REAL *8 x(n),fvec,TOLF,TOLMIN,TOLX,STPMX PARAMETER (NP=40,MAXITS=200,TOLF=1.e-16,TOLMIN=1.e-16,TOLX=1.e-16, *STPMX=100.) COMMON /newtv/ fvec(NP),nn SAVE /newtv/ CU USES fdjac,fmin,lnsrch,lubksb,ludcmp INTEGER i,its,j,indx(NP) REAL *8 d,den,f,fold,stpmax,sum,temp,test,fjac(NP,NP),g(NP),p(NP), *xold(NP),fmin EXTERNAL fmin

nn=n f=fmin(x) test=0.

do 11 i=1,n if (abs(fvec(i)) .gt. test)test=abs(fvec(i))

11

continue if (test .lt. .01*TOLF) then check= .false.

Made with FlippingBook Publishing Software