ararpy 0.2.2__py3-none-any.whl → 0.2.3__py3-none-any.whl

ararpy/__init__.py

@@ -16,7 +16,7 @@ from . import calc, smp, files, thermo, test
 """ Information """
 
 name = 'ararpy'
-version = '0.2.2'
+version = '0.2.3'
 __version__ = version
 full_version = version
 last_update = '2026-01-01'

ararpy/calc/corr.py

@@ -272,6 +272,7 @@ def Monte_Carlo_F(ar40m: Tuple[float, float], ar39m: Tuple[float, float], ar38m:
                   G40: Tuple[float, float], G39: Tuple[float, float], G38: Tuple[float, float],
                   G37: Tuple[float, float], G36: Tuple[float, float],
                   stand_time_year: float, JNFactor: Tuple[float, float],
+                  KCaFactor: Tuple[float, float], KClFactor: Tuple[float, float],
                   **options):
     """
 
@@ -359,6 +360,9 @@ def Monte_Carlo_F(ar40m: Tuple[float, float], ar39m: Tuple[float, float], ar38m:
     R38v39k = random_normal_relative(*R38v39k, size=Monte_Carlo_Size)
     R36v38clp = random_normal_relative(*R36v38clp, size=Monte_Carlo_Size)
 
+    KCa = random_normal_relative(*KCaFactor, size=Monte_Carlo_Size)
+    KCl = random_normal_relative(*KClFactor, size=Monte_Carlo_Size)
+
     JNFactor = random_normal_absolute(*JNFactor, size=Monte_Carlo_Size)
 
     def do_simulation():
@@ -380,31 +384,23 @@ def Monte_Carlo_F(ar40m: Tuple[float, float], ar39m: Tuple[float, float], ar38m:
             _ar39 = (ar39m[i] / G39[i] - ar39b[i] / (G39[i] if blank_gain_corr else 1)) * P39Mdf * P39Decay
             _ar40 = (ar40m[i] / G40[i] - ar40b[i] / (G40[i] if blank_gain_corr else 1)) * P40Mdf
 
+            _ar37ca = _ar37
+
             if force_to_zero:
-                _ar36 = max(_ar36, 0)
-                _ar37 = max(_ar37, 0)
-                _ar38 = max(_ar38, 0)
-                _ar39 = max(_ar39, 0)
-                _ar40 = max(_ar40, 0)
+                _ar37ca = max(_ar37, 0)
 
-            _ar37ca = _ar37
             _ar39ca = _ar37ca * R39v37ca[i]
             _ar38ca = _ar37ca * R38v37ca[i]
             _ar36ca = _ar37ca * R36v37ca[i]
 
+            _ar39k = _ar39 - _ar39ca
+
             if force_to_zero:
-                _ar39ca = min(_ar39ca, _ar39)
-                _ar38ca = min(_ar38ca, _ar38)
-                _ar36ca = min(_ar36ca, _ar36)
+                _ar39k = max(_ar39k, 0)
 
-            _ar39k = _ar39 - _ar39ca
             _ar38k = _ar39k * R38v39k[i]
             _ar40k = _ar39k * R40v39k[i]
 
-            if force_to_zero:
-                _ar38k = min(_ar38k, _ar38 - _ar38ca)
-                _ar40k = min(_ar40k, _ar40)
-
             _ar38res = _ar38 - _ar38k - _ar38ca
             _ar36res = _ar36 - _ar36ca
 
@@ -433,6 +429,8 @@ def Monte_Carlo_F(ar40m: Tuple[float, float], ar39m: Tuple[float, float], ar38m:
 
             factor = JNFactor[i]
 
+            _CaK = _ar37ca / _ar39k / KCa[i]
+
             i += 1
 
             yield _ar36a, 0, _ar36ca, _ar36cl, _ar37ca, _ar38cl, _ar38a, 0, _ar38k, _ar38ca, _ar39k, _ar39ca, _ar40r, _ar40a, 0, _ar40k, \
@@ -442,7 +440,7 @@ def Monte_Carlo_F(ar40m: Tuple[float, float], ar39m: Tuple[float, float], ar38m:
                   (_ar39k * factor) / _ar40r, _ar38cl / _ar40r, \
                   _ar38cl / (_ar39k * factor), _ar40r / (_ar39k * factor), \
                   _ar36res / _ar40ar, _ar38res / _ar40ar, (_ar39k * factor) / _ar40ar, \
-                  _ar36, _ar37, _ar38, _ar39, _ar40
+                  _ar36, _ar37, _ar38, _ar39, _ar40, _CaK
 
     data = np.array(list(do_simulation()))
     means = data.mean(axis=0)
@@ -467,8 +465,9 @@ def Monte_Carlo_F(ar40m: Tuple[float, float], ar39m: Tuple[float, float], ar38m:
     corrected_data = np.column_stack((means[29:34], stds[29:34])).ravel()
 
     F = np.append(means[25], stds[25])
+    CaK = np.append(means[34], stds[34])
 
-    return *F, *degas_data, *nor_iso, *inv_iso, *cl1_iso, *cl2_iso, *cl3_iso, *threed_iso, *corrected_data
+    return *F, *degas_data, *nor_iso, *inv_iso, *cl1_iso, *cl2_iso, *cl3_iso, *threed_iso, *corrected_data, *CaK
 
     # print("F = {0} ± {1}".format(np.mean(F), np.std(F)))
     #

ararpy/files/calc_file.py

@@ -81,7 +81,7 @@ def open_252(data: pd.DataFrame, logs01: pd.DataFrame, logs02: pd.DataFrame):
         162, 163, 164, 165, 166, 167, 168, 169,  # 40 r, a, c, k    26-33
     ]
     publish_values_index = [
-        101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111,
+        188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205,
     ]
     apparent_age_values_index = [
         198, 199, 200, 201, -999, -999, 202, 203,  # f, sf, ages, s, s, s, 40Arr%, 39Ar%
@@ -244,7 +244,9 @@ def open_240(data: pd.DataFrame, logs01: pd.DataFrame, logs02: pd.DataFrame):
         148, -999, 149, -999, 150, -999, 151, -999,
     ]
     publish_values_index = [
-        99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109,
+        -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
+        -1, -1, -1, -1,
+        -1, -1, -1, -1,
     ]
     apparent_age_values_index = [
         # f, sf, ages, s, s, s, 39Ar%

ararpy/smp/basic.py

@@ -119,7 +119,7 @@ def calc_apparent_ages(smp: Sample):
         raise TypeError(f"Sample type is not supported: {smp.Info.sample.type}")
 
     smp.ApparentAgeValues[2:6] = v[0:4]
-    smp.PublishValues[5:7] = copy.deepcopy(v[0:2])
+    smp.PublishValues[12:14] = copy.deepcopy(v[0:2])
 
 
 def calc_j(ar40ar39=None, params: dict = None, smp: Sample = None, index: list = None):
@@ -591,7 +591,7 @@ def get_diff_smp(backup: (dict, Sample), smp: (dict, Sample)):
             if _res != {}:
                 res.update({name: _res})
             continue
-        if str(backup[name]) == str(attr) or backup[name] == attr:
+        if str(backup[name]) == str(attr):
             continue
         res.update({name: attr})
     return res

ararpy/smp/corr.py

@@ -55,7 +55,6 @@ def corr_blank(sample: Sample):
     for i in range(0, 10, 2):
         blank_corrected[i] = [blank_corrected[i][index] if sample.TotalParam[102][index] else j for index, j in enumerate(sample.SampleIntercept[i])]
         blank_corrected[i + 1] = [blank_corrected[i + 1][index] if sample.TotalParam[102][index] else j for index, j in enumerate(sample.SampleIntercept[i + 1])]
-        blank_corrected[i] = [0 if j < 0 and sample.TotalParam[101][index] else j for index, j in enumerate(blank_corrected[i])]
     sample.BlankCorrected = blank_corrected
 
 
@@ -147,9 +146,6 @@ def corr_decay(sample: Sample):
             *sample.MassDiscrCorrected[2:4], t1, t2, t3, *sample.TotalParam[44:46], isRelative=True)
         decay_corrected[6:8] = calc.corr.decay(
             *sample.MassDiscrCorrected[6:8], t1, t2, t3, *sample.TotalParam[42:44], isRelative=True)
-        # Negative number set to zero in decay correction
-        decay_corrected[2] = [0 if i < 0 else i for i in decay_corrected[2]]
-        decay_corrected[6] = [0 if i < 0 else i for i in decay_corrected[6]]
     except Exception as e:
         print(traceback.format_exc())
         raise ValueError(f'Decay correction error: {str(e)}')
@@ -162,6 +158,10 @@ def corr_decay(sample: Sample):
     sample.CorrectedValues[6] = [val if corrDecay39[idx] else 0 for idx, val in enumerate(decay_corrected[6])]
     sample.CorrectedValues[7] = [val if corrDecay39[idx] else 0 for idx, val in enumerate(decay_corrected[7])]
 
+    data = np.array(sample.CorrectedValues)
+    data[1:10:2] = data[1:10:2] = np.abs(np.divide(data[1:10:2], data[0:10:2])) * 100
+    sample.PublishValues[0:10] = copy.deepcopy(data.tolist())
+
 
 # =======================
 # Degas Calcium derived 37Ar 36Ar 38Ar 39Ar
@@ -189,8 +189,10 @@ def calc_degas_ca(sample: Sample):
     if not validate_params(**params_to_check):
         return
 
+    set_negative_zero = sample.TotalParam[101]
     corrDegasCa = sample.TotalParam[106]
-    ar37ca = sample.CorrectedValues[2:4]
+    ar37ca = copy.deepcopy(sample.CorrectedValues[2:4])
+    ar37ca[0] = [0 if val < 0 and set_negative_zero[i] else val for i, val in enumerate(ar37ca[0])]
     ar39ca = calc.arr.mul_factor(ar37ca, sample.TotalParam[8:10], isRelative=True)
     ar38ca = calc.arr.mul_factor(ar37ca, sample.TotalParam[10:12], isRelative=True)
     ar36ca = calc.arr.mul_factor(ar37ca, sample.TotalParam[12:14], isRelative=True)
@@ -201,7 +203,6 @@ def calc_degas_ca(sample: Sample):
     sample.DegasValues[19] = [val if corrDegasCa[idx] else 0 for idx, val in enumerate(ar38ca[1])]
     sample.DegasValues[22] = [val if corrDegasCa[idx] else 0 for idx, val in enumerate(ar39ca[0])]  # 39Ca
     sample.DegasValues[23] = [val if corrDegasCa[idx] else 0 for idx, val in enumerate(ar39ca[1])]
-    sample.PublishValues[1] = copy.deepcopy(sample.DegasValues[8])
 
 
 # =======================
@@ -242,7 +243,6 @@ def calc_degas_k(sample: Sample):
     sample.DegasValues[17] = [val if corrDecasK[idx] else 0 for idx, val in enumerate(ar38k[1])]
     sample.DegasValues[30] = [val if corrDecasK[idx] else 0 for idx, val in enumerate(ar40k[0])]
     sample.DegasValues[31] = [val if corrDecasK[idx] else 0 for idx, val in enumerate(ar40k[1])]
-    sample.PublishValues[3] = copy.deepcopy(sample.DegasValues[20])
 
 
 # =======================
@@ -351,7 +351,6 @@ def calc_degas_cl(sample: Sample):
 
     sample.DegasValues[6:8] = copy.deepcopy(ar36cl)
     sample.DegasValues[10:12] = copy.deepcopy(ar38cl)
-    sample.PublishValues[2] = copy.deepcopy(sample.DegasValues[10])
 
 
 # =======================
@@ -408,7 +407,6 @@ def calc_degas_atm(sample: Sample):
     sample.DegasValues[13] = [val if corrDecasAtm[idx] else 0 for idx, val in enumerate(ar38a[1])]
     sample.DegasValues[26] = [val if corrDecasAtm[idx] else 0 for idx, val in enumerate(ar40a[0])]  # Ar40a
     sample.DegasValues[27] = [val if corrDecasAtm[idx] else 0 for idx, val in enumerate(ar40a[1])]
-    sample.PublishValues[0] = copy.deepcopy(sample.DegasValues[ 0])
 
 
 # =======================
@@ -439,7 +437,6 @@ def calc_degas_r(sample: Sample):
     ar40r = calc.arr.sub(ar40ar, sample.DegasValues[26:28])
     ar40r[0] = [0 if item < 0 and sample.TotalParam[101][index] else item for index, item in enumerate(ar40r[0])]
     sample.DegasValues[24:26] = copy.deepcopy(ar40r)
-    sample.PublishValues[4] = copy.deepcopy(sample.DegasValues[24])
 
 
 def calc_degas_c(sample: Sample):
@@ -588,7 +585,8 @@ def calc_ratio(sample: Sample):
         calc.arr.rec_factor(sample.TotalParam[20:22], isRelative=True))
 
     sample.ApparentAgeValues[6] = ar40r_percent
-    sample.PublishValues[7:11] = copy.deepcopy([*sample.ApparentAgeValues[6:8], *CaK])
+    sample.PublishValues[10:12] = copy.deepcopy(sample.ApparentAgeValues[0:2])
+    sample.PublishValues[14:18] = copy.deepcopy([*sample.ApparentAgeValues[6:8], *CaK])
 
     sample.IsochronValues[0:5], sample.IsochronValues[6:11] = calc_nor_inv_isochrons(sample)
     sample.IsochronValues[12:17], sample.IsochronValues[18:23], sample.IsochronValues[24:29] = \
@@ -629,13 +627,11 @@ def calc_ratio_monte_carlo(sample: Sample):
         # corrected
         sample.CorrectedValues = res[2 + 32 + 39:2 + 32 + 39 + 10]
         # publish
-        sample.PublishValues[0] = copy.deepcopy(sample.DegasValues[ 0])
-        sample.PublishValues[1] = copy.deepcopy(sample.DegasValues[ 8])
-        sample.PublishValues[2] = copy.deepcopy(sample.DegasValues[10])
-        sample.PublishValues[3] = copy.deepcopy(sample.DegasValues[20])
-        sample.PublishValues[4] = copy.deepcopy(sample.DegasValues[24])
-        sample.PublishValues[5:7] = copy.deepcopy(sample.ApparentAgeValues[2:4])
-        sample.PublishValues[7:9] = copy.deepcopy(sample.ApparentAgeValues[6:8])
+        data = np.array(sample.CorrectedValues)
+        data[1:10:2] = data[1:10:2] = np.abs(np.divide(data[1:10:2], data[0:10:2])) * 100
+        sample.PublishValues[0:10] = copy.deepcopy(data.tolist())
+        sample.PublishValues[10:14] = copy.deepcopy(sample.ApparentAgeValues[0:4])
+        sample.PublishValues[14:16] = copy.deepcopy(sample.ApparentAgeValues[6:8])
 
 
 def monte_carlo_f(sample: Sample):
@@ -731,6 +727,10 @@ def monte_carlo_f(sample: Sample):
     R38v36a = np.transpose(sample.TotalParam[4:6])
     R36v38clp = np.transpose(sample.TotalParam[56:58])
 
+    R36v38cl = np.transpose(sample.TotalParam[18:20])  # ?
+    KCaFactor = np.transpose(sample.TotalParam[20:22])
+    KClFactor = np.transpose(sample.TotalParam[22:24])
+
     stand_time_year = np.transpose(sample.TotalParam[32])
     JNFactor = np.transpose(sample.TotalParam[136:138])
 
@@ -764,6 +764,8 @@ def monte_carlo_f(sample: Sample):
             L37ar=L37ar[i], L39ar=L39ar[i], L36cl=L36cl[i],
             MDF=MDF[i], stand_time_year=stand_time_year[i],
             JNFactor=JNFactor[i],
+            KCaFactor=KCaFactor[i],
+            KClFactor=KClFactor[i],
             blank_gain_corr=sample.TotalParam[111][i],
             MDF_method=sample.TotalParam[100][i],
             force_to_zero=sample.TotalParam[101][i],

ararpy/smp/export.py

@@ -1219,7 +1219,7 @@ class WritingWorkbook:
             [(4, 11, 1), list(map(lambda x: 2*x, self.sample.ApparentAgeValues[3])), {'num_format': '± ' + self.default_fmt_prop['num_format']}],
             [(4, 12, 1), self.sample.ApparentAgeValues[6], {}],
             [(4, 13, 1), self.sample.ApparentAgeValues[7], {}],
-            [(4, 14, 1), self.sample.PublishValues[9], {}],
+            [(4, 14, 1), self.sample.PublishValues[16], {}],
             [(4 + num_step, 0, 0), [''] * 15, {'bold': 1, 'top': 1}],
 
             [(5 + num_step, 0), "Table 2. 40Ar/39Ar age summary", {'bold': 1, 'align': 'left'}],

ararpy/smp/initial.py

@@ -128,6 +128,8 @@ def initial(smp: Sample):
     smp.BlankIntercept = arr.create_arr((len(samples.BLANK_INTERCEPT_HEADERS) - 3, 0))
     smp.SampleIntercept = arr.create_arr((len(samples.SAMPLE_INTERCEPT_HEADERS) - 3, 0))
     smp.PublishValues = arr.create_arr((len(samples.PUBLISH_TABLE_HEADERS) - 3, 0))
+    smp.BlankCorrected = arr.create_arr((10, 0))
+    smp.MassDiscrCorrected = arr.create_arr((10, 0))
     smp.DecayCorrected = arr.create_arr((10, 0))
     smp.CorrectedValues = arr.create_arr((len(samples.CORRECTED_HEADERS) - 3, 0))
     smp.DegasValues = arr.create_arr((len(samples.DEGAS_HEADERS) - 3, 0))

ararpy/smp/sample.py

@@ -52,24 +52,36 @@ DEGAS_HEADERS = [
 ]
 PUBLISH_TABLE_HEADERS = [
     'Sequence', '', 'Mark',  # 0-2
-    '\u00B3\u2076Ar[a]', '\u00B3\u2077Ar[Ca]', '\u00B3\u2078Ar[Cl]',  # 2-4
-    '\u00B3\u2079Ar[K]', '\u2074\u2070Ar[r]',  # 5-6
-    'Apparent Age', '1\u03C3', '\u2074\u2070Ar[r]%', '\u00B3\u2079Ar[K]%',  # 7-10
-    'Ca/K', '1\u03C3',  # 11-12
+    '\u00B3\u2076Ar', '%1\u03C3',
+    '\u00B3\u2077Ar', '%1\u03C3',
+    '\u00B3\u2078Ar', '%1\u03C3',
+    '\u00B3\u2079Ar', '%1\u03C3',
+    '\u2074\u2070Ar', '%1\u03C3', # 3-12
+    '\u2074\u2070Ar/\u00B3\u2076Ar', '1\u03C3',  # 13-14
+    'Apparent Age', '1\u03C3', '\u2074\u2070Ar[r]%', '\u00B3\u2079Ar[K]%',  # 15-18
+    'Ca/K', '1\u03C3',  # 19-20
 ]
 PUBLISH_TABLE_HEADERS_AIR = [
     'Sequence', '', 'Mark',  # 0-2
-    '\u00B3\u2076Ar[a]', '\u00B3\u2077Ar[Ca]', '\u00B3\u2078Ar[Cl]',  # 2-4
-    '\u00B3\u2079Ar[K]', '\u2074\u2070Ar[r]',  # 5-6
-    'MDF', '1\u03C3', '\u2074\u2070Ar[r]%', '\u00B3\u2076Ar%',  # 7-10
-    'Ca/K', '1\u03C3',  # 11-12
+    '\u00B3\u2076Ar', '%1\u03C3',
+    '\u00B3\u2077Ar', '%1\u03C3',
+    '\u00B3\u2078Ar', '%1\u03C3',
+    '\u00B3\u2079Ar', '%1\u03C3',
+    '\u2074\u2070Ar', '%1\u03C3',  # 3-12
+    '\u2074\u2070Ar/\u00B3\u2076Ar', '1\u03C3',  # 13-14
+    'MDF', '1\u03C3', '\u2074\u2070Ar[r]%', '\u00B3\u2076Ar%',  # 15-18
+    'Ca/K', '1\u03C3',  # 19-20
 ]
 PUBLISH_TABLE_HEADERS_STD = [
     'Sequence', '', 'Mark',  # 0-2
-    '\u00B3\u2076Ar[a]', '\u00B3\u2077Ar[Ca]', '\u00B3\u2078Ar[Cl]',  # 2-4
-    '\u00B3\u2079Ar[K]', '\u2074\u2070Ar[r]',  # 5-6
-    'J', '1\u03C3', '\u2074\u2070Ar[r]%', '\u00B3\u2076Ar%',  # 7-10
-    'Ca/K', '1\u03C3',  # 11-12
+    '\u00B3\u2076Ar', '%1\u03C3',
+    '\u00B3\u2077Ar', '%1\u03C3',
+    '\u00B3\u2078Ar', '%1\u03C3',
+    '\u00B3\u2079Ar', '%1\u03C3',
+    '\u2074\u2070Ar', '%1\u03C3',  # 3-12
+    '\u2074\u2070Ar/\u00B3\u2076Ar', '1\u03C3',  # 13-14
+    'J', '1\u03C3', '\u2074\u2070Ar[r]%', '\u00B3\u2076Ar%',  # 15-18
+    'Ca/K', '1\u03C3',  # 19-20
 ]
 SPECTRUM_TABLE_HEADERS = [
     'Sequence', '', 'Mark',  # 0-2
@@ -200,9 +212,13 @@ DEGAS_SHORT_HEADERS = [
 ]
 PUBLISH_TABLE_SHORT_HEADERS = [
     'Seq', 'Label', 'Mark',  # 0-2
-    'Ar36[a]', 'Ar37[Ca]', 'Ar38[Cl]',
-    'Ar39[K]', 'Ar40[r]', 'Apparent Age', '1s', 'Ar40r%',
-    'Ar39K%', 'Ca/K', '1s'
+    'Ar36', '%1s',
+    'Ar37', '%1s',
+    'Ar38', '%1s',
+    'Ar39', '%1s',
+    'Ar40', '%1s',
+    'Apparent Age', '1s', 'Ar40r%', 'Ar39K%',
+    'Ca/K', '1s'
 ]
 SPECTRUM_TABLE_SHORT_HEADERS = [
     'Seq', 'Label', 'Mark',  # 0-2

ararpy/smp/table.py

@@ -12,6 +12,7 @@
 import ast
 import re
 import copy
+import numpy as np
 from .. import calc
 from . import (sample as samples, basic)
 
@@ -46,13 +47,11 @@ def update_table_data(smp: Sample, only_table: str = None):
         elif key == '4':
             data = [smp.SequenceName, smp.SequenceValue, smp.IsochronMark, *smp.DegasValues]
         elif key == '5':
-            smp.PublishValues[0] = copy.deepcopy(smp.DegasValues[ 0])
-            smp.PublishValues[1] = copy.deepcopy(smp.DegasValues[ 8])
-            smp.PublishValues[2] = copy.deepcopy(smp.DegasValues[10])
-            smp.PublishValues[3] = copy.deepcopy(smp.DegasValues[20])
-            smp.PublishValues[4] = copy.deepcopy(smp.DegasValues[24])
-            smp.PublishValues[5:7] = copy.deepcopy(smp.ApparentAgeValues[2:4])
-            smp.PublishValues[7:9] = copy.deepcopy(smp.ApparentAgeValues[6:8])
+            data = np.array(smp.CorrectedValues)
+            data[1:10:2] = np.abs(np.divide(data[1:10:2], data[0:10:2])) * 100
+            smp.PublishValues[0:10] = copy.deepcopy(data.tolist())
+            smp.PublishValues[10:14] = copy.deepcopy(smp.ApparentAgeValues[0:4])
+            smp.PublishValues[14:16] = copy.deepcopy(smp.ApparentAgeValues[6:8])
             data = [smp.SequenceName, smp.SequenceValue, smp.IsochronMark, *smp.PublishValues]
         elif key == '6':
             data = [smp.SequenceName, smp.SequenceValue, smp.IsochronMark, *smp.ApparentAgeValues]

{ararpy-0.2.2.dist-info → ararpy-0.2.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ararpy
-Version: 0.2.2
+Version: 0.2.3
 Summary: A project for Ar-Ar geochronology
 Home-page: https://github.com/wuyangchn/ararpy.git
 Author: Yang Wu

{ararpy-0.2.2.dist-info → ararpy-0.2.3.dist-info}/RECORD

@@ -5,13 +5,13 @@ ararpy/Example - Show MDD results.py,sha256=YFkiQual60lyCClsfvivr4REY6waSYuomf0u
 ararpy/Example - Show all Kfs age spectra.py,sha256=1edPtBpFfS0lC7vLnLu34mDHt_xjbc8ptk6rZ--pf-I,12304
 ararpy/Example - Show random walk results.py,sha256=8WWvbAI7ySiMR-XwtFe5kTQihW3mcqzL3S7EUZdoxYo,15477
 ararpy/Example - Tc calculation.py,sha256=sD9pu3IaZ8mBV95rV2wOEhlQUexFNqBUoBqnNMdUBis,19617
-ararpy/__init__.py,sha256=wz69XzhfcI0QMPU2QJ8gZoWEsG2vHeSE7b5E7GytwI0,6767
+ararpy/__init__.py,sha256=dCyvAK_sTG0KYdq42pz1cfnl8pXamf04nr46e2fDKK0,6767
 ararpy/test.py,sha256=4F46-JJ1Ge12HGae0qO44Qc6kiEMHBgn2MsY_5LlHDo,3973
 ararpy/calc/__init__.py,sha256=kUjRuLE8TLuKOv3i976RnGJoEMj23QBZDu37LWs81U4,322
 ararpy/calc/age.py,sha256=WOZs70zXiBWDIEhXJLIaNiYTOFJNk0NDbH5e5zCbCks,5435
 ararpy/calc/arr.py,sha256=x6z16N2H6QFYAPY11lnPRFZsZfwWu2tSEo87q7M9Ym4,14869
 ararpy/calc/basic.py,sha256=AMQrhkNHhHtiguqTe4GPOYLuV8TEJBXpRJXJAoOwD9w,4634
-ararpy/calc/corr.py,sha256=qh77w1uJ98x6b_N-0_u1seWc7YLtCRClpxG6fuYHx-M,21667
+ararpy/calc/corr.py,sha256=e9HGsD1GCS_Kv7X_dFYachBxAJgeKg2GP--fvZHedns,21609
 ararpy/calc/err.py,sha256=63LtprqjemlIb1QGDst4Ggcv5KMSDHdlAIL-nyQs1eA,2691
 ararpy/calc/histogram.py,sha256=0GVbDdsjd91KQ1sa2B7NtZ4KGo0XpRIJapgIrzAwQUo,5777
 ararpy/calc/isochron.py,sha256=ej9G2e68k6yszonWHsLcEubh3TA7eh1upTJP_X0ttAA,5726
@@ -42,32 +42,32 @@ ararpy/examples/sample-default.smp,sha256=YNkoQGgPrsL_fXS7ZHxfRtLQWekCDqT9czS6vB
 ararpy/files/__init__.py,sha256=l5B5ZQ01WdtvjjN0aMkyAFNgpwANdM_1I0tQbqnRuEY,69
 ararpy/files/arr_file.py,sha256=pD5MxkAydL7cNq2wmKFUaOU4jHhc3MzTYrwbxZ3f46w,881
 ararpy/files/basic.py,sha256=k2GXgZjhqSmKvpXQLjsXDksS_ZLvqD7AWW54fXnYvTI,1228
-ararpy/files/calc_file.py,sha256=CmZsgUj2t2FnLcCE38pbOqej7M-CbhQs6KMFXE3PbQA,29233
+ararpy/files/calc_file.py,sha256=-CZXmlHVf65hKMLhwp6eSewbxzxywXcnpycoEeLDhZA,29304
 ararpy/files/new_file.py,sha256=efblARIBROVLWS2w3-98BxLX5VZ8grRpiTkJFtf_rAk,214
 ararpy/files/raw_file.py,sha256=X2_RjlNOUbdkA9MM2ZMNlZmGL7jc1-eV-Wy-bIu2bb0,20592
 ararpy/files/xls.py,sha256=DVcZ_yRnc19p-m4leGGjt-YPDpSa2udYKmGyrM0qub0,640
 ararpy/smp/EXPORT_TO_PDF_DATA_PROPERTIES.py,sha256=baDM437tu6hsPv0uYfod0TREXlPd6kvMBFT1S9ZZlkk,3024
 ararpy/smp/__init__.py,sha256=k6_fa27UJsQK7K7oC5GYlwMo6l0Xd8af3QtOrZz2XJk,478
-ararpy/smp/basic.py,sha256=uJJjQajTzg4X7tg4loTgg8ilvWJqBhZIYmIXY_F-VUQ,28601
+ararpy/smp/basic.py,sha256=VLjNIlWh7mwVH9GBkOMOe7R7FX3fcTPXYSxF-SlMlt4,28579
 ararpy/smp/calculation.py,sha256=4Vg0HNbAGdTVdYIt17WeoKPmrTjmq9EqaP44P8AxT9U,2874
 ararpy/smp/consts.py,sha256=XIdjdz8cYxspG2jMnoItdlUsxr3hKbNFJjMZJh1bpzw,393
-ararpy/smp/corr.py,sha256=FBovb2T-pkrjq1OGVh0znYijQYw3BdRLTqoeHzUDHF4,38049
+ararpy/smp/corr.py,sha256=foNC_7bRc0UBZtDOntu_73W7DvtvTchhufnEHqC2jAA,37896
 ararpy/smp/diffusion_funcs.py,sha256=X6LC1XY-vub5BpAuZMlyA5k3DQcU8Qqf5uukwGjJOxw,186690
-ararpy/smp/export.py,sha256=nJhyPtrbOYR30JncePNKR87IpCzVgy4wyo6FZ0kJMEY,120438
+ararpy/smp/export.py,sha256=XnLKjWk4LFqtqXdy2qCUkOLy0vkjMq-X-tE9A545BVE,120439
 ararpy/smp/info.py,sha256=tpEIjrE4nR-GAYTricBk9gq0LuHh6F1Bt7HPo1rS2HM,497
-ararpy/smp/initial.py,sha256=hSDU3nptrECErDDgGcrLB0g7UmcnRshY_rv1YG0jdmw,19270
+ararpy/smp/initial.py,sha256=5YHvhMBnpXtcTH-cIV7ENhy0tzP03oKIjNMA4fg_Dc0,19374
 ararpy/smp/json.py,sha256=zfJCC_2LCDckqC8Fpu10jEA6Knl3UtKO31I5g4fvsBE,2273
 ararpy/smp/plots.py,sha256=hU3VyAyEkR_3enWvMHfuL_Ajjg-cK9f7SUHBKwHkmVQ,32577
 ararpy/smp/raw.py,sha256=77J1dEYL4ZSeftp4tyog_Cy1Y5bwHNaLqJK4i4KLOSY,6500
-ararpy/smp/sample.py,sha256=tCDCAS_KgJy4qN1nCGWTikZfsBu5EV1giccLI3gRY5g,60066
+ararpy/smp/sample.py,sha256=AqH3HJ0oQsUmbdPz7BFhzdt9UdpLf4FnlHMTNvqvgsU,60436
 ararpy/smp/style.py,sha256=kP1YnyAoVbKXfEARsrCiz0lTYipW3xEKqR4vIfuVNe4,8206
-ararpy/smp/table.py,sha256=9D1oGlYtoVgRte5lxOWYZbEManyudO5un-t7V6LQ-qg,7351
+ararpy/smp/table.py,sha256=doL4t2KnlZd9dTFjRz4r1s5sj76CoGSJ8UiHfJd_FV8,7218
 ararpy/thermo/__init__.py,sha256=6VBuqTRFl403PVqOuMkVrut0nKaQsAosBmfW91X1dMg,263
 ararpy/thermo/arrhenius.py,sha256=Ass1ichHfqIAtpv8eLlgrUc1UOb3Urh1qzr1E3gLB4U,233
 ararpy/thermo/atomic_level_random_walk.py,sha256=ncw9DtxRfS6zlQbLVLNX7WNoO9sX_nSomwAsTH0_O3k,25910
 ararpy/thermo/basic.py,sha256=JJRZbYmvXlpRAV2FeFPwLhrig4ZhNQmJnWqgOjo-1YQ,11508
-ararpy-0.2.2.dist-info/licenses/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
-ararpy-0.2.2.dist-info/METADATA,sha256=wvnw6khN1TWEJEGRbAunuJ3A95geGtqgvmAwCN7wGAc,24726
-ararpy-0.2.2.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-ararpy-0.2.2.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
-ararpy-0.2.2.dist-info/RECORD,,
+ararpy-0.2.3.dist-info/licenses/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
+ararpy-0.2.3.dist-info/METADATA,sha256=Zxw4VyknZZI6_KomWl7OlURNNpeLAyrI-vxPwGcRaqc,24726
+ararpy-0.2.3.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+ararpy-0.2.3.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
+ararpy-0.2.3.dist-info/RECORD,,