search

gwmod / nlmod

Python package to build, run and visualize MODFLOW 6 groundwater models in the Netherlands.
https://nlmod.readthedocs.io
MIT License
30 stars 2 forks source link

calculate_resistance calculates resistance between layers #310

Open rubencalje opened 8 months ago

rubencalje commented 8 months ago

Right now calculate_resistance calculates the resistance between layers (also see issue #134), and not the resistance of layers themselves (if that is correct english...).

I propose to rename the current calculate_resistance to calculate_resistance_between_layers, and create a new calculate_reistance with the formula:

c = nlmod.layers.calculate_thickness(ds) / ds['kv'].

Any objections?

rubencalje commented 8 months ago

Or add an optional argument to calculate_resistance named between_layers, which is set to False by default (which would change the current outcome of calculate_resistance). I prefer this option.

OnnoEbbens commented 8 months ago

Or add an optional argument to calculate_resistance named between_layers, which is set to False by default (which would change the current outcome of calculate_resistance). I prefer this option.

I also prefer this option