In solve_rate_cool.F rate k22 is used to determine a time step even if MultiSpecies is set to 1 (just the 6 species solver). This caused NaN b/c the rate was (correctly) not initialized for MultiSpecies=1. One line change.
Yup, great catch -- my fault, no doubt. Thank you!