From 12c25a0742e37b5e95da6f33531809b1e18772f2 Mon Sep 17 00:00:00 2001 From: rstiskalek Date: Wed, 3 Jul 2024 09:48:13 +0100 Subject: [PATCH] Update to linear interpolation --- scripts/field_los.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/scripts/field_los.py b/scripts/field_los.py index 1d393a6..6d734f8 100644 --- a/scripts/field_los.py +++ b/scripts/field_los.py @@ -297,7 +297,8 @@ def interpolate_field(pos, simname, nsim, MAS, grid, dump_folder, rmax, density = get_field(simname, nsim, "density", MAS, grid) rdist, finterp = csiborgtools.field.evaluate_los( density, sky_pos=pos, boxsize=boxsize, rmax=rmax, dr=dr, - smooth_scales=smooth_scales, verbose=verbose) + smooth_scales=smooth_scales, verbose=verbose, + interpolation_method="linear") print(f"Writing temporary file `{fname_out}`.") with File(fname_out, 'w') as f: @@ -314,7 +315,8 @@ def interpolate_field(pos, simname, nsim, MAS, grid, dump_folder, rmax, rdist, finterp = csiborgtools.field.evaluate_los( velocity[0], velocity[1], velocity[2], sky_pos=pos, boxsize=boxsize, rmax=rmax, dr=dr, - smooth_scales=smooth_scales, verbose=verbose) + smooth_scales=smooth_scales, verbose=verbose, + interpolation_method="linear") with File(fname_out, 'a') as f: f.create_dataset("velocity", data=finterp)