def concentration(self, pM=False):
concentr = (self.num_particles / NA) / self.box.volume_L<EOL>if pM:<EOL><INDENT>concentr *= <NUM_LIT><EOL><DEDENT>return concentr<EOL>
Return the concentration (in Moles) of the particles in the box.