diff --git a/kinetics.cpp b/kinetics.cpp index 2f6f6592..8b427856 100644 --- a/kinetics.cpp +++ b/kinetics.cpp @@ -1880,13 +1880,13 @@ set_reaction(int i, int use_mix, int use_kinetics) /* * Find surface */ - if (use.Get_surface_in() && use.Get_kinetics_in() && use.Get_kinetics_ptr() && !use.Get_kinetics_ptr()->Get_use_cvode() && reaction_step > 1) - { - // use.Set_surface_ptr(Utilities::Rxn_find(Rxn_surface_map, i)); - // appt: we may come here with zero kinetic reaction, but surface may have to keep DONNAN_DL - } - else - dl_type_x = cxxSurface::NO_DL; + //if (use.Get_surface_in() && use.Get_kinetics_in() && use.Get_kinetics_ptr() && !use.Get_kinetics_ptr()->Get_use_cvode() && reaction_step > 1) + //{ + // // use.Set_surface_ptr(Utilities::Rxn_find(Rxn_surface_map, i)); + // // appt: we may come here with zero kinetic reaction, but surface may have to keep DONNAN_DL + //} + //else + // dl_type_x = cxxSurface::NO_DL; if (use.Get_surface_in() == TRUE) { use.Set_surface_ptr(Utilities::Rxn_find(Rxn_surface_map, i)); diff --git a/mainsubs.cpp b/mainsubs.cpp index 27cb9ad2..ba317af0 100644 --- a/mainsubs.cpp +++ b/mainsubs.cpp @@ -1444,7 +1444,8 @@ xsurface_save(int n_user) temp_surface.Set_n_user(n_user); temp_surface.Set_n_user_end(n_user); temp_surface.Set_new_def(false); - temp_surface.Set_dl_type(dl_type_x); + //temp_surface.Set_dl_type(dl_type_x); + temp_surface.Set_dl_type(use.Get_surface_ptr()->Get_dl_type()); temp_surface.Set_solution_equilibria(false); temp_surface.Set_n_solution(-999);