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