Aquila-Consortium/csiborgtools
0
0
Fork 0
mirror of https://github.com/Richard-Sti/csiborgtools.git synced 2025-04-19 06:40:54 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Update H0 values for the void

Browse source
This commit is contained in:
rstiskalek 2024-09-19 11:19:00 +01:00
parent 947ac009a1
commit b3f295122a
1 changed files with 11 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

17
scripts/field_prop/field_los_indranil_void.py
View file

@ -28,7 +28,13 @@ from tqdm import trange
def interpolate_indranil_void(kind, nsims, RA, dec, rmax, dr, dump_folder, def interpolate_indranil_void(kind, nsims, RA, dec, rmax, dr, dump_folder,
catalogue): catalogue):
if kind not in ["exp", "gauss", "mb"]: if kind == "mb":
h = 0.7615
elif kind == "gauss":
h = 0.7724
elif kind == "exp":
h = 0.7725
else:
raise ValueError(f"Unknown void kind: `{kind}`.") raise ValueError(f"Unknown void kind: `{kind}`.")
kind = kind.upper() kind = kind.upper()
@ -48,11 +54,10 @@ def interpolate_indranil_void(kind, nsims, RA, dec, rmax, dr, dump_folder,
fname = join(fdir, f"v_pec_{kind}profile_rLG_{nsim}.dat") fname = join(fdir, f"v_pec_{kind}profile_rLG_{nsim}.dat")
data = np.loadtxt(fname) data = np.loadtxt(fname)
# The grid is in Mpc # The grid is in Mpc from 0 to 251 but we convert to Mpc / h
r_grid = np.arange(0, 251) r_grid = np.arange(0, 251) * h
phi_grid = np.arange(0, 181) phi_grid = np.arange(0, 181)
# The input is in Mpc/h, so we need to convert to Mpc r_eval = np.arange(0, rmax, dr).astype(float)
r_eval = np.arange(0, rmax, dr).astype(float) / 0.674
model_axis = SkyCoord(l=117, b=4, frame='galactic', unit='deg').icrs model_axis = SkyCoord(l=117, b=4, frame='galactic', unit='deg').icrs
coords = SkyCoord(ra=RA, dec=dec, unit='deg').icrs coords = SkyCoord(ra=RA, dec=dec, unit='deg').icrs
@ -76,7 +81,7 @@ def interpolate_indranil_void(kind, nsims, RA, dec, rmax, dr, dump_folder,
# Write the output, homogenous density. # Write the output, homogenous density.
density = np.ones_like(result) density = np.ones_like(result)
with File(fname_out, 'a') as f_out: with File(fname_out, 'a') as f_out:
f_out.create_dataset(f"rdist_{k}", data=r_eval * 0.674) f_out.create_dataset(f"rdist_{k}", data=r_eval)
f_out.create_dataset(f"density_{k}", data=density) f_out.create_dataset(f"density_{k}", data=density)
f_out.create_dataset(f"velocity_{k}", data=result) f_out.create_dataset(f"velocity_{k}", data=result)