Error during grid.show_dashboard

[itroitskaya]

Hi Gully, I am a PhD student working with @benjaminpope the The University of Queensland. I am trying to use interactive fitting with PHOENIX dashboard and get this error during grid.show_dashboard. Could you please have a look?

from gollum.phoenix import PHOENIXGrid
from specutils import Spectrum1D
import pandas as pd
import astropy.units as u
import numpy as np
from readmultispec import readmultispec

from IPython.display import HTML
from IPython.display import Image

data = readmultispec('./98Tau_2023-02-10_04h48m13s_cb.spec.fits')
data_path = '~/Dev/HiResFITS/PHOENIX-ACES-AGSS-COND-2011'

bdss_spectrum = Spectrum1D(spectral_axis=data['wavelen'][0]*u.micron,
                           flux=data['flux'][0]*u.ct)

wl_lo, wl_hi = (bdss_spectrum.wavelength.value.min(),
                bdss_spectrum.wavelength.value.max())

grid = PHOENIXGrid(wl_lo=wl_lo, wl_hi=wl_hi, path=data_path)

grid.show_dashboard(data=bdss_spectrum)

ERROR:tornado.application:Uncaught exception GET /autoload.js?bokeh-autoload-element=p1002&bokeh-absolute-url=http://localhost:54194&resources=none (::1)
HTTPServerRequest(protocol='http', host='localhost:54194', method='GET', uri='/autoload.js?bokeh-autoload-element=p1002&bokeh-absolute-url=http://localhost:54194&resources=none', version='HTTP/1.1', remote_ip='::1')
Traceback (most recent call last):
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/tornado/web.py", line 1786, in _execute
    result = await result
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/bokeh/server/views/autoload_js_handler.py", line 62, in get
    session = await self.get_session()
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/bokeh/server/views/session_handler.py", line 145, in get_session
    session = await self.application_context.create_session_if_needed(session_id, self.request, token)
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/bokeh/server/contexts.py", line 242, in create_session_if_needed
    self._application.initialize_document(doc)
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/bokeh/application/application.py", line 192, in initialize_document
    h.modify_document(doc)
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/bokeh/application/handlers/function.py", line 142, in modify_document
    self._func(doc)
  File "/Users/uqitroit/Dev/gollum/src/gollum/phoenix.py", line 302, in create_interact_ui
    scalar_norm = np.percentile(self[0].flux.value, 95)
  File "<__array_function__ internals>", line 200, in percentile
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/numpy/lib/function_base.py", line 4205, in percentile
    return _quantile_unchecked(
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/numpy/lib/function_base.py", line 4473, in _quantile_unchecked
    return _ureduce(a,
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/numpy/lib/function_base.py", line 3752, in _ureduce
    r = func(a, **kwargs)
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/numpy/lib/function_base.py", line 4639, in _quantile_ureduce_func
    result = _quantile(arr,
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/numpy/lib/function_base.py", line 4745, in _quantile
    take(arr, indices=-1, axis=DATA_AXIS)
  File "<__array_function__ internals>", line 200, in take
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/numpy/core/fromnumeric.py", line 190, in take
    return _wrapfunc(a, 'take', indices, axis=axis, out=out, mode=mode)
  File "/Users/uqitroit/anaconda3/envs/gollum_dev/lib/python3.8/site-packages/numpy/core/fromnumeric.py", line 57, in _wrapfunc
    return bound(*args, **kwds)
IndexError: cannot do a non-empty take from an empty axes.