On this ellipsoidal Fermi surface model, the element of Fermi surface area is expressed as \(dS_F = d\theta d\phi \sin {\theta }\sqrt{A^2B^2\cos ^2{\theta }+C^2(B^2\cos ^2{\phi }+A^2\sin ^2{\phi })\sin ^2{\theta }}\)  [1]}. The ratio of anisotropic parameters \(A\) , \(B\) , and \(C\)  are assumed from the experimental values of the directional dependence of the averaged Fermi velocity or the Fermi surface structure. Here we consider three cases: (I) Isotropic FS: The soft X-ray ARPES measurements observed the isotropic 3D Fermi surface [2]}. (II) Ellipsoidal FS (Fermi Surface): From the temperature gradient of the upper critical field, \((A,B,C)=(E_F/\hbar )\) (1/11000 m/s, 1/5500 m/s, 1/9500 m/s) [3]}. (III) Open FS: The vacuum UV synchrotron ARPES measurements fitted the cylindrical light electron band with \((A,B,C)\) =(0.313 Å\(^{-1}\) , 1.103 Å\(^{-1}\) , 0.729 Å\(^{-1}\) ) and the isotropic heavy electron pocket with \((A,B,C)\) =(0.229 Å\(^{-1}\) , 0.229 Å\(^{-1}\) , 0.229 Å\(^{-1}\) ) [4]}. The open light electron band is calculated by integrating the ellipsoidal model over the Brillouin zone (\(|k_z|<0.45\)  Å\(^{-1}\) ). In the cases I and II, the superfluid density only depends on the ratio of \(A\) , \(B\) , and \(C\)  because of the normalization of \(T\)  (see Eq. REF ). We note that there are gaps in diagonal directions in measured Fermi surfaces. These gaps in the diagonal directions affect the integral over the Fermi surface, but in our model there is no any node near these directions. Thus here we consider this effect as small enough to be ignored.
