velcube extraction, bugfix

This commit is contained in:
Noe Brucy
2022-08-23 10:11:28 +02:00
parent cdafe5cb61
commit 5d2e61c657
+1 -1
View File
@@ -35,7 +35,7 @@ def write_data(filename, vel, dens):
def extract_from_pp(pp): def extract_from_pp(pp):
vel = get_velocity_cubes(pp, unit=U.km_s) vel = get_velocity_cubes(pp, unit=U.km_s)
dens = get_density_cube(pp, unit=U.H_cc) dens = get_density_cube(pp, unit=U.H_cc)
write_data(f"{pp.run}_{pp.num}_velrho.data", vel, dens) write_data(f"{pp.path_out}/{pp.run}_{pp.num}_velrho.data", vel, dens)
def extract(path, snap_number): def extract(path, snap_number):