ararpy 0.1.25__py3-none-any.whl → 0.1.27__py3-none-any.whl

ararpy/__init__.py

@@ -16,10 +16,10 @@ from . import calc, smp, files, thermo, test
 """ Information """
 
 name = 'ararpy'
-version = '0.1.25'
+version = '0.1.26'
 __version__ = version
 full_version = version
-last_update = '2024-05-20'
+last_update = '2024-05-22'
 
 """ ArArPy Functions """
 

ararpy/calc/basic.py

@@ -129,3 +129,11 @@ def monte_carlo(func, inputs, confidence_level, **kwargs):
     limits = values[:, [l_range, r_range]]
 
     return np.concatenate((means, limits), axis=1), cov
+
+
+def parser(v: str, t: type):
+    if t == bool:
+        if str(v).lower().rstrip() in ["true", "1"]:
+            return True
+    return t(v)
+

ararpy/files/calc_file.py

@@ -22,7 +22,7 @@ import os
 import msoffcrypto
 
 from . import xls
-from ..calc.basic import get_datetime
+from ..calc.basic import get_datetime, parser
 from ..calc import arr, err, isochron
 
 
@@ -107,9 +107,8 @@ def open_252(data: pd.DataFrame, logs01: pd.DataFrame, logs02: pd.DataFrame):
         -999, -999, -999, -999, -999, -999, -999, -999, -999, -999,  # 90-99
         -999, -999, -999, -999, -999, -999, -999, -999, -999, -999,  # 100-109
         -999, -999, -999, -999, -1,  # 110-114
-        -999, -999, -999, -999, -999, -999, -999, -999,  # 115-122
-        -999, -999, -999, -999, -999, -999, -999, -999,  # 123-130
-        -999, -999, -999, -999, -999,  # 131-135
+        -999, -999, -999, -999, -999, -999, -999, -999, -999, -999, -999,  # 115-125
+        999, -999, 999, -999, 999, -999, 999, -999, 999, -999,  # 126-135
     ]
 
     # double transpose to remove keys
@@ -128,6 +127,13 @@ def open_252(data: pd.DataFrame, logs01: pd.DataFrame, logs02: pd.DataFrame):
 
     # adjustment
     isochron_mark.replace({4.0: '1'}, inplace=True)
+
+    total_param[71] = 35.96754628
+    total_param[73] = 36.9667759
+    total_param[75] = 37.9627322
+    total_param[77] = 38.964313
+    total_param[79] = 39.962383123
+
     total_param[83] = logs01[1][9]  # No
     total_param[84] = logs01[1][10]  # %SD
     total_param[81] = logs01[1][11]  # K mass
@@ -150,9 +156,9 @@ def open_252(data: pd.DataFrame, logs01: pd.DataFrame, logs02: pd.DataFrame):
     total_param[51] = logs01[1][27]  # %SD
     total_param[52] = logs01[1][28]  # 40K(\beta-) activities param
     total_param[53] = logs01[1][29]  # %SD
-    total_param[48] = logs01[1][26] + logs01[1][28]  # 40K(EC) activities param
-    total_param[49] = 100 / total_param[48] * pow(
-        (logs01[1][27] / 100 * logs01[1][26]) ** 2 + (logs01[1][29] / 100 * logs01[1][28]) ** 2, 0.5)  # %SD
+    total_param[48] = logs01[1][26] + logs01[1][28]  # 40K total activities param
+    total_param[49] = round(100 / total_param[48] * pow(
+        (logs01[1][27] / 100 * logs01[1][26]) ** 2 + (logs01[1][29] / 100 * logs01[1][28]) ** 2, 0.5), 6)  # %SD
 
     total_param[34] = logs01[1][36]  # decay constant of 40K total
     total_param[35] = logs01[1][37]  # %SD
@@ -169,30 +175,29 @@ def open_252(data: pd.DataFrame, logs01: pd.DataFrame, logs02: pd.DataFrame):
     total_param[44] = logs01[1][40]  # decay constant of 37Ar
     total_param[45] = logs01[1][41]  # %SD
 
-    # logs01[1][49]]  # 1 for using weighted YORK isochron regression, 2 for unweighted
-    # logs01[1][50]]  # True for including errors on irradiation constants, False for excluding those
-    # [logs01[1][67]]  # MDF method, 'LIN'
+    # logs01[1][49]  # 1 for using weighted YORK isochron regression, 2 for unweighted
+    # logs01[1][50]  # True for including errors on irradiation constants, False for excluding those
+    # logs01[1][67]  # MDF method, 'LIN'
     total_param[93] = logs01[1][70]  # 40Ar/36Ar air
     total_param[94] = logs01[1][71]  # %SD
-
-    total_param[97] = 'York-2'  # Fitting method
+    total_param[97] = 'York-2' if int(logs01[1][49]) == 1 else ""   # Fitting method
     total_param[98] = logs01[1][0]  # Convergence
     total_param[99] = logs01[1][1]  # Iterations number
     total_param[100] = logs01[1][67]  # MDF method
-    total_param[101] = logs01[1][65]  # force negative to zero when correct blank
+    total_param[101] = parser(logs01[1][65], bool)  # force negative to zero when correct blank
     total_param[102] = True  # Apply 37Ar decay
     total_param[103] = True  # Apply 39Ar decay
-    total_param[104] = logs01[1][6]  # Apply 37Ar decay
-    total_param[105] = logs01[1][7]  # Apply 39Ar decay
-    total_param[106] = True  # Apply K degas
-    total_param[107] = True  # Apply Ca degas
-    total_param[108] = True  # Apply Air degas
-    total_param[109] = logs01[1][8]  # Apply Cl degas
+    total_param[104] = parser(logs01[1][6], bool)  # Apply 37Ar decay
+    total_param[105] = parser(logs01[1][7], bool)  # Apply 39Ar decay
+    total_param[106] = True  # Apply Ca degas
+    total_param[107] = True  # Apply K degas
+    total_param[108] = parser(logs01[1][8], bool)  # Apply Cl degas
+    total_param[109] = True  # Apply Trap degas
 
     total_param[110] = True if logs01[1][5] == 'MIN' else False  # Calculating ages using Min equation [True] or conventional equation [False]
-    total_param[111] = logs01[1][89]  # Use primary standard or not
-    total_param[112] = logs01[1][91]  # Use standard age or not
-    total_param[113] = logs01[1][92]  # Use primary ratio or not
+    total_param[111] = parser(logs01[1][89], bool)  # Use primary standard or not
+    total_param[112] = parser(logs01[1][91], bool)  # Use standard age or not
+    total_param[113] = parser(logs01[1][92], bool)  # Use primary ratio or not
 
     # irradiation time
     total_param = general_adjustment(
@@ -310,12 +315,12 @@ def open_240(data: pd.DataFrame, logs01: pd.DataFrame, logs02: pd.DataFrame):
     )
     # Do adjustment
     isochron_mark.replace({4.0: '1'}, inplace=True)
-    total_param[44] = logs01[1][40]
-    total_param[45] = logs01[1][41]
-    total_param[42] = logs01[1][38]
-    total_param[43] = logs01[1][39]
     total_param[34] = logs01[1][36]
     total_param[35] = logs01[1][37]
