table powerlaw
save heating "heating.het"
hden 3
stop zone 1
constant temperature 1e4 K
ionization parameter -2
