-
Notifications
You must be signed in to change notification settings - Fork 15
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
Create constant infiltration rates based on lulc and reclass_table #103
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks @roeldegoede for this PR!
I've made a few small changes (see also comments), let me know if you agree.
It would be good to also update the csv tables - can you still do that?
Otherwise good to be merged.
hydromt_sfincs/sfincs.py
Outdated
elif lulc is not None: | ||
# landuse/landcover should always be combined with mapping | ||
if reclass_table is None and isinstance(lulc, str): | ||
reclass_table = join(DATADIR, "qinf", f"{lulc}_mapping.csv") |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I changed this to point the the "data/lulc" folder. Should we also modify the tables in there to contain a qinf column?
hydromt_sfincs/sfincs.py
Outdated
if ( | ||
not os.path.isfile(reclass_table) | ||
and reclass_table not in self.data_catalog | ||
): | ||
raise IOError(f"Infiltration mapping file not found: {reclass_table}") |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This doesn't work if reclass_table is a pandas DataFrame
In addition to the existing functionality, we can now also set up constant infiltration rates based on landuse/landcover and a reclass table. Before we could only use spatiallly varying infiltration fields.