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 = 5.8d-8
Q = 4500d+00
l = 0.5d-3
Hs = 0.03d+00
Eb1=42000.0d+00
Eb2=47000.0d+00
Eb3=35000.0d+00
Nl= 2.01d29
Nt1=3.6d25
Nt2=3.6d25
Nt3=0.0d25
/

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 = 21600.0d+00
ag_time = 900.0d+00
des_time = 180.0d+00
dt = 0.2d+00
T = 300.0d+00
Tend = 573.0d+00
dTdt = 0.013888888889d+00
/

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 = 50
/
