Geometry

Created on Wed May 26 11:55:30 2021

@author: sam

nupyprop.geometry.PREMgramVSang(beta_deg, idepth)
nupyprop.geometry.cdtot(x_v, beta_deg, idepth)
Parameters:

  • x_v (None) – Lambda function (integration variable).

  • beta_deg (float) – Earth emergence angle (beta), in degrees.

  • idepth (int) – Depth of water layer in km.

Returns:

cdtot_val – Integrated column depth.

Return type:

float

nupyprop.geometry.columndepth(beta_deg, idepth)
Parameters:

  • beta_deg (float) – Earth emergence angle (beta), in degrees.

  • idepth (int) – Depth of water layer in km.

Returns:

columndepth_val – Column depth, in g/cm^2.

Return type:

float

nupyprop.geometry.create_traj_table(idepth)
Parameters:

idepth (int) – Depth of water layer in km.

Returns:

Adds trajectory table lookup entries in lookup_table.h5.

Return type:

None

nupyprop.geometry.find_interface(idepth)
Parameters:

idepth (int) – Depth of water layer in km.

Returns:

The (Earth emergence) angle at which the transition from water to rock occurs, in degrees.

Return type:

float

nupyprop.geometry.gen_col_trajs(idepth)
Parameters:

  • idepth (int) – Depth of water layer in km.

  • Returns'beta' ('Earth emergence angle, in degrees',) – ‘xalong’:’Distance in water, in km’, ‘cdalong’:’Column depth at xalong, in g/cm^2’})

  • ------- –

  • betad_fix (ndarray) – 1D array containing 13500 entries, with repeating (x100) entries from 0.1 deg to 90 deg.

  • xalong (ndarray) – 1D array containing distance in water, in km.

  • cdalong (where x_v limits go from 0 to) – 1D array containing column depth at xalong, in g/cm^2.

  • Essentially –

  • integral(cdtot(x_v (cdalong =) –

  • beta –

  • idepth)) –

  • cdalong –

  • value. (for each beta) –

nupyprop.geometry.gen_water_trajs(idepth)
Parameters:

idepth (int) – Depth of water layer in km.

Returns:

  • beta_arr (ndarray) – 1D array containing Earth emergence angles from 0.1 deg to 90 deg.

  • chord (ndarray) – 1D array containing chord length, in km.

  • water (ndarray) – 1D array containing final water layer distance, in km.

nupyprop.geometry.pathlength(tnadir)
Parameters:

tnadir (float) – Nadir angle, in radians.

Returns:

pathlength – Path length, in km.

Return type:

float

nupyprop.geometry.sagitta_deg(beta_deg)
Parameters:

beta_deg (float) – Earth emergence angle (beta), in degrees.

Returns:

sagitta – The sagitta (distance from the center of the arc to the center of its base), in km.

Return type:

float

nupyprop.geometry.trajlength(beta_deg)
Parameters:

beta_deg (float) – Earth emergence angle (beta), in degrees.

Returns:

trajlength – Trajectory length, in km.

Return type:

float