Attention!
To reproduce the example case shown in the reference, you are recommended to go into Auto_chemistry.f90， modify the expression for YDOT(15) and YDOT(16) to 0.d0 in the subroutine Fex_Plasma.