example_long_selafin.py

[garyghgao]

I was trying to run this example for TELEMAC, there are errors, the error message are below, has anyone have the same problem before or can help on this problem, Thanks.
@Traceback (most recent call last):
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\models\basemodel\environment.py", line 649, in get_environment
reader.get_variables_interpolated(
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\basereader\variables.py", line 889, in get_variables_interpolated
env, env_profiles = self.get_variables_interpolated_xy(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\basereader\variables.py", line 755, in get_variables_interpolated_xy
env, env_profiles = self.get_variables_interpolated(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\basereader\unstructured.py", line 61, in get_variables_interpolated
env = self.get_variables(variables, time, reader_x, reader_y, z)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\reader_telemac_selafin.py", line 230, in get_variables
pm = self.extractslf(frames, duration, self.altitude_ID, idx_3D)[0]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^^^
IndexError: index 0 is out of bounds for axis 0 with size 0
12:24:59 DEBUG opendrift.models.basemodel.environment:675: Traceback (most recent call last):
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\models\basemodel\environment.py", line 649, in get_environment
reader.get_variables_interpolated(
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\basereader\variables.py", line 889, in get_variables_interpolated
env, env_profiles = self.get_variables_interpolated_xy(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\basereader\variables.py", line 755, in get_variables_interpolated_xy
env, env_profiles = self.get_variables_interpolated(
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\basereader\unstructured.py", line 61, in get_variables_interpolated
env = self.get_variables(variables, time, reader_x, reader_y, z)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "C:\Users\G.Gao.conda\envs\opendrift\Lib\site-packages\opendrift\readers\reader_telemac_selafin.py", line 230, in get_variables
pm = self.extractslf(frames, duration, self.altitude_ID, idx_3D)[0]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^^^
IndexError: index 0 is out of bounds for axis 0 with size 0@

[knutfrode]

Hi,

I do not have much insight into how this specific reader works, as I am normally not working with TELEMAC/Selafin output files.
Thus I also have no idea why this example has stopped working.
Perhaps its creator @Boorhin has an idea?

However, as discussed in #1293 it might be a better option to rather read and convert these SLF-files outside of OpenDrift (i.e. not using the specific selafin reader), and rather pass an Xarray object to the generic unstructured reader.