diff --git a/Cooling/05_POWER_estimation_water_cooling.py b/Cooling/05_POWER_estimation_water_cooling.py
index e5ce273..8a707bb 100644
--- a/Cooling/05_POWER_estimation_water_cooling.py
+++ b/Cooling/05_POWER_estimation_water_cooling.py
@@ -19,9 +19,9 @@ def Q_heat(flow,d_T):
 def main():
 
 
-    d_T = 1.5
+    d_T = 2
 
-    flow = 0.3/5#/5 #m/s
+    flow = 1#/5 #m/s
     #flow *=2
     print(f"flow = {flow}m/s")
     V_t = 4.8 * 3.2 * 1e-6 * flow
diff --git a/Data_Coils/Data.py b/Data_Coils/Data.py
new file mode 100644
index 0000000..8d9444e
--- /dev/null
+++ b/Data_Coils/Data.py
@@ -0,0 +1,64 @@
+# %%
+import matplotlib.pyplot as plt
+import RigolWFM.wfm as rigol
+
+# %%
+
+hhpion = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhpion.wfm'
+hhpioff = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhpioff.wfm'
+hhon = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhon.wfm'
+hhoff = 'C:/Users/Joschka/Desktop/Coil_Data/New/hhoff.wfm'
+
+scope = 'DS1104Z-S'
+
+# %%
+hhpion = rigol.Wfm.from_file(hhpion, scope)
+hhpion.plot()
+plt.show()
+# %%
+hhpioff = rigol.Wfm.from_file(hhpioff, scope)
+hhpioff.plot()
+plt.show()
+
+# %%
+hhoff = rigol.Wfm.from_file(hhoff, scope)
+hhoff.plot()
+plt.show()
+
+# %%
+hhon = rigol.Wfm.from_file(hhon, scope)
+hhon.plot()
+plt.show()
+
+
+# %%
+print(w.channels[0].times)
+print(w.channels[0].volts)
+
+# %%
+
+ahhpion = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhpion.wfm'
+ahhpioff = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhpioff.wfm'
+ahhon = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhon.wfm'
+ahhoff = 'C:/Users/Joschka/Desktop/Coil_Data/New/ahhoff.wfm'
+
+scope = 'DS1104Z-S'
+
+# %%
+ahhpion = rigol.Wfm.from_file(ahhpion, scope)
+ahhpion.plot()
+plt.show()
+# %%
+ahhpioff = rigol.Wfm.from_file(ahhpioff, scope)
+ahhpioff.plot()
+plt.show()
+
+# %%
+ahhoff = rigol.Wfm.from_file(ahhoff, scope)
+ahhoff.plot()
+plt.show()
+
+# %%
+ahhon = rigol.Wfm.from_file(ahhon, scope)
+ahhon.plot()
+plt.show()
\ No newline at end of file
diff --git a/Test_class.py b/Test_class.py
index de2ba2e..a13d985 100644
--- a/Test_class.py
+++ b/Test_class.py
@@ -10,177 +10,16 @@ from src import B_field_calculation as bf
 
 from src import coil_class as BC
 
-from IPython import get_ipython
-get_ipython().run_line_magic('matplotlib', 'qt')
-#get_ipython().run_line_magic('matplotlib', 'inline')
 
-#set up axis
-x_m = np.linspace(-0.05, 0.05, 101)
-z_m = np.linspace(-0.05, 0.05, 101)
+# %%
 
+res = BC.BCoil.resistivity_copper(20)
+res2 = BC.BCoil.resistivity_copper(50)
+l = 64 * 2*np.pi *48e-3
+A = (0.25e-3)**2 * np.pi
+R1 = res * l/A
+print(R1)
+R2 = res2 *l/A
+print(R2)
 
