diff --git a/analysis/power_spectrum.py b/analysis/power_spectrum.py
index 64dacb5..93bb8ee 100644
--- a/analysis/power_spectrum.py
+++ b/analysis/power_spectrum.py
@@ -15,10 +15,10 @@ def get_power_spectrum(field, kmin=kmin, kmax=kmax, Nk=Nk, G=None):
             field.L0,
             field.L1,
             field.L2,
-            field.N0,
-            field.N1,
-            field.N2,
-            k_modes=np.concat([PowerSpectrum(field.L0,field.L1,field.L2,field.N0,field.N1,field.N2,).FourierGrid.k_modes[:10],np.logspace(
+            int(field.N0),
+            int(field.N1),
+            int(field.N2),
+            k_modes=np.concat([PowerSpectrum(field.L0,field.L1,field.L2,int(field.N0),int(field.N1),int(field.N2),).FourierGrid.k_modes[:10],np.logspace(
                 np.log10(kmin),
                 np.log10(kmax),
                 Nk,
@@ -43,10 +43,10 @@ def get_cross_correlations(field_A, field_B, kmin=kmin, kmax=kmax, Nk=Nk, G=None
             field_A.L0,
             field_A.L1,
             field_A.L2,
-            field_A.N0,
-            field_A.N1,
-            field_A.N2,
-            k_modes=np.concat([PowerSpectrum(field_A.L0,field_A.L1,field_A.L2,field_A.N0,field_A.N1,field_A.N2,).FourierGrid.k_modes[:10],np.logspace(
+            int(field_A.N0),
+            int(field_A.N1),
+            int(field_A.N2),
+            k_modes=np.concat([PowerSpectrum(field_A.L0,field_A.L1,field_A.L2,int(field_A.N0),int(field_A.N1),int(field_A.N2),).FourierGrid.k_modes[:10],np.logspace(
                 np.log10(kmin),
                 np.log10(kmax),
                 Nk,