docs: def ldcoeffs(teff,logg=4.5,feh=0):
      """
      Returns limb-darkening coefficients in Kepler band.
      """
      teffs = np.atleast_1d(teff)
      loggs = np.atleast_1d(logg)
  
      Tmin,Tmax = (LDPOINTS[:,0].min(),LDPOINTS[:,0].max())
      gmin,gmax = (LDPOINTS[:,1].min(),LDPOINTS[:,1].max())
  
      teffs[(teffs < Tmin)] = Tmin + 1
      teffs[(teffs > Tmax)] = Tmax - 1
      loggs[(loggs < gmin)] = gmin + 0.01
      loggs[(loggs > gmax)] = gmax - 0.01
  
      u1,u2 = (U1FN(teffs,loggs),U2FN(teffs,loggs))
      return u1,u2