Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Problem with gateLowerEdgeLevel of orifice #45

Open
rhutten opened this issue Nov 24, 2022 · 2 comments
Open

Problem with gateLowerEdgeLevel of orifice #45

rhutten opened this issue Nov 24, 2022 · 2 comments
Assignees

Comments

@rhutten
Copy link
Collaborator

rhutten commented Nov 24, 2022

Problem
The gateLowerEdgeLevel of an orifice is the position of the gate door's lower edge [m AD]. This is an absolute value in D-Flow FM.

However, the gateLowerEdgeLevels of orifices resulting from the D-HyDAMO converter is a relative value with respect to the crestLevel (laagstedoorstroomhoogte). This results from the 'hoogteopening' [m] being assigned as gateLowerEdgeLevel. More information about the D-HyDAMO onderspuier (orifice) can be found here.

The code below shows how the hoogteopening is assigned as gateLowerEdgeLevel:

  1. generate_weirs function in delft3dfmpy/converters/hydamo_to_dflowfm.py: image
  2. weirs_from_hydamo function in delft3dfmpy/io/dfmreader.py:
    image

Proposed solution

Adjust line 158 in generate_weirs function in delft3dfmpy/converters/hydamo_to_dflowfm.py to
orifices_dfm.at[idx, 'schuifhoogte'] = str(weir_opening.laagstedoorstroomhoogte + weir_mandev.hoogteopening)

@rhutten rhutten assigned rhutten and RuudHurkmans and unassigned rhutten Nov 24, 2022
@rhutten
Copy link
Collaborator Author

rhutten commented Nov 24, 2022

@RuudHurkmans, I checked the code in more detail. See above my proposed solution.

@RuudHurkmans
Copy link
Collaborator

This is correct, 'hoogteopening' was treated as an absolute level, but is in the HyDAMO format relative to the crest. In the code on github, the above solution has been implemented and it will be in the next release of the installable version.

As a workaround, in the meantime, an orifice can be added in the notebook as follows:

image

Here, an absolute value (m+NAP) should be provided, such as D-HYDRO expects it.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants