table power law
c 
ionization parameter -2
constant temperature 1e4 K
hden 0
stop zone 1
set dr 0
save heating "heat.het"
save cooling "heat.col"
