Hoeijmakers / tayph

Analysis of high resolution spectroscopic time-series of exoplanets
9 stars 8 forks source link

DS selection results in out of bounds #137

Open bprinoth opened 3 months ago

bprinoth commented 3 months ago

For partial transits only/no baseline: selection too close to the edge (upper or lower) in phase results in interpolation range out of bounds. Instead of killing the run, ask to reselect points instead.

Selected points:
   [38.806451612903174, 29.058823529411768]
   [41.50134408602145, 12.367647058823527]
Traceback (most recent call last):
  File "tayph_start.py", line 9, in <module>
    run.start_run('2023.dat', parallel=False, xcor_parallel=False)
  File "/usr/local/src/tayph/tayph/run.py", line 206, in start_run
    run_instance(params,parallel=parallel,xcor_parallel=xcor_parallel)
  File "/usr/local/src/tayph/tayph/run.py", line 977, in run_instance
    shadow.construct_doppler_model(rv,ccf_nn,dp,shadowname,xrange=[-200,200],Nxticks=20.0,
  File "/usr/local/src/tayph/tayph/shadow.py", line 600, in construct_doppler_model
    model_callback = fit_doppler_model(fig,ax,rv_sel,z,primer,dp,shadowname)
  File "/usr/local/src/tayph/tayph/shadow.py", line 301, in __init__
    p1 = float(p_i(primer[0][1]))
  File "/home/bibi/.local/lib/python3.8/site-packages/scipy/interpolate/_polyint.py", line 80, in __call__
    y = self._evaluate(x)
  File "/home/bibi/.local/lib/python3.8/site-packages/scipy/interpolate/_interpolate.py", line 752, in _evaluate
    below_bounds, above_bounds = self._check_bounds(x_new)
  File "/home/bibi/.local/lib/python3.8/site-packages/scipy/interpolate/_interpolate.py", line 786, in _check_bounds
    raise ValueError("A value ({}) in x_new is above "
ValueError: A value (29.058823529411768) in x_new is above the interpolation range's maximum value (29.0).