code: def ns_g_mass_to_ns_b_mass(ns_g_mass, ns_sequence):
      
      x = ns_sequence[:,0]
      y = ns_sequence[:,1]
      f = scipy.interpolate.interp1d(x, y)
  
      return f(ns_g_mass)
docs: Determines the baryonic mass of an NS given its gravitational
      mass and an NS equilibrium sequence.
  
      Parameters
      -----------
      ns_g_mass: float
          NS gravitational mass (in solar masses)
      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)
  
      Returns
      ----------
      float
          The NS baryonic mass (in solar massesr**3*(r**2*(r-6)+chi**2*(3*r+4))+chi**4*(3*r*(r-2)+chi**2))