c----------------------
c D0 = diffusion coefficient, m/s2
c Q = diffusion activation energy, J
c l = half thickness, m
c Hs = hydrogen concentration at the surface, ppm
c Eb = trap binding energy, J/mol
c Nl = lattice density, /m3
c Nt = trap density, /m3
c----------------------

&specimen
D0 = 4.79e-7
Q = 51590
l = 0.5e-3
Hs = 23
Eb1=6100.0
Eb2=40000.0
Eb3=40000.0
Nl= 9.6568e28
Nt1=1.2071e28
Nt2=0.0
Nt3=0.0
/

c-------------------
c mode0 = kinetic
c mode1 = local eq.
c ch_time = charging time, s
c ag_time = aging time, s
c des_time = desorption measuring time interval, s
c dt = increase time interval, s
c T = initial temperature, K
c Tend = target temperature, K
c dTdt = heating rate, K/s
c-------------------

&analysis
mode = 1
ch_time = 259200.0
ag_time = 1800.0
des_time = 180.0
dt = 0.2
T = 300.0
Tend = 673.0
dTdt = 0.0277777777778
/

c----------------------
c icharge1 = simulate charging / icharge0 = read 'charging.txt'
c iage1 = simulate aging / iage0 = read 'aging.txt'
c ishape0 = flat
c ishape1 = cylinder
c----------------------

&option
icharge = 1
iage = 1
ishape = 0
grid=500
/