やっと

非平衡の化学反応とガスの冷却/加熱が解けるようになった. 密度の低いところの reionization 時の振る舞いが厳しくて, そこをちゃんと解くには, implicit に化学反応を解いた後に, implicit に冷却/加熱を解くのではなく, 化学反応 -> 冷却/加熱 -> 化学反応 -> 冷却/加熱 -> ... と 1 つの timestep の中をこいつらの timescale で subcycling しながら解かないとひどいことになってしまう, ということを思い知ったのだった.

多分これで simulation に使える状態になったはず.