title check that grains equilibriate at correct temp in ste limit
no buffering
c this also tests that the code works when H He ionization is ZERO
black 3 lte
init file "hheonly.ini"
stop zone 1
grain abund 1 type 1
grain abund 1 type 2
grain abund 1 type 3
grain abund 1 type 4
grain abund 1 type 5
grain abund 1 type 6
grain abund 1 type 7
grain abund 1 type 8
assert grain temperature 1 3 error 0.005
assert grain temperature 2 3 error 0.005
assert grain temperature 3 3 error 0.005
assert grain temperature 4 3 error 0.005
assert grain temperature 5 3 error 0.005
assert grain temperature 6 3 error 0.005
assert grain temperature 7 3 error 0.005
assert grain temperature 8 3 error 0.005
punch asserts last "grainlte.asr"
hden 5
eden 5
iterate
print last iteration
c the equilibrium temp would not be correct since current H_2 molecule
c does not go to ste, being an equivalent two-level system
constant temperature 3
c grainlte.in