urdme_model = _translate(self.model, param_values, dif0, y0, initial_dist)
urdme_model.timespan(tspan)
urdme_model.mesh = mesh
result = pyurdme.run(number_of_trajectories = n_runs, solver = solver, seed = seed )// output time points (in case they aren"t the same tspan, which is possible in BNG)
self.tout = result.get_timespan()
// species