-z = z_m*1e3
-x = x_m*1e3 #for plotting in mm
-
-#Import Values from simulation
-
-################# My simulation #########################
-I = 5
-HH = 1
-d_coils = 44
-R_mid = 44
-R_inner = 44-3*1.7
-
-layers = 6
-windings = 2
-wire_width = 1.7
-wire_height = 2.6
-
-
-HH_Coil1 = BC.BCoil(HH, d_coils ,R_mid, layers, windings, wire_width, wire_height)
-#HH_Coil1.print_info()
-#B_z_sim, B_x_sim = HH_Coil1.B_field(5, x, z)
-
-#Bz, B_x = bf.B_multiple_raster(I,HH,R_inner,d_coils,layers,windings,wire_width, wire_height, x_m,z_m)
-
-#B_test = B_field_ideal_AHH(layers*windings,I,R_inner*1e-3,d_coils*1e-3,z_m)
-
-#B_x = np.concatenate((-np.flip(B_r),B_r[1:len(B_r)]))
-
-
-HH_Coil1.B_quick_plot(I, x, z)
-
-
-
-
-
-#Calculate gradients/curvature
-B_z_sim_grad = np.gradient(np.gradient(B_z_sim,z_m),z_m)/1e4
-B_x_sim_grad = np.gradient(B_x_sim,x_m)/100
-#B_z_grad = np.gradient(np.gradient(Bz,z_m),z_m)/1e4
-B_z_grad = np.gradient(B_z,z_m)/100
-
-B_z_sim_grad = np.gradient(B_z_grad,z_m)/100
-
-B_x_grad = np.gradient(B_x,x_m)/100
-
-
-#Calculate relative differences in permille
-rel_diff_Bz = (B_z-B_z_sim)/np.mean(B_z)
-#rel_diff_Bx = (B_x-B_x_sim)/np.mean(B_x)
-rel_diff_Bz_grad = (B_z_grad-B_z_sim_grad)/np.mean(B_z_grad)
-
-rel_diff_Bz_grad_mean = (B_z_grad-B_z_sim_grad)/np.mean(B_z_grad)
-#rel_diff_Bx_grad = (B_x_grad-B_x_sim_grad)/np.mean(B_x_grad)
-
-#Plotting
-plt.figure(1,figsize=(20,18))
-
-plt.rcParams.update({'font.size': 15})
-plt.suptitle("Helmholtz coil field Bz along z-axis, comparison of simulations", fontsize=30)
-
-
-#Field plot
-##########################
-plt.subplot(3,2,1)
-plt.plot(z,B_z,linestyle = "solid", label = r"$Bz$: Result via elliptic integrals")
-plt.plot(z,B_z_sim,linestyle = "dashdot", label = r"$B_{z, sim}$: Numerical Matlab simulation")
-plt.plot(z,(B_z-B_z_sim), label = r"$Bz - B_{z, sim}$")
-#plt.xlim(-0.01,0.01)
-plt.title("B-field" ,fontsize = 30)
-
-plt.ylabel(r"$Bz$ [G]")
-plt.xlabel("z-axis [mm]")
-plt.legend()
-
-
-#############################
-plt.subplot(3,2,3)
-plt.plot(z,(B_z-B_z_sim), label = r"$Bz - B_{z, sim}$")
-plt.ylabel("absolute deviation [G]")
-plt.xlabel("z-axis [mm]")
-plt.legend()
-
-#############################
-plt.subplot(3,2,5)
-plt.plot(z,1000*rel_diff_Bz, label = "$(Bz - B_{z, sim}) / Bz$")
-plt.ylabel("relative deviation [‰]")
-plt.xlabel("z-axis [mm]")
-plt.legend()
-
-
-######################Gradient plot############################
-
-################
-plt.subplot(3,2,2)
-plt.plot(z,B_z_grad,linestyle = "solid", label = r"$\nabla_z^2 Bz$: Result via elliptic integrals")
-plt.plot(z,B_z_sim_grad,linestyle = "dashdot", label = r"$\nabla_z^2 B_{z, sim}$: Numerical Matlab sim.")
-plt.plot(z,(B_z_grad-B_z_sim_grad), label = r"$\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}$")
-
-plt.ylabel(r"$\nabla_z^2 Bz [G/cm^2]$")
-plt.xlabel("z-axis [mm]")
-plt.title("Curvature of B-field",fontsize = 30)
-plt.legend(loc='lower right')
-
-
-#################
-
-plt.subplot(3,2,4)
-plt.plot(z,(B_z_grad-B_z_sim_grad), label = r"$\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}$")
-plt.ylabel(r"absolute deviation $[G/cm^2]$")
-plt.xlabel("z-axis [mm]")
-plt.legend()
-
-#####################
-plt.subplot(3,2,6)
-plt.plot(z,1000*rel_diff_Bz_grad, label = r"$(\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}) / \nabla_z^2 Bz$")
-#plt.ylim(-57,10)
-plt.ylabel("relative deviation [‰]")
-plt.xlabel("z-axis [mm]")
-plt.legend()
-
-
-plt.savefig("output/HH_benchmark_5A_6x2.pdf")
-plt.show()
-
-############### relative deviation with averaging by the mean not the individual value ########################################
-plt.figure(2)
-
-plt.plot(z,1000*rel_diff_Bz_grad_mean, label = r"$(\nabla_z^2 Bz - \nabla_z^2 B_{z, sim}) / mean(\nabla_z^2 Bz)$")
-#plt.ylim(-57,10)
-plt.ylabel("relative deviation [‰]")
-plt.xlabel("z-axis [mm]")
-plt.legend()
-plt.savefig("output/HH_benchmark_5A_6x2_rel_deviation_via_mean.pdf")
-plt.show()
-
-##################### x-Axis #########################################################
-
-plt.figure(3)
-
-plt.rcParams.update({'font.size': 15})
-plt.suptitle("Helmholtz coil field B_x along x-axis, comparison of simulations", fontsize=30)
-
-
-#Field plot
-##########################
-
-plt.plot(x,B_x,linestyle = "solid", label = r"$B_x$: Result via elliptic integrals")
-plt.plot(x,B_x_sim,linestyle = "dashdot", label = r"$B_{x, sim}$: Numerical Matlab simulation")
-plt.plot(x,(B_x-B_x_sim), label = r"$B_x - B_{x, sim}$")
-#plt.xlim(-0.01,0.01)
-plt.title("B-field" ,fontsize = 30)
-
-plt.ylabel(r"$B_x$ [G]")
-plt.xlabel("x-axis [mm]")
-plt.legend()
-
-
-
-
-#################
-
-
-plt.savefig("output/HH_benchmark_5A_6x2_x-axis.pdf")
-plt.show()
-
+print((R2-R1)/R1)
\ No newline at end of file
diff --git a/src/physical_constants.py b/src/physical_constants.py
index cd26612..e05dd99 100644
--- a/src/physical_constants.py
+++ b/src/physical_constants.py
@@ -26,7 +26,7 @@ f_626 = 478.839e12
 lambda_626 = 626.082e-9
 
 #specific constants
-rho_copper_20 = 1.68e-8
+rho_copper_20 = 1.72e-8
 density_copper = 8.940e3
 
 # water