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