ararpy 0.1.20__tar.gz → 0.1.21__tar.gz

{ararpy-0.1.20 → ararpy-0.1.21}/PKG-INFO

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

{ararpy-0.1.20 → ararpy-0.1.21}/ararpy/__init__.py

@@ -16,7 +16,7 @@ from . import calc, smp, files, thermo, test
 """ Information """
 
 name = 'ararpy'
-version = '0.1.20'
+version = '0.1.21'
 __version__ = version
 full_version = version
 last_update = '2024-03-28'

{ararpy-0.1.20 → ararpy-0.1.21}/ararpy/calc/age.py

@@ -41,11 +41,13 @@ def calc_age_min(F, sF, **kwargs) -> tuple:
     J = arr.array_as_float(kwargs.pop('J'))
     sJ = arr.array_as_float(kwargs.pop('sJ') * J / 100)
     A = arr.array_as_float(kwargs.pop('A'))
-    sA = arr.array_as_float(kwargs.pop('sA') * A / 100)
+    sA = arr.array_as_float(kwargs.pop('sA') * A / 100)  # total A, A = Aec + (Ab+) + (Ab-). Ab- for Ca
     Ae = arr.array_as_float(kwargs.pop('Ae'))
-    sAe = arr.array_as_float(kwargs.pop('sAe') * Ae / 100)
+    sAe = arr.array_as_float(kwargs.pop('sAe') * Ae / 100)  # Aec
     Ab = arr.array_as_float(kwargs.pop('Ab'))
-    sAb = arr.array_as_float(kwargs.pop('sAb') * Ab / 100)
+    sAb = arr.array_as_float(kwargs.pop('sAb') * Ab / 100)  # Ab-
+    Abp = arr.array_as_float(kwargs.pop('Abp'))
+    sAbp = arr.array_as_float(kwargs.pop('sAbp') * Ab / 100)  # Ab+
     W = arr.array_as_float(kwargs.pop('W'))
     sW = arr.array_as_float(kwargs.pop('sW') * W / 100)
     Y = arr.array_as_float(kwargs.pop('Y'))
@@ -62,21 +64,21 @@ def calc_age_min(F, sF, **kwargs) -> tuple:
     sLb = arr.array_as_float(kwargs.pop('sLb') * Lb / 100)
     t = arr.array_as_float(kwargs.pop('t'))
     st = arr.array_as_float(kwargs.pop('st'))
-    # Ap = arr.array_as_float(kwargs.pop('Ap'))
-    # sAp = arr.array_as_float(kwargs.pop('sAp') * Ap / 100)
-    # Kp = arr.array_as_float(kwargs.pop('Kp'))
-    # sKp = arr.array_as_float(kwargs.pop('sKp') * Kp / 100)
+    Ap = arr.array_as_float(kwargs.pop('Ap'))
+    sAp = arr.array_as_float(kwargs.pop('sAp') * Ap / 100)
+    Kp = arr.array_as_float(kwargs.pop('Kp'))
+    sKp = arr.array_as_float(kwargs.pop('sKp') * Kp / 100)
 
-    # recalculation using Min et al.(2000) equation
+    # calculating using Min et al.(2000) equation
     # lmd = A * W * Y / (f * No)
     V = f * No / ((Ab + Ae) * W * Y)
-    sf = 0
     sV = pow((V / f * sf) ** 2 + (V / No * sNo) ** 2 + (V / (Ab + Ae)) ** 2 * (sAb ** 2 + sAe ** 2) +
              (V / W * sW) ** 2 + (V / Y * sY) ** 2, 0.5)
+
     # standard age in year, change to Ma
     t = t * 1000000
     st = st * 1000000
-    # back-calculating Ar40/Ar39 ration for the standard
+    # back-calculating Ar40/Ar39 of the standard
     stdR = (np.exp(t * L) - 1) / J
     # errors of standard age and decay constants were not applied
     sStdR = pow((stdR / J) ** 2 * sJ ** 2, 0.5)
