code: def load_ns_sequence(eos_name):
      
      ns_sequence = []
  
      if eos_name == '2H':
          ns_sequence_path = os.path.join(pycbc.tmpltbank.NS_SEQUENCE_FILE_DIRECTORY, 'equil_2H.dat')
          #ns_sequence_path = os.path.join(NS_SEQUENCE_FILE_DIRECTORY, 'equil_2H.dat')
          ns_sequence = np.loadtxt(ns_sequence_path)
      else:
          print('Only the 2H EOS is currently supported!')
          print('If you plan to use a different NS EOS, be sure not to filter')
          print('too many templates!\n')
          raise Exception('Unsupported EOS!')
  
      max_ns_g_mass = max(ns_sequence[:,0])
  
      return (ns_sequence, max_ns_g_mass)
docs: Load the data of an NS non-rotating equilibrium sequence
      generated using the equation of state (EOS) chosen by the
      user.  [Only the 2H 2-piecewise polytropic EOS is currently
      supported.  This yields NSs with large radiss (15-16km).]
  
      Parameters
      -----------
      eos_name: string
          NS equation of state label ('2H' is the only supported
          choice at the moment)
  
      Returns
      ----------
      ns_sequence: 3D-array
          contains the sequence data in the form NS gravitational
           mass (in solar masses), NS baryonic mass (in solar
           masses), NS compactness (dimensionless)
      max_ns_g_mass: float
          the maximum NS gravitational mass (in solar masses) in
          the sequence (this is the mass of the most massive stable
          NS)