Skip to content
Snippets Groups Projects
Commit b496da41 authored by Cain, Payton David's avatar Cain, Payton David
Browse files

Fix comments, imports, and spacing. Add diagnostics to testing

parent 3759d870
No related branches found
No related tags found
No related merge requests found
...@@ -43,3 +43,30 @@ for filename in os.listdir("etc/residual"): ...@@ -43,3 +43,30 @@ for filename in os.listdir("etc/residual"):
err_msg="Baselines not within 4 decimals", err_msg="Baselines not within 4 decimals",
verbose=True, verbose=True,
) )
# gather original and resulting diagnostics
o_diagnostics = original.diagnostics
r_diagnostics = original.diagnostics
# test mean mark values
assert_almost_equal(
o_diagnostics.mean_mark,
r_diagnostics.mean_mark,
decimal=2,
err_msg="Baselines not within 4 decimals",
verbose=True,
)
# test magnetic azimuth values
assert_almost_equal(
o_diagnostics.magnetic_azimuth,
r_diagnostics.magnetic_azimuth,
decimal=2,
err_msg="Baselines not within 4 decimals",
verbose=True,
)
# test meridian values
assert_almost_equal(
o_diagnostics.meridian,
r_diagnostics.meridian,
decimal=2,
err_msg="Baselines not within 4 decimals",
verbose=True,
)
...@@ -8,14 +8,9 @@ class Diagnostics(BaseModel): ...@@ -8,14 +8,9 @@ class Diagnostics(BaseModel):
Attributes Attributes
---------- ----------
meridian: claculated from declination measurements meridian: calculated from declination measurements
mean_mark: average mark angles from measurements mean_mark: average mark angles from measurements
magnetic_azimuh: after adjustment magnetic_azimuh: after adjustment
mark_azimuth: from metadata
declination: from declination measurements
inclination: from inclination measurements
h_component: H baseline from calculations
z_component: Z baseline from calculations
""" """
meridian: float = None meridian: float = None
......
...@@ -6,11 +6,11 @@ from obspy.core import UTCDateTime ...@@ -6,11 +6,11 @@ from obspy.core import UTCDateTime
import openpyxl import openpyxl
from .Absolute import Absolute from .Absolute import Absolute
from .Calculation import DECLINATION_TYPES, MARK_TYPES, average_measurement
from .Diagnostics import Diagnostics
from .Measurement import Measurement from .Measurement import Measurement
from .MeasurementType import MeasurementType as mt from .MeasurementType import MeasurementType as mt
from .Reading import Reading from .Reading import Reading
from .Diagnostics import Diagnostics
from .Calculation import DECLINATION_TYPES, MARK_TYPES, average_measurement
from . import Angle from . import Angle
...@@ -414,7 +414,6 @@ class SpreadsheetAbsolutesFactory(object): ...@@ -414,7 +414,6 @@ class SpreadsheetAbsolutesFactory(object):
meridian = average_measurement(measurements, DECLINATION_TYPES).angle meridian = average_measurement(measurements, DECLINATION_TYPES).angle
magnetic_azimuth = mean_mark - meridian magnetic_azimuth = mean_mark - meridian
if meridian > 180: if meridian > 180:
magnetic_azimuth = mean_mark - (meridian - 90) magnetic_azimuth = mean_mark - (meridian - 90)
if mean_mark > 180: if mean_mark > 180:
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment