docs: def read_simulation_temps(pathname,NumTemps):
      """Reads in the various temperatures from each TEMP#/simul.output file by knowing
          beforehand the total number of temperatures (parameter at top)
      """
  
      print("--Reading temperatures from %s/..." % pathname)
  
      # Initialize return variable
      temps_from_file = numpy.zeros(NumTemps, numpy.float64)
  
      for k in range(NumTemps):
          infile = open(os.path.join(pathname,'TEMP'+ str(k), 'simul'+str(k)+'.output'), 'r')
          lines = infile.readlines()
          infile.close()
  
          for line in lines:
              if (line[0:11] == 'Temperature'):
                  vals = line.split(':')
                  break
          temps_from_file[k] = float(vals[1])
  
      return temps_from_file