Merge DFSUs with different mesh files but same area

[RenewAndPlay]

I try to merge 2 DFSU files. Essentially they cover the same area, although 1 of the files cover a greater area. In the overlapping area the elements are of completely same size. However elements in both files will have different element_ids.

I want to merge the 2 dfsu files with all items from both of them into a new dfsu file.

1 DFSU file covers a greater area and consists of a time series for current speed, current direction

1 DFSU file covers a smaller area (which can also be found in above listed file) but only contains a single time step with items of extreme values for RP 1, 10, 25 and 50 year for current speeds.

The merged DFSU file should contain the cropped DFSU file with time series of current speed and direction and for every time step items of extreme values for RP 1, 10, 25 and 50 year for current speeds should be added- they will be constant throughout the time series.

I have tried the following without any results so far, since I don't know how to add a DFSU with time series to another mesh, even though they cover the same area. I don't want to interpolate, but instead I want to select discrete values of the time series file.

Can you help me lead me on the way in the task combining 2 DFSUs?

Br Andreas

#%% Load dfsu files
folder = 'HD'
# normal_dfsu = 'HD_EINS_MDD_normalCS.dfsu'
extreme_dfsu = 'HD_EINS_MDD_extremeCS.dfsu'

main_dfsu = 'M21HD_NSEI_area_2010.dfsu'
script_dir = os.path.dirname(os.path.abspath(__file__))

main_dfsu = os.path.join(script_dir, folder, main_dfsu)
extreme_dfsu = os.path.join(script_dir, folder, extreme_dfsu)

main_dfsu_read = mikeio.read(main_dfsu)
extreme_dfsu_read = mikeio.read(extreme_dfsu)

#%% Mesh handling

Msh_Extreme = mikeio.Mesh(extreme_dfsu)
Msh_Main = mikeio.Mesh(main_dfsu)

Msh_Extreme.plot()
Msh_Main.plot()

# Find elements with matching coordinates of Msh_Extreme in Msh_Main
El_Extreme = Msh_Extreme.element_coordinates
El_Main = Msh_Main.element_coordinates

# Find elements with matching coordinates of Msh_Extreme in Msh_Main
matching_elements = []

for idx_extreme, coord_extreme in enumerate(El_Extreme):
    for idx_main, coord_main in enumerate(El_Main):
        if all(coord_extreme[0:2] == coord_main[0:2]):
            matching_elements.append((Msh_Extreme.element_ids[idx_extreme], Msh_Main.element_ids[idx_main]))
            break

#%% Plot only matching elements in Msh_Main

# Assuming matching_elements contains pairs of element IDs (Extreme, Main)
element_ids_main = [main_element_id for (extreme_element_id, main_element_id) in matching_elements]

# Plot elements in Msh_Main using the extracted element IDs
Msh_Main.plot(elements=element_ids_main)

# Save new dataset with only elements of "element_ids_main"
main_dfsu_reduc = main_dfsu_read.isel(element=element_ids_main)

#%% Append main_dfsu values to extreme_dfsu

time_new = main_dfsu_reduc.time
geometry_new = extreme_dfsu_read.geometry
# Make new dataset of "extreme_dfsu_read". This dataset only has 1 timestep, but instead all items in dataset should be copied to new dataset for each time step of "time"

items_comb = extreme_dfsu_read.items+main_dfsu_reduc.items

# Initialize a new Dataset based on mesh from extreme_dfsu_read
new_extreme = extreme_dfsu_read.copy()

# Create new dataarray 
time_new = np.array(time_new)

# Create a new DataArray with the correct dimensions
da = mikeio.DataArray(data=np.tile(extreme_dfsu_read[0], (len(time_new), 1)),
                      time=time_new,
                      item=mikeio.ItemInfo("Current speed", mikeio.EUMType.Current_Speed),
                      geometry=geometry_new,
                      dims=("time", "element"))



#mikeio.Dataset.concat([main_dfsu_reduc,extreme_dfsu_read])