coincident_profile_functions module

coincident_profile_functions.compare_adjacent_profiles(profiles, coincident_profiles, approach)[source]

Calculate and compare bulk properties of overlapping/coincident velocity profiles.

This function takes velocity profiles and coincident profiles, and calculates and compares various properties between them.

Parameters:

  • profiles (list of dictionaries) – Velocity profiles containing location, velocity, density, etc. Output from run_convert_velocities.

  • coincident_profiles (list of lists) – Output from get_coincident_profiles. It’s a list in the same format as v_profiles with empty entries where no coincident profile exists and an entry where there is a profile that coincides with a v_profile with a corresponding index.

  • approach (str) – Method used to calculate velocities. Required since there may not be a calculated vp profile. It can be ‘brocher’ or ‘stephenson’.

Returns:

  • vp_vs_vpcalc (numpy.ndarray) – Comparison between vp, vs, and calculated vp. If calculated vp is not available, it writes np.nan.

  • rho_rho (numpy.ndarray) – Comparison of density estimates between vp and vs profiles.

  • moho_moho (numpy.ndarray) – Comparison of moho depth estimates between vp and vs profiles.

  • lon_lat (numpy.ndarray) – Locations of coincident profiles.

coincident_profile_functions.get_coincident_profiles(profiles1, profiles2, buff)[source]

Find all surveys that are in the same place.

Parameters:

  • profiles1 (list) – List of vp profiles.

  • profiles2 (list) – List of vs profiles.

  • buff (float) – Buffer distance in km (i.e. search distance between different profiles to look for coincident profiles).

Returns:

  • coincident_profiles_all (list) – List of profiles where Vp and Vs are coincident. If no coincident profiles are available then None is written as the list entry.

  • profiles1 (list) – List of vp profiles.

  • profiles2 (list) – List of vs profiles.