+    total_param[42] = logs01[1][38]
+    total_param[43] = logs01[1][39]
+    total_param[44] = logs01[1][40]
+    total_param[45] = logs01[1][41]
 
     return [
         sample_values, blank_values, corrected_values, degas_values, publish_values,
@@ -399,6 +404,8 @@ def general_adjustment(
     total_param[32] = [i / (3600 * 24 * 365.242) for i in stand_time_second]  # stand year
 
     # initial ratios & error display settings
+    total_param[95] = logs02[9][irradiation_index]  # 38Ar/36Ar air
+    total_param[96] = logs02[10][irradiation_index]  # %SD
     total_param[115] = 0
     total_param[116] = 298.56
     total_param[117] = 0.31
@@ -624,10 +631,10 @@ class ArArCalcFile:
         # list(range(1, 26, 2)) irradiation correction constants
         # list(range(39, 58, 2)) decay constants
         # list(range(68, 97, 2)) J, MDF, other constants
-        for column in list(range(1, 26, 2)) + list(range(68, 71, 2)):
+        for column in list(range(1, 26, 2)) + list(range(68, 71, 2)):  # irradiation constants and J, MDF uncertainties
             self.content.loc[:, ('pam', column)] = \
-                (self.content['pam', column].astype("float") / (
-                    self.content['pam', column - 1].astype("float")) * 100).replace(np.nan, 0)
+                round(self.content['pam', column].astype("float") / (
+                    self.content['pam', column - 1].astype("float")) * 100, 6).replace(np.nan, 0)
 
         # sample info
         self.sample_info = {

ararpy/files/raw_file.py

@@ -25,7 +25,43 @@ from ..calc.basic import utc_dt
 
 """ Open raw data file """
 
-DEFAULT_SAMPLE_INFO = {}
+DEFAULT_SAMPLE_INFO = {
+    "ExpName": "en",        # "Experiment Name"
+    "StepName": "sn",       # "Step Name"
+    "StepType": "st",       # "Step Type"
+    "StepLabel": "sl",      # "Step Label"
+    "ZeroYear": "date",     # "Zero Date Year"
+    "ZeroHour": "time",     # "Zero Time Hour"
+    "ZeroMon": "mon",       # "Zero Date Month"
+    "ZeroMin": "min",       # "Zero Time Minute"
+    "ZeroDay": "day",       # "Zero Date Day"
+    "ZeroSec": "sec",       # "Zero Time Second"
+    "SmpName": "smpn",      # "Sample Name"
+    "SmpLoc": "smpl",       # "Sample location"
+    "SmpMatr": "smpm",      # "Sample Material"
+    "ExpType": "et",        # "Experiment Type"
+    "SmpWeight": "smpw",    # "Sample Weight"
+    "Stepunit": "su",       # "Step unit"
+    "HeatingTime": "ht",    # "Heating time"
+    "InstrName": "inn",     # "Instrument name"
+    "Researcher": "re",     # "Researcher"
+    "Analyst": "an",        # "Analyst"
+    "Lab": "la",            # "Laboratory"
+    "Jv": "jv",             # "J value"
+    "Jsig": "je",           # "J value error"
+    "MDF": "mdf",           # "MDF"
+    "MDFSig": "mdfe",       # "MDF error"
+    "CalcName": "cp",       # "Calc params"
+    "IrraName": "in",       # "Irra name"
+    "IrraLabel": "il",      # "Irra label"
+    "IrraPosH": "ih",       # "Irra position H"
+    "IrraPosX": "ix",       # "Irra position X"
+    "IrraPosY": "iy",       # "Irra position Y"
+    "StdName": "stdn",      # "Standard name"
+    "StdAge": "stda",       # "Standard age"
+    "StdAgeSig": "stde",    # "Standard age error"
+}
+
 
 
 def open_file(file_path: str, input_filter: List[Union[str, int, bool]]):
@@ -61,7 +97,6 @@ def open_file(file_path: str, input_filter: List[Union[str, int, bool]]):
                    'Qtegra Exported XLS': open_qtegra_exported_xls, 'Seq': open_raw_seq}[
             ['txt', 'excel', 'Qtegra Exported XLS', 'Seq'][int(input_filter[1])]]
     except KeyError:
-        print(traceback.format_exc())
         raise FileNotFoundError("Wrong File.")
     return handler(file_path, input_filter)
 
@@ -217,48 +252,27 @@ def get_raw_data(file_contents: List[List[Union[int, float, str, bool, list]]],
     idx = step_index = 0
 
     header = input_filter[5]
-    isotopic_data_index = input_filter[8:28]
-    sample_info_index = input_filter[33:65]
-    optional_info_index = input_filter[37:-6]
-    check_box_index = input_filter[-6:]
+    isotope_index = input_filter[8:28]
+    data_index = input_filter[4:33]
+    strings_index = input_filter[33:41]         # strings: filename, parsing information, datetime, timezone
+    optional_info_index = input_filter[41:-7]       # from Exp Name to Std Age Error
+    check_box_index = input_filter[-7:]
 
-    timezone = sample_info_index[3] if sample_info_index[3] != "" else "utc"
     while True:  # measurment steps sloop
 
-        # all text information
+        # ============ all text information ============
         options = get_sample_info(file_contents, optional_info_index, default="", base=[1, 1 - idx, 1])
 
-        # Zero datetime
-        try:
-            if check_box_index[2]:  # input_filter[134]: date in one string
-                zero_date = datetime_parse(options.get('ZeroYear'), sample_info_index[1])
-            else:
-                zero_date = datetime(year=options.get('ZeroYear'), month=options.get('ZeroMon'),
-                                     day=options.get('ZeroDay'))
-
-            if check_box_index[3]:  # input_filter[135]: time in one string
-                zero_time = datetime_parse(options.get('ZeroHour'), sample_info_index[2])
-            else:
-                zero_time = datetime(year=2020, month=12, day=31, hour=options.get('ZeroHour'),
-                                     minute=options.get('ZeroMin'), second=options.get('ZeroSec'))
-
-            zero_datetime = datetime(
-                zero_date.year, zero_date.month, zero_date.day, zero_time.hour, zero_time.minute, zero_time.second)
-            # adjust to UTC
-            zero_datetime = utc_dt(zero_datetime, tz=timezone).isoformat(timespec='seconds')
-        except (TypeError, ValueError, IndexError):
-            # print(f"Cannot parse zero datetime")
-            zero_datetime = datetime(1970, 1, 1, 0, 0, 0).isoformat(timespec='seconds')
-
-        # Step name
+        # ============ Step name ============
         try:
             step_name = options.get('StepName')
             experiment_name = options.get('ExpName')
-            if check_box_index[1] and sample_info_index[0].strip() != "":
-                _res = string_parser(sample_info_index[0], file_name)
-                if _res is not None:
-                    experiment_name = _res.named.get("en", experiment_name)
-                    step_index = _res.named.get("sn", step_name)
+            if check_box_index[1] and strings_index[0].strip() != "":
+                res = string_parser(strings_index[0], file_name)
+                if res is not None and "en" in res.named.keys():
+                    experiment_name = res.named.get("en")
+                if res is not None and "sn" in res.named.keys():
+                    step_index = res.named.get("sn")
                     if step_index.isnumeric():
                         step_name = f"{experiment_name}-{int(step_index):02d}"
                     else:
@@ -273,33 +287,71 @@ def get_raw_data(file_contents: List[List[Union[int, float, str, bool, list]]],
             options.update({'StepName': step_name})
             options.update({'ExpName': experiment_name})
 
+        # ============ Step information ============
+        try:
+            if check_box_index[2]:
+                string = get_item(file_contents, strings_index[1:4], default="", base=[1, 1 - idx, 1])
+                res = string_parser(strings_index[4], string)
+                if res is not None:
+                    options.update(dict(zip(DEFAULT_SAMPLE_INFO.keys(), [res.named.get(value, options.get(key)) for key, value in DEFAULT_SAMPLE_INFO.items()])))
+        except (TypeError, ValueError, IndexError):
+            break
+        else:
+            pass
+
+        # ============ Zero datetime ============
+        timezone = strings_index[7] if strings_index[7] != "" else "utc"
+        options.update({'TimeZone': timezone})
+        try:
+            if check_box_index[3]:  # date in one string
+                zero_date = datetime_parse(options.get('ZeroYear'), strings_index[5])
+            else:
+                zero_date = datetime(year=options.get('ZeroYear'), month=options.get('ZeroMon'),
+                                     day=options.get('ZeroDay'))
+
+            if check_box_index[4]:  # time in one string
+                zero_time = datetime_parse(options.get('ZeroHour'), strings_index[6])
+            else:
+                zero_time = datetime(year=2020, month=12, day=31, hour=options.get('ZeroHour'),
+                                     minute=options.get('ZeroMin'), second=options.get('ZeroSec'))
+
+            zero_datetime = datetime(
+                zero_date.year, zero_date.month, zero_date.day, zero_time.hour, zero_time.minute, zero_time.second)
+            # adjust to UTC
+            zero_datetime = utc_dt(zero_datetime, tz=timezone).isoformat(timespec='seconds')
+            options.update({'ZeroDT': zero_datetime})
+        except (TypeError, ValueError, IndexError) as e:
+            print(traceback.format_exc())
+            raise ValueError(f"Failed to parse zero datetime: {e}")
+            # zero_datetime = datetime(1970, 1, 1, 0, 0, 0).isoformat(timespec='seconds')
         current_step = [[step_name, zero_datetime, options]]
 
+        # ============ isotope data ============
         break_num = 0
         cycle_num = 0
-        f = float(input_filter[31])  # Intensity Scale Factor
-        data_content = file_contents[input_filter[4] - 1 if input_filter[4] != 0 else 0]
+        f = float(data_index[27])  # Intensity Scale Factor
+        data_content = file_contents[data_index[0] - 1 if data_index[0] != 0 else 0]
         for i in range(2000):  # measurement cycle sloop
-            if break_num < input_filter[29]:
+            if break_num < data_index[25]:
                 break_num += 1
                 continue
             break_num = 0
-            if int(input_filter[6]) == 0:  # == 0, vertical
-                start_row = input_filter[28] * cycle_num + input_filter[29] * cycle_num + header + idx - 1
+            if int(data_index[2]) == 0:  # == 0, vertical
+                start_row = data_index[24] * cycle_num + data_index[25] * cycle_num + header + idx - 1
                 try:
                     current_step.append([
                         str(cycle_num + 1),
                         # in sequence: Ar36, Ar37, Ar38, Ar39, Ar40
-                        float(data_content[start_row + isotopic_data_index[18]][isotopic_data_index[19] - 1]),
-                        float(data_content[start_row + isotopic_data_index[16]][isotopic_data_index[17] - 1]) * f,
-                        float(data_content[start_row + isotopic_data_index[14]][isotopic_data_index[15] - 1]),
-                        float(data_content[start_row + isotopic_data_index[12]][isotopic_data_index[13] - 1]) * f,
-                        float(data_content[start_row + isotopic_data_index[10]][isotopic_data_index[11] - 1]),
-                        float(data_content[start_row + isotopic_data_index[ 8]][isotopic_data_index[ 9] - 1]) * f,
-                        float(data_content[start_row + isotopic_data_index[ 6]][isotopic_data_index[ 7] - 1]),
-                        float(data_content[start_row + isotopic_data_index[ 4]][isotopic_data_index[ 5] - 1]) * f,
-                        float(data_content[start_row + isotopic_data_index[ 2]][isotopic_data_index[ 3] - 1]),
-                        float(data_content[start_row + isotopic_data_index[ 0]][isotopic_data_index[ 1] - 1]) * f,
+                        float(data_content[start_row + isotope_index[18]][isotope_index[19] - 1]),
+                        float(data_content[start_row + isotope_index[16]][isotope_index[17] - 1]) * f,
+                        float(data_content[start_row + isotope_index[14]][isotope_index[15] - 1]),
+                        float(data_content[start_row + isotope_index[12]][isotope_index[13] - 1]) * f,
+                        float(data_content[start_row + isotope_index[10]][isotope_index[11] - 1]),
+                        float(data_content[start_row + isotope_index[ 8]][isotope_index[ 9] - 1]) * f,
+                        float(data_content[start_row + isotope_index[ 6]][isotope_index[ 7] - 1]),
+                        float(data_content[start_row + isotope_index[ 4]][isotope_index[ 5] - 1]) * f,
+                        float(data_content[start_row + isotope_index[ 2]][isotope_index[ 3] - 1]),
+                        float(data_content[start_row + isotope_index[ 0]][isotope_index[ 1] - 1]) * f,
                     ])
                 except (ValueError, IndexError):
                     print(f"Cannot parse isotope data")
@@ -307,23 +359,23 @@ def get_raw_data(file_contents: List[List[Union[int, float, str, bool, list]]],
                     current_step.append([
                         str(cycle_num + 1), None, None, None, None, None, None, None, None, None, None,
                     ])
-            elif int(input_filter[6]) == 1:  # == 1, horizontal
-                start_row = input_filter[5] + idx
-                col_inc = input_filter[28] * cycle_num + input_filter[29] * cycle_num - 1
+            elif int(data_index[2]) == 1:  # == 1, horizontal
+                start_row = data_index[1] + idx
+                col_inc = data_index[24] * cycle_num + data_index[25] * cycle_num - 1
                 try:
                     current_step.append([
                         str(cycle_num + 1),
                         # Ar36, Ar37, Ar38, Ar39, Ar40
-                        float(data_content[start_row][isotopic_data_index[19] + col_inc]),
-                        float(data_content[start_row][isotopic_data_index[17] + col_inc]) * f,
-                        float(data_content[start_row][isotopic_data_index[15] + col_inc]),
-                        float(data_content[start_row][isotopic_data_index[13] + col_inc]) * f,
-                        float(data_content[start_row][isotopic_data_index[11] + col_inc]),
-                        float(data_content[start_row][isotopic_data_index[ 9] + col_inc]) * f,
-                        float(data_content[start_row][isotopic_data_index[ 7] + col_inc]),
-                        float(data_content[start_row][isotopic_data_index[ 5] + col_inc]) * f,
-                        float(data_content[start_row][isotopic_data_index[ 3] + col_inc]),
-                        float(data_content[start_row][isotopic_data_index[ 1] + col_inc]) * f,
+                        float(data_content[start_row][isotope_index[19] + col_inc]),
+                        float(data_content[start_row][isotope_index[17] + col_inc]) * f,
+                        float(data_content[start_row][isotope_index[15] + col_inc]),
+                        float(data_content[start_row][isotope_index[13] + col_inc]) * f,
+                        float(data_content[start_row][isotope_index[11] + col_inc]),
+                        float(data_content[start_row][isotope_index[ 9] + col_inc]) * f,
+                        float(data_content[start_row][isotope_index[ 7] + col_inc]),
+                        float(data_content[start_row][isotope_index[ 5] + col_inc]) * f,
+                        float(data_content[start_row][isotope_index[ 3] + col_inc]),
+                        float(data_content[start_row][isotope_index[ 1] + col_inc]) * f,
                     ])
                 except (ValueError, IndexError):
                     print(f"Cannot parse isotope data")
@@ -332,18 +384,21 @@ def get_raw_data(file_contents: List[List[Union[int, float, str, bool, list]]],
                         str(cycle_num + 1), None, None, None, None, None, None, None, None, None, None,
                     ])
             else:
-                raise ValueError(f"{input_filter[6]} not in [0, 1]")
+                raise ValueError(f"{data_index[2]} not in [0, 1]")
 
             cycle_num += 1
-            if cycle_num >= input_filter[7]:
+            if cycle_num >= data_index[3]:
                 break
 
         step_list.append(current_step)
-        idx = input_filter[32] * len(step_list)
+        idx = data_index[28] * len(step_list)
 
         if not check_box_index[0] or len(step_list) >= 500:  # check_box_index[0]: multiple sequences
             break
 
+    if not step_list:
+        raise ValueError("Failed to read the original file. It might be because the names of the experiments or steps were not recognized.")
+
     return step_list
 
 
@@ -366,12 +421,12 @@ def get_sample_info(file_contents: list, index_list: list, default="", base: Uni
         "StepName": get_item(file_contents, index_list[3:6], default=default, base=base),
         "StepType": get_item(file_contents, index_list[6:9], default=default, base=base),
         "StepLabel": get_item(file_contents, index_list[9:12], default=default, base=base),
-        "ZeroYear": get_item(file_contents, index_list[12:15], default=default, base=base),  # year
-        "ZeroHour": get_item(file_contents, index_list[15:18], default=default, base=base),  # hour
-        "ZeroMon": get_item(file_contents, index_list[18:21], default=default, base=base),  # month
-        "ZeroMin": get_item(file_contents, index_list[21:24], default=default, base=base),  # minute
-        "ZeroDay": get_item(file_contents, index_list[24:27], default=default, base=base),  # day
-        "ZeroSec": get_item(file_contents, index_list[27:30], default=default, base=base),  # second
+        "ZeroYear": get_item(file_contents, index_list[12:15], default=default, base=base),
+        "ZeroHour": get_item(file_contents, index_list[15:18], default=default, base=base),
+        "ZeroMon": get_item(file_contents, index_list[18:21], default=default, base=base),
+        "ZeroMin": get_item(file_contents, index_list[21:24], default=default, base=base),
+        "ZeroDay": get_item(file_contents, index_list[24:27], default=default, base=base),
+        "ZeroSec": get_item(file_contents, index_list[27:30], default=default, base=base),
         "SmpName": get_item(file_contents, index_list[30:33], default=default, base=base),
         "SmpLoc": get_item(file_contents, index_list[33:36], default=default, base=base),
         "SmpMatr": get_item(file_contents, index_list[36:39], default=default, base=base),
@@ -385,46 +440,16 @@ def get_sample_info(file_contents: list, index_list: list, default="", base: Uni
         "Lab": get_item(file_contents, index_list[60:63], default=default, base=base),
         "Jv": get_item(file_contents, index_list[63:66], default=default, base=base),
         "Jsig": get_item(file_contents, index_list[66:69], default=default, base=base),
-        "CalcName": get_item(file_contents, index_list[69:72], default=default, base=base),
-        "IrraName": get_item(file_contents, index_list[72:75], default=default, base=base),
-        "IrraLabel": get_item(file_contents, index_list[75:78], default=default, base=base),
-        "IrraPosH": get_item(file_contents, index_list[78:81], default=default, base=base),
-        "IrraPosX": get_item(file_contents, index_list[81:84], default=default, base=base),
-        "IrraPosY": get_item(file_contents, index_list[84:87], default=default, base=base),
-        "StdName": get_item(file_contents, index_list[87:90], default=default, base=base),
-        "StdAge": get_item(file_contents, index_list[90:93], default=default, base=base),
-        "StdAgeSig": get_item(file_contents, index_list[93:96], default=default, base=base),
-        # "Experiment Name": get_item(file_contents, sample_info_index[0:3], default=default, base=base),
-        # "Step Name": get_item(file_contents, sample_info_index[3:6], default=default, base=base),
-        # "Step Type": get_item(file_contents, sample_info_index[6:9], default=default, base=base),
-        # "Step Label": get_item(file_contents, sample_info_index[9:12], default=default, base=base),
-        # "Zero Date Year": get_item(file_contents, sample_info_index[12:15], default=default, base=base),
-        # "Zero Time Hour": get_item(file_contents, sample_info_index[15:18], default=default, base=base),
-        # "Zero Date Month": get_item(file_contents, sample_info_index[18:21], default=default, base=base),
-        # "Zero Time Minute": get_item(file_contents, sample_info_index[21:24], default=default, base=base),
-        # "Zero Date Day": get_item(file_contents, sample_info_index[24:27], default=default, base=base),
-        # "Zero Time Second": get_item(file_contents, sample_info_index[27:30], default=default, base=base),
-        # "Sample Name": get_item(file_contents, sample_info_index[30:33], default=default, base=base),
-        # "Sample location": get_item(file_contents, sample_info_index[33:36], default=default, base=base),
-        # "Sample Material": get_item(file_contents, sample_info_index[36:39], default=default, base=base),
-        # "Experiment Type": get_item(file_contents, sample_info_index[39:42], default=default, base=base),
-        # "Sample Weight": get_item(file_contents, sample_info_index[42:45], default=default, base=base),
-        # "Step unit": get_item(file_contents, sample_info_index[45:48], default=default, base=base),
-        # "Heating time": get_item(file_contents, sample_info_index[48:51], default=default, base=base),
-        # "Instrument name": get_item(file_contents, sample_info_index[51:54], default=default, base=base),
-        # "Researcher": get_item(file_contents, sample_info_index[54:57], default=default, base=base),
-        # "Analyst": get_item(file_contents, sample_info_index[57:60], default=default, base=base),
-        # "Laboratory": get_item(file_contents, sample_info_index[60:63], default=default, base=base),
-        # "J value": get_item(file_contents, sample_info_index[63:66], default=default, base=base),
-        # "J value error": get_item(file_contents, sample_info_index[66:69], default=default, base=base),
-        # "Calc params": get_item(file_contents, sample_info_index[69:72], default=default, base=base),
-        # "Irra name": get_item(file_contents, sample_info_index[72:75], default=default, base=base),
-        # "Irra label": get_item(file_contents, sample_info_index[75:78], default=default, base=base),
-        # "Irra position H": get_item(file_contents, sample_info_index[78:81], default=default, base=base),
-        # "Irra position X": get_item(file_contents, sample_info_index[81:84], default=default, base=base),
-        # "Irra position Y": get_item(file_contents, sample_info_index[84:87], default=default, base=base),
-        # "Standard name": get_item(file_contents, sample_info_index[87:90], default=default, base=base),
-        # "Standard age": get_item(file_contents, sample_info_index[90:93], default=default, base=base),
-        # "Standard age error": get_item(file_contents, sample_info_index[93:96], default=default, base=base),
+        "MDF": get_item(file_contents, index_list[69:72], default=default, base=base),
+        "MDFSig": get_item(file_contents, index_list[72:75], default=default, base=base),
+        "CalcName": get_item(file_contents, index_list[75:78], default=default, base=base),
+        "IrraName": get_item(file_contents, index_list[78:81], default=default, base=base),
+        "IrraLabel": get_item(file_contents, index_list[81:84], default=default, base=base),
+        "IrraPosH": get_item(file_contents, index_list[84:87], default=default, base=base),
+        "IrraPosX": get_item(file_contents, index_list[87:90], default=default, base=base),
+        "IrraPosY": get_item(file_contents, index_list[90:93], default=default, base=base),
+        "StdName": get_item(file_contents, index_list[93:96], default=default, base=base),
+        "StdAge": get_item(file_contents, index_list[96:99], default=default, base=base),
+        "StdAgeSig": get_item(file_contents, index_list[99:102], default=default, base=base),
     })
     return sample_info

ararpy/smp/corr.py

@@ -266,7 +266,7 @@ def calc_degas_cl(sample: Sample):
             d6 = (a1 - a2 + (a3 + a4) / a7 - a5 / a7) / a7 / (1 - a6 / a7) ** 2
             d7 = -(a1 - a2) * (a6) / (a7 - a6) ** 2 - (a3 + a4 - a5) / (a7 - a6) ** 2
 
-            v1 = (a1 - a2 + (a3 + a4) / a7 - a5 / a7) / (1 - a6 / a7)
+            v1 = (a1 - a2 + (a3 + a4 - a5) / a7) / (1 - a6 / a7)
             s1 = (d1**2*s1**2 + d2**2*s2**2 + d3**2*s3**2 + d4**2*s4**2+ d5**2*s5**2 + d6**2*s6**2 + d7**2*s7**2) ** .5
             s2 = calc.err.mul((v1, s1), (vDecay[i], sDecay[i]))
             v1 = 0 if (ar36acl[0][i] - v1 < 0 or v1 < 0) and set_negative_zero[i] else v1

ararpy/smp/export.py

@@ -19,6 +19,7 @@ import pickle
 import traceback
 import pdf_maker as pm
 import numpy as np
+from math import log
 from decimal import Decimal
 
 from ..calc import arr, isochron, spectra, err
@@ -63,19 +64,19 @@ def get_cv_from_dict(data: dict, **kwargs):
         scale = (float(scale[0]), float(scale[1]), float(scale[2]), float(scale[3]))
         # create plot area based on axis scale
         plot_area = (kwargs.get("pt_left", 0.15), kwargs.get("pt_bottom", 0.15), kwargs.get("pt_width", 0.8), kwargs.get("pt_height", 0.8))
-        pt = cv.add_plot_area(name=f"PlotArea{i}", plot_area=plot_area, plot_scale=scale, show_frame=True)
+        pt = cv.add_plot_area(name=f"PlotArea{i}", plot_area=plot_area, plot_scale=scale, **data['xAxis'][i])
         for stick in data['xAxis'][i]['interval']:
             start = pt.scale_to_points(float(stick), scale[2])
             end = pt.scale_to_points(float(stick), scale[2])
             end = (end[0], end[1] - 5)
-            if pt.line(start=start, end=end, width=1, line_style="solid", y_clip=False, coordinate="pt", z_index=100):
+            if pt.line(start=start, end=end, width=data['xAxis'][i]['line_width'], line_style="solid", y_clip=False, coordinate="pt", z_index=100):
                 pt.text(x=start[0], y=end[1] - 15, text=f"{stick}", clip=False, size=int(data['xAxis'][i]['label_size']),
                         coordinate="pt", h_align="middle", z_index=150)
         for stick in data['yAxis'][i]['interval']:
             start = pt.scale_to_points(scale[0], float(stick))
             end = pt.scale_to_points(scale[0], float(stick))
             end = (end[0] - 5, end[1])
-            if pt.line(start=start, end=end, width=1, line_style="solid", x_clip=False, coordinate="pt", z_index=100):
+            if pt.line(start=start, end=end, width=data['xAxis'][i]['line_width'], line_style="solid", x_clip=False, coordinate="pt", z_index=100):
                 pt.text(x=end[0] - 5, y=end[1], text=f"{stick}", clip=False, size=int(data['yAxis'][i]['label_size']),
                         coordinate="pt", h_align="right", v_align="center", z_index=150)
         # axis titles
@@ -107,13 +108,13 @@ def get_cv_from_dict(data: dict, **kwargs):
                     start=data[index - 1], end=data[index], width=se.get('line_width', 1),
                     line_style=se.get('line_style', 'solid'), name=se['name'],
                     color=se.get('color', 'black'), clip=True, line_caps=se.get('line_caps', 'none'),
-                    z_index=se.get('z_index', 9))
+                    z_index=int(se.get('z_index', 9)))
         if 'scatter' in se['type'] and se['name'] != 'Text':
             for each in data:
                 pt.scatter(
                     each[0], each[1], fill_color=se.get('fill_color', 'black'), size=se.get('size', 5),
                     stroke_color=se.get('stroke_color', se.get('color', 'black')),
-                    z_index=se.get('z_index', 9)
+                    z_index=int(se.get('z_index', 9))
                 )
         if 'scatter' in se['type'] and se['name'] == 'Text' or 'text' in se['type']:
             for each in data:
@@ -181,6 +182,12 @@ def export_chart_to_pdf(cvs, file_name: str = "", file_path: str = "", **kwargs)
     -------
 
     """
+    # capital letter for label
+    def get_label(l="", num=0):
+        if num // 26 > 0:
+            return get_label(get_label(l, num // 26 - 1), num % 26)
+        return l + chr(num % 26 + 65)
+
     # write pdf
     show_label = kwargs.get('show_label', False)
     num_cols = kwargs.get('num_columns', 2)
@@ -195,21 +202,12 @@ def export_chart_to_pdf(cvs, file_name: str = "", file_path: str = "", **kwargs)
         file.del_obj(i)
     file.add_page(size=page_size[kwargs.get('page_size', 'a4').lower()], unit='mm')
     for page_index, page in enumerate(cvs):
-
-        # for x in range(0, 600, 10):
-        #     file.line(page_index, start=(x, 0), end=(x, 900), style='dash', color='grey')
-        # for y in range(0, 900, 10):
-        #     file.line(page_index, start=(0, y), end=(800, y), style='dash', color='grey')
-
-        numbers = ['A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'J', 'K', 'L', 'M', 'N']
-
         for cv_index, cv in enumerate(page):
             left = 10 + cv_index % num_cols * cv.width()
             top = 20 + cv_index // num_cols * cv.height()
             if show_label:
-                # file.text(page=page_index, x=left + 65, y=810 - top, text=numbers[cv_index])
                 cv._plot_areas[1].text(
-                    x=5, y=95, text=numbers[cv_index],
+                    x=5, y=95, text=f"({get_label(num=cv_index)})",
                     id=f'text-label-{page_index}-{cv_index}',
                     name=f'text-label-{page_index}-{cv_index}',
                     color='black', size=title_size, axis_index=1, h_align='middle', v_align='center'
@@ -327,28 +325,36 @@ def _get_additional_text_series(name):
 
 def _get_plot_data_age_spectra_series(smp: sample, **options):
     color = options.get('color', 'black')
-    sigma = options.get('sigma', 1)
+    sigma = options.get('sigma', smp.Info.preference['confidenceLevel'])
+    unit = options.get('unit', smp.Info.preference['ageUnit'])
     xAxis = options.get('xAxis', [{}])
     yAxis = options.get('yAxis', [{}])
+    line_width = options.get('line_width', 1)
+    line_caps = options.get('line_caps', 'square')
     series = []
+    # apparent spectra lines
     age = smp.ApparentAgeValues[2:4]
     ar = smp.DegasValues[20]
-    data = spectra.get_data(*age, ar, cumulative=False)
+    data = spectra.get_data(*age, ar, cumulative=False, sigma=sigma)
     series.append({
         'type': 'series.line', 'id': f'line-{smp.name()}-{get_random_digits()}', 'name': f'line-{smp.name()}-{get_random_digits()}',
-        'color': color, 'fill_color': color, 'line_width': 1, 'line_style': 'solid', 'z_index': 9,
-        'data': np.transpose([data[0], data[1]]).tolist(), 'line_caps': 'square',
+        'color': color, 'fill_color': color, 'line_width': line_width, 'line_style': 'solid', 'z_index': 19,
+        'data': np.transpose([data[0], data[1]]).tolist(),
+        'line_caps': line_caps,
+        # 'line_caps': 'none',
         'axis_index': 0,
     })
     series.append({
         'type': 'series.line', 'id': f'line-{smp.name()}-{get_random_digits()}', 'name': f'line-{smp.name()}-{get_random_digits()}',
-        'color': color, 'fill_color': color, 'line_width': 1, 'line_style': 'solid', 'z_index': 9,
-        'data': np.transpose([data[0], data[2]]).tolist(), 'line_caps': 'square',
+        'color': color, 'fill_color': color, 'line_width': line_width, 'line_style': 'solid', 'z_index': 19,
+        'data': np.transpose([data[0], data[2]]).tolist(),
+        'line_caps': line_caps,
+        # 'line_caps': 'none',
         'axis_index': 0,
     })
+    # initial values corrected spectra lines and rect
     text1 = smp.AgeSpectraPlot.text1
     text2 = smp.AgeSpectraPlot.text2
-    sigma = 1
     set_data = [spectraData2Sgima(smp.AgeSpectraPlot.set1.data, sigma),
                 spectraData2Sgima(smp.AgeSpectraPlot.set2.data, sigma)]
     for index, text in enumerate([text1, text2]):
@@ -358,56 +364,63 @@ def _get_plot_data_age_spectra_series(smp: sample, **options):
             plateau_data = np.array(set_data[index])
             # spectra_fill_data = list(map(lambda each: [data[0][each - 1], min(data[1][each], data[2][each]), data[0][each + 1] - data[0][each], sigma * abs(data[1][each] - data[2][each])], selection))
             spectra_fill_data = getRectFromSpectra(plateau_data)
+            print(spectra_fill_data)
             series.append({
                 'type': 'series.rect', 'id': f'rect-{smp.name()}-{get_random_digits()}', 'name': f'rect-{smp.name()}-{get_random_digits()}',
-                'color': ['red', 'blue'][index], 'fill_color': ['red', 'blue'][index], 'line_width': 0, 'line_style': 'solid', 'z_index': 9,
-                'data': spectra_fill_data, 'line_caps': 'square', 'fill': True,
+                'color': ['red', 'blue'][index], 'fill_color': ['red', 'blue'][index], 'line_width': 0, 'line_style': 'none', 'z_index': 9,
+                'data': spectra_fill_data, 'line_caps': 'none', 'fill': True,
                 'axis_index': 0,
             })
             series.append({
                 'type': 'series.line', 'id': f'line-{smp.name()}-{get_random_digits()}',
                 'name': f'line-{smp.name()}-{get_random_digits()}',
-                'color': ['red', 'blue'][index], 'line_width': 1, 'line_style': 'solid', 'z_index': 9,
-                'data': plateau_data[:,[0,1]].tolist(), 'line_caps': 'square',
+                'color': ['red', 'blue'][index], 'line_width': line_width, 'line_style': 'solid', 'z_index': 9,
+                'data': plateau_data[:,[0,1]].tolist(),
+                'line_caps': line_caps,
+                # 'line_caps': 'none',
                 'axis_index': 0,
             })
             series.append({
                 'type': 'series.line', 'id': f'line-{smp.name()}-{get_random_digits()}',
                 'name': f'line-{smp.name()}-{get_random_digits()}',
-                'color': ['red', 'blue'][index], 'line_width': 1, 'line_style': 'solid', 'z_index': 9,
-                'data': plateau_data[:,[0,2]].tolist(), 'line_caps': 'square',
+                'color': ['red', 'blue'][index], 'line_width': line_width, 'line_style': 'solid', 'z_index': 9,
+                'data': plateau_data[:,[0,2]].tolist(),
+                'line_caps': line_caps,
+                # 'line_caps': 'none',
                 'axis_index': 0,
             })
             series.append({
                 'type': 'text', 'id': f'text-{smp.name()}-{get_random_digits()}', 'name': f'text-{smp.name()}-{get_random_digits()}',
                 'color': ['red', 'blue'][index], 'fill_color': ['red', 'blue'][index],
-                'text': f"WMPA = {res['age']:.2f} ± {(sigma * res['s3']):.2f} Ma<r>"
+                'text': f"WMPA = {res['age']:.2f} ± {(sigma * res['s3']):.2f} {unit} ({sigma}s)<r>"
                         f"MSWD = {res['MSWD']:.2f}, n = {res['Num']:.0f}<r>"
-                        f"<sup>39</sup>Ar = {res['Ar39']:.2f}%",
+                        f"<sup>39</sup>Ar = {res['Ar39']*100:.2f}%",
                 'size': title_size,
-                'data': [text.pos],
-                'axis_index': 1,
+                # 'data': [text.pos],
+                'data': [[50, 50]],
+                'axis_index': 1, 'z_index': 29,
                 'h_align': "middle",
                 'v_align': "center",
             })
+    # TGA
     tga = smp.Info.results.age_spectra['TGA']
     if not np.isnan(tga['age']):
         series.append({
             'type': 'text', 'id': f'text-{smp.name()}-{get_random_digits()}', 'name': f'text-{smp.name()}-{get_random_digits()}',
             'color': color, 'fill_color': color,
-            'text': f"TGA = {tga['age']:.2f} ± {(sigma * tga['s3']):.2f} Ma",
+            'text': f"TGA = {tga['age']:.2f} ± {(sigma * tga['s3']):.2f} {unit} ({sigma}s)",
             'size': title_size,
             'data': [[50, 3]],
             'axis_index': 1,
             'h_align': "middle",
             'v_align': "bottom",
         })
-
+    # title
     series.append({
         'type': 'text', 'id': f'text-{smp.name()}-{get_random_digits()}',
         'name': f'text-{smp.name()}-{get_random_digits()}', 'color': color,
         'text': smp.name(), 'size': title_size,
-        'data': [[50, 95]], 'axis_index': 1,
+        'data': [[50, 95]], 'axis_index': 1, 'z_index': 29,
         'h_align': "middle", 'v_align': "center",
     })
 
@@ -416,50 +429,60 @@ def _get_plot_data_age_spectra_series(smp: sample, **options):
 
 def _get_plot_data_age_spectra_axis(smp: sample, **options):
     color = options.get('color', 'black')
+    line_width = options.get('line_width', 1)
     xAxis, yAxis = [], []
     age = smp.ApparentAgeValues[2:4]
     ar = smp.DegasValues[20]
     data = spectra.get_data(*age, ar, cumulative=False)
 
-    y_low_limit, y_up_limit, stick_num, inc = get_axis_scale(np.array(data[1:3]).flatten())
-    y_extent = [y_low_limit, y_up_limit]
-    y_interval = [y_low_limit + i * inc for i in range(stick_num + 1)]
-    y_interval.extend([y_up_limit * 10] * 3)
-    dp = int(np.log10(abs(y_up_limit - y_low_limit)) - 1)
-    dp = 0 if dp >= 0 else abs(dp)
-    print(f"{dp = }")
-    y_interval = [f"{i:.{dp}f}" for i in y_interval]
-
-    # y_extent = [4, 12]
-    # y_interval = ['0', '10', '20', '30', '40', '50', '60', '70', '80', '90', '100']
-    # y_interval = ['0', '4', '6', '8', '10', '12', '60', '70', '80', '90', '100']
+    # y_low_limit, y_up_limit, stick_num, inc = get_axis_scale(np.array(data[1:3]).flatten())
+    # y_extent = [y_low_limit, y_up_limit]
+    # y_interval = [y_low_limit + i * inc for i in range(stick_num + 1)]
+    # y_interval.extend([y_up_limit * 10] * 3)
+    # dp = int(np.log10(abs(y_up_limit - y_low_limit)) - 1)
+    # dp = 0 if dp >= 0 else abs(dp)
+    # y_interval = [f"{i:.{dp}f}" for i in y_interval]
+    try:
+        y_low_limit = float(smp.AgeSpectraPlot.yaxis.min)
+        y_up_limit = float(smp.AgeSpectraPlot.yaxis.max)
+        y_extent = [y_low_limit, y_up_limit]
+        stick_num = int(smp.AgeSpectraPlot.yaxis.split_number)
+        inc = float(smp.AgeSpectraPlot.yaxis.interval)
+        y_interval = [y_low_limit + i * inc for i in range(stick_num + 1)]
+        dp = int(np.log10(abs(y_up_limit - y_low_limit)) - 1)
+        dp = 0 if dp >= 0 else abs(dp)
+        y_interval = [f"{i:.{dp}f}" for i in y_interval]
+    except (BaseException, Exception):
+        y_extent = [0, 250]
+        y_interval = ['0', '50', '100', '150', '200', '250', '600', '700', '800', '900', '1000']
 
     xAxis.append({
         'extent': [0, 100], 'interval': [0, 20, 40, 60, 80, 100], 'id': 0, 'show_frame': True,
-        'title': 'Cumulative <sup>39</sup>Ar Released (%)', 'name_location': 'middle',
-        'line_width': 1, 'line_style': 'solid', 'z_index': 9, 'label_size': label_size, 'title_size': title_size,
+        'title': 'Cumulative <sup>39</sup>Ar released (%)', 'name_location': 'middle',
+        'line_width': line_width, 'line_style': 'solid', 'z_index': 9, 'label_size': label_size, 'title_size': title_size,
     })
     xAxis.append({
         'extent': [0, 100], 'interval': [], 'id': 1, 'show_frame': False,
         'title': '', 'name_location': 'middle',
-        'line_width': 1, 'line_style': 'solid', 'z_index': 0, 'label_size': label_size, 'title_size': title_size,
+        'line_width': 0, 'line_style': 'none', 'z_index': 0, 'label_size': label_size, 'title_size': title_size,
     })
     yAxis.append({
         'extent': y_extent, 'interval': y_interval,
         'id': 0, 'show_frame': True,
-        'title': 'Apparent Age (Ma)', 'name_location': 'middle',
-        'line_width': 1, 'line_style': 'solid', 'z_index': 9, 'label_size': label_size, 'title_size': title_size,
+        'title': 'Apparent age (Ma)', 'name_location': 'middle',
+        'line_width': line_width, 'line_style': 'solid', 'z_index': 9, 'label_size': label_size, 'title_size': title_size,
     })
     yAxis.append({
         'extent': [0, 100], 'interval': [], 'id': 1, 'show_frame': False,
         'title': '', 'name_location': 'middle',
-        'line_width': 1, 'line_style': 'solid', 'z_index': 0, 'label_size': label_size, 'title_size': title_size,
+        'line_width': 0, 'line_style': 'none', 'z_index': 0, 'label_size': label_size, 'title_size': title_size,
     })
     return xAxis, yAxis
 
 
 def _get_plot_data_inv_isochron_axis(smp: sample, **options):
     color = options.get('color', 'black')
+    line_width = options.get('line_width', 1)
     xAxis, yAxis = [], []
     data = np.array(smp.InvIsochronPlot.data)
     set1 = smp.InvIsochronPlot.set1.data
@@ -478,10 +501,12 @@ def _get_plot_data_inv_isochron_axis(smp: sample, **options):
         'interval': interval,
         'id': 0, 'show_frame': True, 'z_index': 9, 'label_size': label_size, 'title_size': title_size,
         'title': '<sup>39</sup>Ar<sub>K</sub> / <sup>40</sup>Ar<sup>*</sup>', 'name_location': 'middle',
+        'line_width': line_width,
     })
     xAxis.append({
         'extent': [0, 100], 'interval': [], 'id': 1, 'show_frame': False,
         'title': '', 'name_location': 'middle', 'z_index': 0, 'label_size': label_size, 'title_size': title_size,
+        'line_width': 0,
     })
 
     low_limit, up_limit, stick_num, inc = get_axis_scale(data[2])
@@ -1164,7 +1189,6 @@ class WritingWorkbook:
         ]:
             for index, _ in enumerate(['Set 1', 'Set 2']):
                 ar39[index] = res[index].get('Ar39', ar39[index] * 100) / 100
-                print(res[index])
                 try:
                     content.extend([
                         [(row + num_step, 0, row + 1 + num_step, 0), _, {'bold': 1, 'align': 'center'}],

{ararpy-0.1.25.dist-info → ararpy-0.1.27.dist-info}/METADATA

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

{ararpy-0.1.25.dist-info → ararpy-0.1.27.dist-info}/RECORD

@@ -1,9 +1,9 @@
-ararpy/__init__.py,sha256=IaTaNQCT51JP-Xn9mL40aPYSjWqEPbZAtHp0nkoS52c,6857
+ararpy/__init__.py,sha256=vsCWb5yzto-kGLLNpOzAb74KvfITyybMQ5RgEJJkmlU,6857
 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=jD1Fd0Cj3xc7NqgnG4cp3VWQWxUlV0qCtPBZZokDG8o,15246
-ararpy/calc/basic.py,sha256=uJCCUFaPd9zvfkggrdbFYSGLl2pt7UJ7ENgXanzHy68,4036
+ararpy/calc/basic.py,sha256=fVpdXmyqcIUU26nhU1bnI-YK4ZQ4Qqk24P3Yi9YwUwI,4187
 ararpy/calc/corr.py,sha256=QDitvLf6y5IBpHtVEJfOUVjWSlIP9qJ4Jhluk9eDEDE,18801
 ararpy/calc/err.py,sha256=63LtprqjemlIb1QGDst4Ggcv5KMSDHdlAIL-nyQs1eA,2691
 ararpy/calc/histogram.py,sha256=0GVbDdsjd91KQ1sa2B7NtZ4KGo0XpRIJapgIrzAwQUo,5777
@@ -35,18 +35,18 @@ ararpy/examples/sample-default.smp,sha256=YNkoQGgPrsL_fXS7ZHxfRtLQWekCDqT9czS6vB
 ararpy/files/__init__.py,sha256=l5B5ZQ01WdtvjjN0aMkyAFNgpwANdM_1I0tQbqnRuEY,69
 ararpy/files/arr_file.py,sha256=KqksGlEA6nmMQofTgi7v45flscQZVtefxaNCKrV3Am4,837
 ararpy/files/basic.py,sha256=nc7Hgo_qLSkdmtKzZmd5SQ8Jy0dhW46ly4gh-oisUDs,2095
-ararpy/files/calc_file.py,sha256=nqv_VfbOzz8ejcnjNoHJafGiYldhYSOLusUGtuZ8jR0,28207
+ararpy/files/calc_file.py,sha256=wbt-cQhPnFJ32QyofXvjCPPehmWSj4CHnidkrP7dFgk,28694
 ararpy/files/new_file.py,sha256=efblARIBROVLWS2w3-98BxLX5VZ8grRpiTkJFtf_rAk,214
-ararpy/files/raw_file.py,sha256=FScTqZzMPinNSCew86ETyDEcTjhbV-FQ3o4mWrbiBbU,21564
+ararpy/files/raw_file.py,sha256=PSJ0WwZHSIRrNDre69UwzpGevZQkI7utqhm9aKAm_Rs,20797
 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=PMvlSoToY8CJEORRrhEo4rYM2-VfSmjg0IuDpGFOGUM,24899
 ararpy/smp/calculation.py,sha256=LCFJWjLVLEKEQ5b7RFUIxsMahEzgLdodW4kCYXV5Z34,2919
 ararpy/smp/consts.py,sha256=XIdjdz8cYxspG2jMnoItdlUsxr3hKbNFJjMZJh1bpzw,393
-ararpy/smp/corr.py,sha256=L4YwCzu14Zot9tm9ijOgYzMbZEYEZGXArACrquIrfDg,26516
+ararpy/smp/corr.py,sha256=TzH7rymToJP0RCI0pQm-SrqnwAZ58HJbxahE8chv2wQ,26511
 ararpy/smp/diffusion_funcs.py,sha256=4-PMMIZWzjk2HOYYWNgSp4GmApygp1MmOxJ2g3xrqWc,175049
-ararpy/smp/export.py,sha256=s89L5B1aHoCeJIjcw6nXA6NtV0j_5XXaWnOETWnomCs,115043
+ararpy/smp/export.py,sha256=VvpqV_WxScbJ1bcmS6lDiBCa2bIn-TuumUOcDo78ItA,116201
 ararpy/smp/info.py,sha256=iKUELm-BuUduDlJKC1d8tKKNHbwwbNmhUg2pi6bcBvA,489
 ararpy/smp/initial.py,sha256=VSNku2fYD1xGkpHxMqnSJCAhSQawfsDCdqYCxtDJryQ,17196
 ararpy/smp/json.py,sha256=BTZCjVN0aj9epc700nwkYEYMKN2lHBYo-pLmtnz5oHY,2300
@@ -59,8 +59,8 @@ ararpy/thermo/__init__.py,sha256=6VBuqTRFl403PVqOuMkVrut0nKaQsAosBmfW91X1dMg,263
 ararpy/thermo/arrhenius.py,sha256=Ass1ichHfqIAtpv8eLlgrUc1UOb3Urh1qzr1E3gLB4U,233
 ararpy/thermo/atomic_level_random_walk.py,sha256=Q97zfe2h2RaxADkoBAqd0uEiP16BFOajrTmXHMkL2EQ,25502
 ararpy/thermo/basic.py,sha256=nBGHI9uK7VdJwThwBIOcKAzdnYqPyQseFoY6s4zKizk,11504
-ararpy-0.1.25.dist-info/licenses/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
-ararpy-0.1.25.dist-info/METADATA,sha256=_jTZMLzib2249ffidHh8OG0E3stzFUz9LXDlIRkQNzk,24516
-ararpy-0.1.25.dist-info/WHEEL,sha256=Nw36Djuh_5VDukK0H78QzOX-_FQEo6V37m3nkm96gtU,91
-ararpy-0.1.25.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
-ararpy-0.1.25.dist-info/RECORD,,
+ararpy-0.1.27.dist-info/licenses/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
+ararpy-0.1.27.dist-info/METADATA,sha256=ekOHHUdF4oyaBG5S7NJR7I4daFUFs-YzFCLo-Z3JXng,24516
+ararpy-0.1.27.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+ararpy-0.1.27.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
+ararpy-0.1.27.dist-info/RECORD,,

{ararpy-0.1.25.dist-info → ararpy-0.1.27.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (80.7.1)
+Generator: setuptools (80.9.0)
 Root-Is-Purelib: true
 Tag: py3-none-any