@@ -85,31 +87,27 @@ def calc_age_min(F, sF, **kwargs) -> tuple:
     sR_1 = np.sqrt((sF / stdR) ** 2 + (F * sStdR / stdR ** 2) ** 2)  # errors of measured 40Ar/39Ar and J value
     sR_2 = np.sqrt((sF / stdR) ** 2)  # error of measured 40Ar/39Ar only
 
-    BB = 1
-    KK = np.exp(t / V) - 1  # 40Arr / 40K Use standard age
+    useStandardAge = True
+    if useStandardAge:
+        # ln part in Min 2000 equation
+        BB = 1
+        KK = np.exp(t / V) - 1  # 40Arr / 40K Use standard age
+    else:  # use Ar40* and K concentrations of primary standard
+        # not fanished, this function is wrong
+        BB = (Ab + Ae) / Ae
+        KK = Ap / Kp / f
+        raise TypeError(f"Not use standard age. Not supported.")
+
     XX = BB * KK * R + 1
     k0 = V * np.log(XX)
-    e1 = (np.log(XX) * V / f - V * BB * KK * R / (f * XX)) ** 2 * sf ** 2  # sF
-    e2 = (np.log(XX) * V / No) ** 2 * sNo ** 2  # sNo
-    e3 = (-1 * np.log(XX) * V / A + BB * KK * R / (A * XX)) ** 2 * sAb ** 2  # sAb
-    e4 = (-1 * np.log(XX) * V / A - Ab * KK * R / (Ae ** 2 * XX)) ** 2 * sAe ** 2  # sAe
-    e5 = (-1 * np.log(XX) * V / W - V * BB * KK * R / (W * XX)) ** 2 * sW ** 2  # sW
-    e6 = (np.log(XX) * V / Y) ** 2 * sY ** 2  # sY
-    e7 = (V * BB * KK / XX) ** 2 * sR_1 ** 2  # sR
-    # e8 = (V * BB * KK * R / (Ap * XX)) ** 2 * sAp ** 2  # sAp
-    # e9 = (V * BB * KK * R / (Kp * XX)) ** 2 * sKp ** 2  # sKp
-    e8, e9 = 0, 0
-    # useDecayConst = False
-    # if useDecayConst:  # k0 = log(L / Le * KK * R + 1) / L
-    #     e1 = (V * BB * KK * R / (f * XX)) ** 2 * sf ** 2
-    #     e2 = 0
-    #     e3 = (-1 * np.log(XX) * V / L + BB * KK * R / (L * XX)) ** 2 * sLb ** 2
-    #     e4 = (-1 * np.log(XX) * V / L - Lb * KK * R / (Le ** 2 * XX)) ** 2 * sLe ** 2
-    #     e5 = (V * BB * KK * R / (W * XX)) ** 2 * sW ** 2
-    #     e6 = 0
-    useStandardAge = True
+
     if useStandardAge:
-        e1, e9 = 0, 0
+        # k0 =  V * ln(exp(t/V) - 1 * R + 1)
+        e1 = (np.log(XX) - R * (KK + 1) * t / (V * XX)) ** 2 * sV ** 2  # sV
+        e2 = (BB * (KK + 1) * R / XX) ** 2 * st ** 2  # st
+        e3 = (V * BB * KK / XX) ** 2 * sR_1 ** 2  # sR
+    else:
+        e1, e2, e3 = 0, 0, 0
 
     # change to Ma
     # analytical error, error of 40Ar/39Ar only
@@ -117,7 +115,7 @@ def calc_age_min(F, sF, **kwargs) -> tuple:
     # internal error, errors of 40Ar/39Ar and J value
     s2 = np.sqrt((V * KK * BB / XX) ** 2 * sR_1 ** 2)
     # total external error
-    s3 = np.sqrt(e1 + e2 + e3 + e4 + e5 + e6 + e7 + e8 + e9)
+    s3 = np.sqrt(e1 + e2 + e3)
     age = k0
     return age, s1, s2, s3
 
@@ -146,7 +144,7 @@ def calc_age_general(F, sF, J, sJ, L, sL, **kwargs):
     v1 = sF ** 2 * (J / (L * (1 + J * F))) ** 2
     v2 = sJ ** 2 * (F / (L * (1 + J * F))) ** 2
     v3 = sL ** 2 * (np.log(1 + J * F) / (L ** 2)) ** 2
