Make station dimension for hisfiles less hard-coded

[veenstrajelmer]

Recently (from DIMRset 2.27.03 onwards, 30 may 2014), the hisfile was refactored and the stations dimension was renamed to station.

That dimension is hardcoded in interpolate_grid2bnd.py:

dfm_tools/dfm_tools/interpolate_grid2bnd.py

Line 517 in 03cabed

hisstations_pd = data_xr_his.stations.to_dataframe() #TODO: add check if stations are strings? (whether preprocess_hisnc was used)

And also in get_nc.py:

dfm_tools/dfm_tools/get_nc.py

Line 538 in 03cabed

raise ValueError(f'ERROR: unexpected number of dimensions in requested squeezed variable ({data_xr_sel[varname].shape}), first use data_xr.isel(stations=int) to select a single station') #TODO: can also have a different cause, improve message/testing?

Try to derive the station dimension automatically just like dfmt.preprocess_hisnc() and prevent hard-coding.