imap-processing 0.11.0__py3-none-any.whl → 0.13.0__py3-none-any.whl

imap_processing/hi/utils.py

@@ -4,7 +4,6 @@ import re
 from collections.abc import Sequence
 from dataclasses import dataclass
 from enum import IntEnum
-from numbers import Number
 from typing import Optional, Union
 
 import numpy as np
@@ -19,10 +18,12 @@ class HIAPID(IntEnum):
     H45_APP_NHK = 754
     H45_SCI_CNT = 769
     H45_SCI_DE = 770
+    H45_DIAG_FEE = 772
 
     H90_APP_NHK = 818
     H90_SCI_CNT = 833
     H90_SCI_DE = 834
+    H90_DIAG_FEE = 836
 
     @property
     def sensor(self) -> str:
@@ -100,7 +101,7 @@ def full_dataarray(
     attrs: dict,
     coords: Optional[dict[str, xr.DataArray]] = None,
     shape: Optional[Union[int, Sequence[int]]] = None,
-    fill_value: Optional[Number] = None,
+    fill_value: Optional[float] = None,
 ) -> xr.DataArray:
     """
     Generate an empty xarray.DataArray with appropriate attributes.
@@ -115,10 +116,14 @@ def full_dataarray(
         Variable name.
     attrs : dict
         CDF variable attributes. Usually retrieved from ImapCdfAttributes.
-    coords : dict
-        Coordinate variables for the Dataset.
-    shape : int or tuple
-        Shape of ndarray data array to instantiate in the xarray.DataArray.
+    coords : dict, optional
+        Coordinate variables for the Dataset. This function will extract the
+        sizes of each dimension defined by the attributes dictionary to determine
+        the size of the DataArray to be created.
+    shape : int or tuple, optional
+        Shape of ndarray data array to instantiate in the xarray.DataArray. If
+        shape is provided, the DataArray created will have this shape regardless
+        of whether coordinates are provided or not.
     fill_value : optional, float
         Override the fill value that the DataArray will be filled with. If not
         supplied, the "FILLVAL" value from `attrs` will be used.
@@ -152,9 +157,9 @@ def full_dataarray(
 
 def create_dataset_variables(
     variable_names: list[str],
-    variable_shape: Union[int, Sequence[int]],
+    variable_shape: Optional[Union[int, Sequence[int]]] = None,
     coords: Optional[dict[str, xr.DataArray]] = None,
-    fill_value: Optional[Number] = None,
+    fill_value: Optional[float] = None,
     att_manager_lookup_str: str = "{0}",
 ) -> dict[str, xr.DataArray]:
     """
