From 2a47693f0eccb34f31a90ab077a64df33de3e14c Mon Sep 17 00:00:00 2001 From: Simon Date: Fri, 24 Jul 2020 15:09:47 +0200 Subject: [PATCH] Put velocity and magnetic field into SI units in Ek_Eb function --- postprocessor.py | 20 +++++++++++--------- 1 file changed, 11 insertions(+), 9 deletions(-) diff --git a/postprocessor.py b/postprocessor.py index af926e5..2730a0c 100644 --- a/postprocessor.py +++ b/postprocessor.py @@ -314,8 +314,6 @@ class PostProcessor(HDF5Container): else: rho_bins = np.linspace(np.min(rho), np.max(rho), bins) - weights = mass_func(self.cells) - # For each cell, bin_number contains the number of the bins it belongs to bin_number = np.zeros(len(B)) @@ -344,13 +342,19 @@ class PostProcessor(HDF5Container): """ self.load_cells() mean_speed = self.save.get_node("/globals/mwa_speed").read() - vel_fluct = (self.cells)["vel"] - mean_speed + mean_speed = mean_speed * self.info["unit_velocity"].express(cst.km_s) + vel_fluct = (self.cells)["vel"] * self.info["unit_velocity"].express( + cst.km_s + ) - mean_speed B_norm = getter_B_int(self.cells) - v_norm = np.sqrt(np.sum((vel_fluct * 10 ** (3)) ** 2, axis=1)) # v [km/s] - print(v_norm) + B_norm = B_norm * self.info["unit_mag"].express(cst.T) + v_norm = np.sqrt( + np.sum((vel_fluct * 10 ** (3)) ** 2, axis=1) + ) # v_norm [m/s] et vel_fluct [km/s] rho = getter_rho(self.cells) + rho_kg_m3 = rho * self.info["unit_density"].express(cst.kg_m3) eb = 0.5 * (B_norm) ** 2 / (4 * np.pi * 10 ** (-7)) # mettre le bon mu - ek = 0.5 * v_norm ** 2 * rho # mettre la masse de la cellule + ek = 0.5 * v_norm ** 2 * rho_kg_m3 rapport = ek / eb if logbins: @@ -360,8 +364,6 @@ class PostProcessor(HDF5Container): else: rho_bins = np.linspace(np.min(rho), np.max(rho), bins) - weights = mass_func(self.cells) - # For each cell, bin_number contains the number of the bins it belongs to bin_number = np.zeros(len(B_norm)) @@ -385,7 +387,7 @@ class PostProcessor(HDF5Container): return ({"rho": centers, "Ek_Eb_rho": ek_eb}, {"logbins": logbins}) def cos_vfluct_B(self): - + "return the cos of the angle between the magnetic field and the velocity fluctuation field" mean_speed = self.save.get_node("/globals/mwa_speed").read() def getter_cos_vfluct_B(dset):