calculate transmissivity (kD) and resistance (c)

[tomvansteijn]

The package already contains a useful function for computing layer thicknesses. Of course it would be nice to have an additional function for calculating transmissivity and resistance as well.

[rubencalje]

Hi Tom, good idea. The formulas are quite simple:

kD = nlmod.layers.calculate_thickness(ds) * ds['kh']

and

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

[OnnoEbbens]

The formula @rubencalje describes calculates the resistance within a modellayer. I added the function calculate_resistance to calculate the resistance between model layers.

This function can be called using:
nlmod.dims.calculate_resistance(ds, kv='kv', thickness='thickness')

[OnnoEbbens]

Added both:

• resistance: nlmod.dims.calculate_resistance(ds)
• transmissivity: nlmod.dims.calculate_transmissivity(ds)