-    s1 = v1 ** .5  # analytical error
-    s2 = (v1 + v2) ** .5  # internal error
-    s3 = (v1 + v2 + v3) ** .5  # full external error
+    s1 = v1 ** .5  # analytical error, F only
+    s2 = (v1 + v2) ** .5  # internal error, F and J
+    s3 = (v1 + v2 + v3) ** .5  # full external error. F, J and L
     return age, s1, s2, s3

{ararpy-0.1.20 → ararpy-0.1.21}/ararpy/smp/basic.py

@@ -153,7 +153,8 @@ def calc_age(ar40ar39=None, params: dict = None, smp: Sample = None, index: list
     """
     params_index_dict = {
         34: 'L', 35: 'sL', 36: 'Le', 37: 'sLe', 38: 'Lb', 39: 'sLb', 48: 'A', 49: 'sA',
-        50: 'Ae', 51: 'sAe', 52: 'Ab', 53: 'sAb', 59: 't', 60: 'st', 67: 'J', 68: 'sJ',
+        50: 'Ae', 51: 'sAe', 52: 'Ab', 53: 'sAb', 54: 'Abp', 55: 'sAbp', 59: 't', 60: 'st',
+        61: 'Ap', 62: 'sAp', 63: 'Kp', 64: 'sKp', 67: 'J', 68: 'sJ',
         81: 'W', 82: 'sW', 83: 'No', 84: 'sNo', 85: 'Y', 86: 'sY', 87: 'f', 88: 'sf', 110: 'Min'
     }
 

{ararpy-0.1.20 → ararpy-0.1.21}/ararpy/smp/corr.py

@@ -545,7 +545,7 @@ def monte_carlo_f(sample: Sample):
     R36v38clp = np.transpose(sample.TotalParam[56:58])
 
     stand_time_year = np.transpose(sample.TotalParam[32])
-    JNFactor = sample.NormalizeFactor
+    JNFactor = np.transpose(sample.NormalizeFactor)
 
     irradiation_cycles = [list(filter(None, re.split(r'[DS]', each_step))) for each_step in sample.TotalParam[27]]
     t1 = [re.findall(r"\d+", i) for i in sample.TotalParam[31]]  # t1: experimental times
@@ -579,7 +579,7 @@ def monte_carlo_f(sample: Sample):
             L37ar=L37ar[i], L39ar=L39ar[i], L36cl=L36cl[i],
             MDFunc=None,
             MDF=MDF[i], stand_time_year=stand_time_year[i],
-            JNFactor=JNFactor,
+            JNFactor=JNFactor[i],
             blank_gain_corr=sample.TotalParam[111][i]
         )
 

{ararpy-0.1.20 → ararpy-0.1.21}/ararpy/smp/sample.py

@@ -248,8 +248,8 @@ TOTAL_PARAMS_SHORT_HEADERS = [
     'DAb-K40', '%1s',  # 54-55
     'DAb+K40', '%1s',  # 56-57
     'Cl36/Cl38P', '%1s',  # 58-59
-    'StdName', 'StdAge', '1s', 'Ar40%', '1s', 'K%', '1s',
-    'F', '1s',  # 60-68
+    'StdName', 'StdAge', '1s', 'Ar40%', '1s', 'K%', '1s',  # 60-66
+    'F', '1s',  # 67-68
     'J', '%1s', 'MDF', '%1s',  # 69-72
     'MassAr36', '%1s', 'MassAr37', '%1s',
     'MassAr38', '%1s', 'MassAr39', '%1s',

{ararpy-0.1.20 → ararpy-0.1.21}/ararpy.egg-info/PKG-INFO

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

{ararpy-0.1.20 → ararpy-0.1.21}/setup.py

@@ -16,7 +16,7 @@ long_description = (here / 'README.md').read_text(encoding='utf-8')
 
 setuptools.setup(
     name='ararpy',  #
-    version='0.1.20',  # version
+    version='0.1.21',  # version
     author='Yang Wu',
     author_email='wuycug@hotmail.com',
     description='A project for Ar-Ar geochronology',  # short description