From 5d9c285c76a272d9aef80760f42222224bc7bf34 Mon Sep 17 00:00:00 2001
From: Jeremy Fee <jmfee@usgs.gov>
Date: Tue, 19 May 2020 13:53:42 -0600
Subject: [PATCH] Use declination mean as reference for null method declination
---
geomagio/residual/Calculation.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/geomagio/residual/Calculation.py b/geomagio/residual/Calculation.py
index 8f51883fd..4001c5f11 100644
--- a/geomagio/residual/Calculation.py
+++ b/geomagio/residual/Calculation.py
@@ -46,7 +46,7 @@ def calculate(reading: Reading, adjust_reference: bool = True) -> Reading:
azimuth=reading.azimuth,
h_baseline=absoluteH.baseline,
measurements=reading.measurements,
- reference=adjust_reference and reference or mean,
+ reference=adjust_reference and reference or None,
)
# calculate scale
scale_value = None
@@ -87,6 +87,7 @@ def calculate_D_absolute(
D Absolute
"""
mean = average_measurement(measurements, DECLINATION_TYPES)
+ reference = reference or mean
# average mark
average_mark = average_measurement(measurements, MARK_TYPES).angle
# adjust based on which is larger
--
GitLab