Some molecules are dissapearing again, check boundary later. Maybe write print number of molecules each timestep vs number of molecules inside boundary