An integration scheme for stiff solid–gas reactor models