code: def print_sizes(self):
          
          float_size = 4
          MB = 1024 * 1024
          size_ = self.n_samples * float_size
          em_size = size_ * self.num_particles / MB
          pos_size = 3 * size_ * self.num_particles / MB
          print("  Number of particles:", self.num_particles)
          print("  Number of time steps:", self.n_samples)
          print("  Emission array - 1 particle (float32): %.1f MB" % (size_ / MB))
          print("  Emission array (float32): %.1f MB" % em_size)
          print("  Position array (float32): %.1f MB " % pos_size)
docs: Print on-disk array sizes required for current set of parameters.