From ed9ec08f9da9e9fac908a0ee5931635be1dd8c07 Mon Sep 17 00:00:00 2001
From: Jean-Baptiste Bayle <j2b.bayle@gmail.com>
Date: Tue, 2 Feb 2021 13:25:56 +0100
Subject: [PATCH] Increase spline order for orbits
---
lisainstrument/instrument.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/lisainstrument/instrument.py b/lisainstrument/instrument.py
index c59a8f7..96095d5 100644
--- a/lisainstrument/instrument.py
+++ b/lisainstrument/instrument.py
@@ -251,11 +251,11 @@ class Instrument:
orbitf = h5py.File(self.orbit_file, 'r')
logging.debug("Interpolating proper pseudo-ranges")
self.pprs = ForEachMOSA(lambda mosa: scipy.interpolate.InterpolatedUnivariateSpline(
- orbitf['tcb']['t'][:], orbitf['tcb'][f'l_{mosa}']['ppr'], k=1, ext='raise')(self.physics_t)
+ orbitf['tcb']['t'][:], orbitf['tcb'][f'l_{mosa}']['ppr'], k=5, ext='raise')(self.physics_t)
)
logging.debug("Interpolating proper pseudo-range derivatives")
self.d_pprs = ForEachMOSA(lambda mosa: scipy.interpolate.InterpolatedUnivariateSpline(
- orbitf['tcb']['t'][:], orbitf['tcb'][f'l_{mosa}']['d_ppr'], k=1, ext='raise')(self.physics_t)
+ orbitf['tcb']['t'][:], orbitf['tcb'][f'l_{mosa}']['d_ppr'], k=5, ext='raise')(self.physics_t)
)
orbitf.close()
elif orbits is None:
--
GitLab