[bugfix] do not call calcsourcespec if incidence angle is outside bounds (for whatever reason)

This commit is contained in:
Marcel Paffrath 2023-08-24 11:27:30 +02:00
parent ba4e6cfe50
commit 2e85d083a3

View File

@ -418,6 +418,10 @@ class MomentMagnitude(Magnitude):
distance = degrees2kilometers(a.distance)
azimuth = a.azimuth
incidence = a.takeoff_angle
if not 0. <= incidence <= 360.:
if self.verbose:
print(f'WARNING: Incidence angle outside bounds - {incidence}')
return
w0, fc = calcsourcespec(scopy, onset, self.p_velocity, distance,
azimuth, incidence, self.p_attenuation,
self.plot_flag, self.verbose)