High Order Finite Difference Methods with Subcell Resolution for Stiff Multispecies Discontinuity Capturing. Commun. Comput. Phys. 2020;17(2):317-336. doi:10.4208/cicp.250214.130814a