@@ -166,10 +171,14 @@ def create_dataset_variables(
     ----------
     variable_names : list[str]
         List of variable names to create.
-    variable_shape : int or sequence of int
-        Shape of the new variables data ndarray.
-    coords : dict
-        Coordinate variables for the Dataset.
+    variable_shape : int or sequence of int, optional
+        Shape of the new variables data ndarray. If not provided the shape will
+        attempt to be derived from the coords dictionary.
+    coords : dict, optional
+        Coordinate variables for the Dataset. If `variable_shape` is not provided
+        the dataset variables created will use this dictionary along with variable
+        attributes from the CdfAttributeManager to determine the shapes of the
+        dataset variables created.
     fill_value : optional, number
         Value to fill the new variables data arrays with. If not supplied,
         the fill value is pulled from the CDF variable attributes "FILLVAL"
@@ -200,3 +209,40 @@ def create_dataset_variables(
             var, attrs, shape=variable_shape, coords=coords, fill_value=fill_value
         )
     return new_variables
+
+
+class CoincidenceBitmap(IntEnum):
+    """IntEnum class for coincidence type bitmap values."""
+
+    A = 2**3
+    B = 2**2
+    C1 = 2**1
+    C2 = 2**0
+
+    @staticmethod
+    def detector_hit_str_to_int(detector_hit_str: str) -> int:
+        """
+        Convert a detector hit string to a coincidence type integer value.
+
+        A detector hit string is a string containing all detectors that were hit
+        for a direct event. Possible detectors include: [A, B, C1, C2]. Converting
+        the detector hit string to a coincidence type integer value involves
+        summing the coincidence bitmap value for each detector hit. e.g. "AC1C2"
+        results in 2**3 + 2**1 + 2**0 = 11.
+
+        Parameters
+        ----------
+        detector_hit_str : str
+            The string containing the set of detectors hit.
+            e.g. "AC1C2".
+
+        Returns
+        -------
+        coincidence_type : int
+            The integer value of the coincidence type.
+        """
+        # Join all detector names with a pipe for use with regex
+        pattern = r"|".join(c.name for c in CoincidenceBitmap)
+        matches = re.findall(pattern, detector_hit_str)
+        # Sum the integer value assigned to the detector name for each match
+        return sum(CoincidenceBitmap[m] for m in matches)

imap_processing/hit/ancillary/imap_hit_l1b-to-l2-sector-dt0-factors_20250219_v002.csv

@@ -0,0 +1,81 @@
+Species, Sector, Lower Energy (MeV), Upper Energy (MeV), Delta E (MeV), Geometry Factor (cm2 sr), Efficiency, b
+H      ,0,1.8,4,2.2,0.41686,1,0
+H      ,1,1.8,4,2.2,0.49253,1,0
+H      ,2,1.8,4,2.2,0.45156,1,0
+H      ,3,1.8,4,2.2,0.35599,1,0
+H      ,4,1.8,4,2.2,0.33312,1,0
+H      ,5,1.8,4,2.2,0.47444,1,0
+H      ,6,1.8,4,2.2,0.49253,1,0
+H      ,7,1.8,4,2.2,0.41686,1,0
+H      ,0,4,6,2,0.42008,1,0
+H      ,1,4,6,2,0.49746,1,0
+H      ,2,4,6,2,0.45515,1,0
+H      ,3,4,6,2,0.35921,1,0
+H      ,4,4,6,2,0.33619,1,0
+H      ,5,4,6,2,0.47817,1,0
+H      ,6,4,6,2,0.49746,1,0
+H      ,7,4,6,2,0.42008,1,0
+H      ,0,6,10,4,0.39177,1,0
+H      ,1,6,10,4,0.47095,1,0
+H      ,2,6,10,4,0.42764,1,0
+H      ,3,6,10,4,0.33171,1,0
+H      ,4,6,10,4,0.31091,1,0
+H      ,5,6,10,4,0.44844,1,0
+H      ,6,6,10,4,0.47095,1,0
+H      ,7,6,10,4,0.39177,1,0
+He4    ,0,4,6,2,0.41995,1,0
+He4    ,1,4,6,2,0.49964,1,0
+He4    ,2,4,6,2,0.45649,1,0
+He4    ,3,4,6,2,0.3584,1,0
+He4    ,4,4,6,2,0.33563,1,0
+He4    ,5,4,6,2,0.47925,1,0
+He4    ,6,4,6,2,0.49964,1,0
+He4    ,7,4,6,2,0.41995,1,0
+He4    ,0,6,12,6,0.42404,1,0
+He4    ,1,6,12,6,0.5045,1,0
+He4    ,2,6,12,6,0.45897,1,0
+He4    ,3,6,12,6,0.36098,1,0
+He4    ,4,6,12,6,0.33793,1,0
+He4    ,5,6,12,6,0.48202,1,0
+He4    ,6,6,12,6,0.5045,1,0
+He4    ,7,6,12,6,0.42404,1,0
+CNO    ,0,4,6,2,0.41976,1,0
+CNO    ,1,4,6,2,0.49594,1,0
+CNO    ,2,4,6,2,0.45356,1,0
+CNO    ,3,4,6,2,0.35698,1,0
+CNO    ,4,4,6,2,0.33409,1,0
+CNO    ,5,4,6,2,0.47645,1,0
+CNO    ,6,4,6,2,0.49594,1,0
+CNO    ,7,4,6,2,0.41976,1,0
+CNO    ,0,6,12,6,0.42026,1,0
+CNO    ,1,6,12,6,0.49871,1,0
+CNO    ,2,6,12,6,0.45567,1,0
+CNO    ,3,6,12,6,0.35809,1,0
+CNO    ,4,6,12,6,0.3351,1,0
+CNO    ,5,6,12,6,0.47866,1,0
+CNO    ,6,6,12,6,0.49871,1,0
+CNO    ,7,6,12,6,0.42026,1,0
+NeMgSi ,0,4,6,2,0.37544,1,0
+NeMgSi ,1,4,6,2,0.43535,1,0
+NeMgSi ,2,4,6,2,0.40251,1,0
+NeMgSi ,3,4,6,2,0.32308,1,0
+NeMgSi ,4,4,6,2,0.30221,1,0
+NeMgSi ,5,4,6,2,0.42339,1,0
+NeMgSi ,6,4,6,2,0.43535,1,0
+NeMgSi ,7,4,6,2,0.37544,1,0
+NeMgSi ,0,6,12,6,0.42046,1,0
+NeMgSi ,1,6,12,6,0.49714,1,0
+NeMgSi ,2,6,12,6,0.45432,1,0
+NeMgSi ,3,6,12,6,0.35764,1,0
+NeMgSi ,4,6,12,6,0.33475,1,0
+NeMgSi ,5,6,12,6,0.47721,1,0
+NeMgSi ,6,6,12,6,0.49714,1,0
+NeMgSi ,7,6,12,6,0.42046,1,0
+Fe     ,0,4,12,8,0.2752,1,0
+Fe     ,1,4,12,8,0.31996,1,0
+Fe     ,2,4,12,8,0.29577,1,0
+Fe     ,3,4,12,8,0.23669,1,0
+Fe     ,4,4,12,8,0.2213,1,0
+Fe     ,5,4,12,8,0.31118,1,0
+Fe     ,6,4,12,8,0.31996,1,0
+Fe     ,7,4,12,8,0.2752,1,0

imap_processing/hit/ancillary/imap_hit_l1b-to-l2-standard-dt0-factors_20250219_v002.csv

@@ -0,0 +1,205 @@
+Species,Lower Energy (MeV),Upper Energy (MeV),Delta E (MeV),Geometry Factor (cm2 sr),Efficiency,b
+H,1.8,2.2,0.4,3.4146,1,0
+H,2.2,2.7,0.5,3.44142,1,0
+H,2.7,3.2,0.5,3.44784,1,0
+H,3.2,3.6,0.4,3.44966,1,0
+H,3.6,4,0.4,3.4596,1,0
+H,4,4.5,0.5,3.45668,1,0
+H,4.5,5,0.5,3.46674,1,0
+H,5,6,1,3.47076,1,0
+H,6,8,2,3.46876,1,0
+H,8,10,2,2.81904,1,0
+H,10,12,2,1.10134,1,0
+H,12,15,3,0.20036,1,0
+He3,2.2,2.7,0.5,3.41332,1,0
+He3,2.7,3.2,0.5,3.45214,1,0
+He3,3.2,3.6,0.4,3.45662,1,0
+He3,3.6,4,0.4,3.4603,1,0
+He3,4,4.5,0.5,3.457,1,0
+He3,4.5,5,0.5,3.46028,1,0
+He3,5,6,1,3.46604,1,0
+He3,6,8,2,3.48172,1,0
+He3,8,10,2,3.49966,1,0
+He3,10,12,2,3.52024,1,0
+He3,12,15,3,3.52832,1,0
+He4,1.8,2.2,0.4,3.44744,1,0
+He4,2.2,2.7,0.5,3.45572,1,0
+He4,2.7,3.2,0.5,3.46184,1,0
+He4,3.2,3.6,0.4,3.45344,1,0
+He4,3.6,4,0.4,3.46182,1,0
+He4,4,4.5,0.5,3.47216,1,0
+He4,4.5,5,0.5,3.46908,1,0
+He4,5,6,1,3.47572,1,0
+He4,6,8,2,3.48828,1,0
+He4,8,10,2,3.50768,1,0
+He4,10,12,2,3.52508,1,0
+He4,12,15,3,3.53526,1,0
+He,2.2,2.7,0.5,3.45528,1,0
+He,2.7,3.2,0.5,3.46172,1,0
+He,3.2,3.6,0.4,3.4535,1,0
+He,3.6,4,0.4,3.4618,1,0
+He,4,4.5,0.5,3.472,1,0
+He,4.5,5,0.5,3.46898,1,0
+He,5,6,1,3.47566,1,0
+He,6,8,2,3.48824,1,0
+He,8,10,2,3.50758,1,0
+He,10,12,2,3.525,1,0
+He,12,15,3,3.53518,1,0
+C,3.2,3.6,0.4,3.41806,1,0
+C,3.6,4,0.4,3.43468,1,0
+C,4,4.5,0.5,3.4472,1,0
+C,4.5,5,0.5,3.4677,1,0
+C,5,6,1,3.46138,1,0
+C,6,8,2,3.46458,1,0
+C,8,10,2,3.47978,1,0
+C,10,12,2,3.474,1,0
+C,12,15,3,3.50882,1,0
+C,15,21,6,3.52714,1,0
+C,21,27,6,3.55204,1,0
+C,27,33,6,3.87152,1,0
+N,3.2,3.6,0.4,2.8743,1,0
+N,3.6,4,0.4,3.4377,1,0
+N,4,4.5,0.5,3.44746,1,0
+N,4.5,5,0.5,3.45346,1,0
+N,5,6,1,3.45678,1,0
+N,6,8,2,3.46598,1,0
+N,8,10,2,3.47534,1,0
+N,10,12,2,3.48836,1,0
+N,12,15,3,3.50548,1,0
+N,15,21,6,3.53182,1,0
+N,21,27,6,3.55572,1,0
+N,27,33,6,3.65516,1,0
+O,3.2,3.6,0.4,1.32412,1,0
+O,3.6,4,0.4,3.40832,1,0
+O,4,4.5,0.5,3.44954,1,0
+O,4.5,5,0.5,3.45086,1,0
+O,5,6,1,3.45518,1,0
+O,6,8,2,3.46514,1,0
+O,8,10,2,3.47316,1,0
+O,10,12,2,3.47588,1,0
+O,12,15,3,3.49978,1,0
+O,15,21,6,3.52392,1,0
+O,21,27,6,3.56288,1,0
+O,27,33,6,3.59158,1,0
+Ne,3.2,3.6,0.4,0.00184,1,0
+Ne,3.6,4,0.4,1.8455,1,0
+Ne,4,4.5,0.5,3.425,1,0
+Ne,4.5,5,0.5,3.43902,1,0
+Ne,5,6,1,3.44606,1,0
+Ne,6,8,2,3.4629,1,0
+Ne,8,10,2,3.47594,1,0
+Ne,10,12,2,3.47526,1,0
+Ne,12,15,3,3.49242,1,0
+Ne,15,21,6,3.51476,1,0
+Ne,21,27,6,3.55948,1,0
+Ne,27,33,6,3.59022,1,0
+Ne,33,40,7,3.70364,1,0
+Na,6,8,2,3.45928,1,0
+Na,8,10,2,3.46536,1,0
+Na,10,12,2,3.47552,1,0
+Na,12,15,3,3.4934,1,0
+Na,15,21,6,3.51534,1,0
+Na,21,27,6,3.55662,1,0
+Na,27,33,6,3.59688,1,0
+Na,33,40,7,3.71056,1,0
+Mg,3.2,3.6,0.4,0.00056,1,0
+Mg,3.6,4,0.4,0.01146,1,0
+Mg,4,4.5,0.5,2.6573,1,0
+Mg,4.5,5,0.5,3.44594,1,0
+Mg,5,6,1,3.44902,1,0
+Mg,6,8,2,3.45872,1,0
+Mg,8,10,2,3.46468,1,0
+Mg,10,12,2,3.48568,1,0
+Mg,12,15,3,3.49012,1,0
+Mg,15,21,6,3.51206,1,0
+Mg,21,27,6,3.55654,1,0
+Mg,27,33,6,3.5903,1,0
+Mg,33,40,7,3.6286,1,0
+Mg,40,52,12,4.05832,1,0
+Al,6,8,2,3.45308,1,0
+Al,8,10,2,3.46306,1,0
+Al,10,12,2,3.47248,1,0
+Al,12,15,3,3.486,1,0
+Al,15,21,6,3.51144,1,0
+Al,21,27,6,3.55256,1,0
+Al,27,33,6,3.58784,1,0
+Al,33,40,7,3.63428,1,0
+Al,40,52,12,4.05398,1,0
+Si,3.2,3.6,0.4,0.00058,1,0
+Si,3.6,4,0.4,0.00166,1,0
+Si,4,4.5,0.5,1.16254,1,0
+Si,4.5,5,0.5,3.38976,1,0
+Si,5,6,1,3.44992,1,0
+Si,6,8,2,3.45184,1,0
+Si,8,10,2,3.46212,1,0
+Si,10,12,2,3.4661,1,0
+Si,12,15,3,3.48014,1,0
+Si,15,21,6,3.51188,1,0
+Si,21,27,6,3.54904,1,0
+Si,27,33,6,3.59342,1,0
+Si,33,40,7,3.64074,1,0
+Si,40,52,12,3.92812,1,0
+S,3.6,4,0.4,0.00148,1,0
+S,4,4.5,0.5,0.13366,1,0
+S,4.5,5,0.5,3.03138,1,0
+S,5,6,1,3.44574,1,0
+S,6,8,2,3.4527,1,0
+S,8,10,2,3.46542,1,0
+S,10,12,2,3.47918,1,0
+S,12,15,3,3.48282,1,0
+S,15,21,6,3.50736,1,0
+S,21,27,6,3.54404,1,0
+S,27,33,6,3.58826,1,0
+S,33,40,7,3.63246,1,0
+S,40,52,12,3.78708,1,0
+Ar,3.6,4,0.4,0.0024,1,0
+Ar,4,4.5,0.5,1.03066,1,0
+Ar,4.5,5,0.5,3.34924,1,0
+Ar,5,6,1,3.44938,1,0
+Ar,6,8,2,3.46096,1,0
+Ar,8,10,2,3.46986,1,0
+Ar,10,12,2,3.4794,1,0
+Ar,12,15,3,3.48722,1,0
+Ar,15,21,6,3.50784,1,0
+Ar,21,27,6,3.54998,1,0
+Ar,27,33,6,3.59038,1,0
+Ar,33,40,7,3.63132,1,0
+Ar,40,52,12,3.78276,1,0
+Ca,3.6,4,0.4,0.00084,1,0
+Ca,4,4.5,0.5,0.00246,1,0
+Ca,4.5,5,0.5,0.96168,1,0
+Ca,5,6,1,3.3616,1,0
+Ca,6,8,2,3.45146,1,0
+Ca,8,10,2,3.46034,1,0
+Ca,10,12,2,3.47002,1,0
+Ca,12,15,3,3.47898,1,0
+Ca,15,21,6,3.5036,1,0
+Ca,21,27,6,3.53976,1,0
+Ca,27,33,6,3.58474,1,0
+Ca,33,40,7,3.62258,1,0
+Ca,40,52,12,3.67846,1,0
+Fe,2.7,3.2,0.5,0.00018,1,0
+Fe,3.2,3.6,0.4,0.00024,1,0
+Fe,3.6,4,0.4,0.00096,1,0
+Fe,4,4.5,0.5,0.00246,1,0
+Fe,4.5,5,0.5,0.3899,1,0
+Fe,5,6,1,3.26502,1,0
+Fe,6,8,2,3.45292,1,0
+Fe,8,10,2,3.45522,1,0
+Fe,10,12,2,3.47062,1,0
+Fe,12,15,3,3.47424,1,0
+Fe,15,21,6,3.4955,1,0
+Fe,21,27,6,3.52856,1,0
+Fe,27,33,6,3.5562,1,0
+Fe,33,40,7,3.5978,1,0
+Fe,40,52,12,3.64532,1,0
+Fe,52,70,18,3.96388,1,0
+Ni,8,10,2,3.45838,1,0
+Ni,10,12,2,3.47046,1,0
+Ni,12,15,3,3.47826,1,0
+Ni,15,21,6,3.49778,1,0
+Ni,21,27,6,3.52796,1,0
+Ni,27,33,6,3.56026,1,0
+Ni,33,40,7,3.59852,1,0
+Ni,40,52,12,3.64296,1,0
+Ni,52,70,18,3.87798,1,0

imap_processing/hit/ancillary/imap_hit_l1b-to-l2-standard-dt1-factors_20250219_v002.csv

@@ -0,0 +1,205 @@
+Species,Lower Energy (MeV),Upper Energy (MeV),Delta E (MeV),Geometry Factor (cm2 sr),Efficiency,b
+H,1.8,2.2,0.4,0.68292,1,0
+H,2.2,2.7,0.5,0.688284,1,0
+H,2.7,3.2,0.5,0.689568,1,0
+H,3.2,3.6,0.4,0.689932,1,0
+H,3.6,4,0.4,0.69192,1,0
+H,4,4.5,0.5,0.691336,1,0
+H,4.5,5,0.5,0.693348,1,0
+H,5,6,1,0.694152,1,0
+H,6,8,2,0.693752,1,0
+H,8,10,2,0.563808,1,0
+H,10,12,2,0.220268,1,0
+H,12,15,3,0.040072,1,0
+He3,2.2,2.7,0.5,0.682664,1,0
+He3,2.7,3.2,0.5,0.690428,1,0
+He3,3.2,3.6,0.4,0.691324,1,0
+He3,3.6,4,0.4,0.69206,1,0
+He3,4,4.5,0.5,0.6914,1,0
+He3,4.5,5,0.5,0.692056,1,0
+He3,5,6,1,0.693208,1,0
+He3,6,8,2,0.696344,1,0
+He3,8,10,2,0.699932,1,0
+He3,10,12,2,0.704048,1,0
+He3,12,15,3,0.705664,1,0
+He4,1.8,2.2,0.4,0.689488,1,0
+He4,2.2,2.7,0.5,0.691144,1,0
+He4,2.7,3.2,0.5,0.692368,1,0
+He4,3.2,3.6,0.4,0.690688,1,0
+He4,3.6,4,0.4,0.692364,1,0
+He4,4,4.5,0.5,0.694432,1,0
+He4,4.5,5,0.5,0.693816,1,0
+He4,5,6,1,0.695144,1,0
+He4,6,8,2,0.697656,1,0
+He4,8,10,2,0.701536,1,0
+He4,10,12,2,0.705016,1,0
+He4,12,15,3,0.707052,1,0
+He,2.2,2.7,0.5,0.691056,1,0
+He,2.7,3.2,0.5,0.692344,1,0
+He,3.2,3.6,0.4,0.6907,1,0
+He,3.6,4,0.4,0.69236,1,0
+He,4,4.5,0.5,0.6944,1,0
+He,4.5,5,0.5,0.693796,1,0
+He,5,6,1,0.695132,1,0
+He,6,8,2,0.697648,1,0
+He,8,10,2,0.701516,1,0
+He,10,12,2,0.705,1,0
+He,12,15,3,0.707036,1,0
+C,3.2,3.6,0.4,0.683612,1,0
+C,3.6,4,0.4,0.686936,1,0
+C,4,4.5,0.5,0.68944,1,0
+C,4.5,5,0.5,0.69354,1,0
+C,5,6,1,0.692276,1,0
+C,6,8,2,0.692916,1,0
+C,8,10,2,0.695956,1,0
+C,10,12,2,0.6948,1,0
+C,12,15,3,0.701764,1,0
+C,15,21,6,0.705428,1,0
+C,21,27,6,0.710408,1,0
+C,27,33,6,0.774304,1,0
+N,3.2,3.6,0.4,0.57486,1,0
+N,3.6,4,0.4,0.68754,1,0
+N,4,4.5,0.5,0.689492,1,0
+N,4.5,5,0.5,0.690692,1,0
+N,5,6,1,0.691356,1,0
+N,6,8,2,0.693196,1,0
+N,8,10,2,0.695068,1,0
+N,10,12,2,0.697672,1,0
+N,12,15,3,0.701096,1,0
+N,15,21,6,0.706364,1,0
+N,21,27,6,0.711144,1,0
+N,27,33,6,0.731032,1,0
+O,3.2,3.6,0.4,0.264824,1,0
+O,3.6,4,0.4,0.681664,1,0
+O,4,4.5,0.5,0.689908,1,0
+O,4.5,5,0.5,0.690172,1,0
+O,5,6,1,0.691036,1,0
+O,6,8,2,0.693028,1,0
+O,8,10,2,0.694632,1,0
+O,10,12,2,0.695176,1,0
+O,12,15,3,0.699956,1,0
+O,15,21,6,0.704784,1,0
+O,21,27,6,0.712576,1,0
+O,27,33,6,0.718316,1,0
+Ne,3.2,3.6,0.4,0.000368,1,0
+Ne,3.6,4,0.4,0.3691,1,0
+Ne,4,4.5,0.5,0.685,1,0
+Ne,4.5,5,0.5,0.687804,1,0
+Ne,5,6,1,0.689212,1,0
+Ne,6,8,2,0.69258,1,0
+Ne,8,10,2,0.695188,1,0
+Ne,10,12,2,0.695052,1,0
+Ne,12,15,3,0.698484,1,0
+Ne,15,21,6,0.702952,1,0
+Ne,21,27,6,0.711896,1,0
+Ne,27,33,6,0.718044,1,0
+Ne,33,40,7,0.740728,1,0
+Na,6,8,2,0.691856,1,0
+Na,8,10,2,0.693072,1,0
+Na,10,12,2,0.695104,1,0
+Na,12,15,3,0.69868,1,0
+Na,15,21,6,0.703068,1,0
+Na,21,27,6,0.711324,1,0
+Na,27,33,6,0.719376,1,0
+Na,33,40,7,0.742112,1,0
+Mg,3.2,3.6,0.4,0.000112,1,0
+Mg,3.6,4,0.4,0.002292,1,0
+Mg,4,4.5,0.5,0.53146,1,0
+Mg,4.5,5,0.5,0.689188,1,0
+Mg,5,6,1,0.689804,1,0
+Mg,6,8,2,0.691744,1,0
+Mg,8,10,2,0.692936,1,0
+Mg,10,12,2,0.697136,1,0
+Mg,12,15,3,0.698024,1,0
+Mg,15,21,6,0.702412,1,0
+Mg,21,27,6,0.711308,1,0
+Mg,27,33,6,0.71806,1,0
+Mg,33,40,7,0.72572,1,0
+Mg,40,52,12,0.811664,1,0
+Al,6,8,2,0.690616,1,0
+Al,8,10,2,0.692612,1,0
+Al,10,12,2,0.694496,1,0
+Al,12,15,3,0.6972,1,0
+Al,15,21,6,0.702288,1,0
+Al,21,27,6,0.710512,1,0
+Al,27,33,6,0.717568,1,0
+Al,33,40,7,0.726856,1,0
+Al,40,52,12,0.810796,1,0
+Si,3.2,3.6,0.4,0.000116,1,0
+Si,3.6,4,0.4,0.000332,1,0
+Si,4,4.5,0.5,0.232508,1,0
+Si,4.5,5,0.5,0.677952,1,0
+Si,5,6,1,0.689984,1,0
+Si,6,8,2,0.690368,1,0
+Si,8,10,2,0.692424,1,0
+Si,10,12,2,0.69322,1,0
+Si,12,15,3,0.696028,1,0
+Si,15,21,6,0.702376,1,0
+Si,21,27,6,0.709808,1,0
+Si,27,33,6,0.718684,1,0
+Si,33,40,7,0.728148,1,0
+Si,40,52,12,0.785624,1,0
+S,3.6,4,0.4,0.000296,1,0
+S,4,4.5,0.5,0.026732,1,0
+S,4.5,5,0.5,0.606276,1,0
+S,5,6,1,0.689148,1,0
+S,6,8,2,0.69054,1,0
+S,8,10,2,0.693084,1,0
+S,10,12,2,0.695836,1,0
+S,12,15,3,0.696564,1,0
+S,15,21,6,0.701472,1,0
+S,21,27,6,0.708808,1,0
+S,27,33,6,0.717652,1,0
+S,33,40,7,0.726492,1,0
+S,40,52,12,0.757416,1,0
+Ar,3.6,4,0.4,0.00048,1,0
+Ar,4,4.5,0.5,0.206132,1,0
+Ar,4.5,5,0.5,0.669848,1,0
+Ar,5,6,1,0.689876,1,0
+Ar,6,8,2,0.692192,1,0
+Ar,8,10,2,0.693972,1,0
+Ar,10,12,2,0.69588,1,0
+Ar,12,15,3,0.697444,1,0
+Ar,15,21,6,0.701568,1,0
+Ar,21,27,6,0.709996,1,0
+Ar,27,33,6,0.718076,1,0
+Ar,33,40,7,0.726264,1,0
+Ar,40,52,12,0.756552,1,0
+Ca,3.6,4,0.4,0.000168,1,0
+Ca,4,4.5,0.5,0.000492,1,0
+Ca,4.5,5,0.5,0.192336,1,0
+Ca,5,6,1,0.67232,1,0
+Ca,6,8,2,0.690292,1,0
+Ca,8,10,2,0.692068,1,0
+Ca,10,12,2,0.694004,1,0
+Ca,12,15,3,0.695796,1,0
+Ca,15,21,6,0.70072,1,0
+Ca,21,27,6,0.707952,1,0
+Ca,27,33,6,0.716948,1,0
+Ca,33,40,7,0.724516,1,0
+Ca,40,52,12,0.735692,1,0
+Fe,2.7,3.2,0.5,0.000036,1,0
+Fe,3.2,3.6,0.4,0.000048,1,0
+Fe,3.6,4,0.4,0.000192,1,0
+Fe,4,4.5,0.5,0.000492,1,0
+Fe,4.5,5,0.5,0.07798,1,0
+Fe,5,6,1,0.653004,1,0
+Fe,6,8,2,0.690584,1,0
+Fe,8,10,2,0.691044,1,0
+Fe,10,12,2,0.694124,1,0
+Fe,12,15,3,0.694848,1,0
+Fe,15,21,6,0.6991,1,0
+Fe,21,27,6,0.705712,1,0
+Fe,27,33,6,0.71124,1,0
+Fe,33,40,7,0.71956,1,0
+Fe,40,52,12,0.729064,1,0
+Fe,52,70,18,0.792776,1,0
+Ni,8,10,2,0.691676,1,0
+Ni,10,12,2,0.694092,1,0
+Ni,12,15,3,0.695652,1,0
+Ni,15,21,6,0.699556,1,0
+Ni,21,27,6,0.705592,1,0
+Ni,27,33,6,0.712052,1,0
+Ni,33,40,7,0.719704,1,0
+Ni,40,52,12,0.728592,1,0
+Ni,52,70,18,0.775596,1,0