watex.utils.geotools.get_s_thicknesses#

watex.utils.geotools.get_s_thicknesses(grouped_z, grouped_s, display_s=True, station=None)[source]#

Compute the thickness of each stratum from the grouped strata from the top to the bottom.

Parameters:

  • grouped_z – depth grouped according its TRES

  • grouped_s – strata grouped according to its TRES

  • s_display – bool, display infos in stdout

Returns:

  • thick : The thickness of each layers

  • strata: name of layers

  • status: check whether the total thickness is equal to the

    depth of investigation(doi). Iftrue: `coverall= 100% otherwise coverall is less which mean there is a missing layer which was not probably taking account.

Example:
>>> import watex.geology.core as GC
>>> obj= GC.quick_read_geomodel()
>>> sg, _, zg, _= fit_stratum_property (obj.fitted_rocks,
...                                    obj.z, obj.nmSites[:, 0]  )
>>> get_s_thicknesses( zg, sg)
... ([13.0, 16.0, 260.0, 240.0, 470.0],
...     ['*i', 'igneous rocks', 'granite', 'igneous rocks', 'granite'],
...     'coverall =100%')