sum_X_s in des_reaction_model.f computed as zero for des phases with SPECIES_EQ == .false. when mmax > 1 and ANY_SOLIDS_SPECIES_EQ == .true. (at least one phase has SPECIES_EQ == .true.)
Fix: Change Line 110 to:
WHERE(PARTICLE_STATE(:MAX_PIP) == NORMAL_PARTICLE .AND. SPECIES_EQ(PIJK(:MAX_PIP,5))) &