Effects of process-generated hydrogen on RPV walls

PROGRAM PWR_CHEM_RK

C#################################################################### C# Program calculates concentrations of species # C# as a function of time for transient processes # C# in Zircaloy crevices in PWRs # C# using a Runge-Kutta method # C# # C# Theory and code developed by Dr. G. Engelhardt, # C# and Dr. D. D. Macdonald # C# Version 1 # C# Date: December 12, 1997 # C####################################################################

IMPLICIT REAL *8(A-H,O-Z) PARAMETER (KMAXX=200,NMAX=50,JMAX = 19 ) COMMON /path/ kmax,kount,dxsav,xp(KMAXX),yp(NMAX,KMAXX) 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 COMMON /RAD/ GAMn,GAMg,GAMa,factor COMMON /LP/ LP(6) DIMENSION Z(25) DIMENSION Y(20),VV(25),QQ(25) LOGICAL check 94 FORMAT (8E12.3) 95 FORMAT (8I5)

OPEN (UNIT = 4,FILE='out.txt',STATUS='unknown') OPEN (UNIT = 3,FILE='FIGURE.dat',STATUS='unknown') OPEN (UNIT = 5,FILE='stat.txt',STATUS='unknown')

WRITE (*,*) '*****************************************************' WRITE (*,*) WRITE (*,*) ' PROGRAM PWR-CHEM-RK ' WRITE (*,*) WRITE (*,*) '*****************************************************' WRITE (*,*)

WRITE (4,*) '*****************************************************' WRITE (4,*) WRITE (4,*) ' PROGRAM PWR-CHEM-RK' WRITE (4,*) WRITE (4,*) '*****************************************************' WRITE (4,*)

cc__________________________ INPUT BLOCK ________________________________