obspy.taup.slowness_model.SlownessModel.layer_number_above

SlownessModel.layer_number_above(depth, is_p_wave)[source]

Find the index of the slowness layer that contains the given depth.

Note that if the depth is a layer boundary, it returns the shallower of the two or possibly more (since total reflections are zero thickness layers) layers.

See also

layer_number_below()

Parameters:

  • depth (float or ndarray) – The depth to find, in km.

  • is_p_wave (bool) – Whether to look at P (True) velocity or S (False) velocity.

Returns:

The slowness layer containing the requested depth.

Return type:

int or ndarray (dtype = int, shape = depth.shape)

Raises:

SlownessModelError – If no layer in the slowness model contains the given depth.