AeroViz 0.1.21__py3-none-any.whl

AeroViz/data/hysplit_example_data.txt

@@ -0,0 +1,101 @@
+     1     1
+    AWRF    24     1    25     0     0
+     4 BACKWARD OMEGA   
+    24     2    27    16   22.630  120.346    10.0
+    24     2    27    16   22.630  120.346   100.0
+    24     2    27    16   22.630  120.346   500.0
+    24     2    27    16   22.630  120.346  1000.0
+     1 PRESSURE
+     1     1    24     2    27    16     0   808     0.0   22.630  120.346     10.0   1018.1
+     2     1    24     2    27    16     0   808     0.0   22.630  120.346    100.0   1007.1
+     3     1    24     2    27    16     0   808     0.0   22.630  120.346    500.0    961.2
+     4     1    24     2    27    16     0   808     0.0   22.630  120.346   1000.0    906.2
+     1     1    24     2    27    15     0   807    -1.0   22.700  120.337     13.4   1016.4
+     2     1    24     2    27    15     0   807    -1.0   22.779  120.325    115.8   1004.4
+     3     1    24     2    27    15     0   807    -1.0   22.615  120.313    471.0    964.8
+     4     1    24     2    27    15     0   807    -1.0   22.535  120.469   1035.2    902.9
+     1     1    24     2    27    14     0   806    -2.0   22.795  120.323     24.3   1015.4
+     2     1    24     2    27    14     0   806    -2.0   23.026  120.322    128.8   1001.6
+     3     1    24     2    27    14     0   806    -2.0   22.606  120.315    471.5    964.7
+     4     1    24     2    27    14     0   806    -2.0   22.484  120.626   1025.3    893.1
+     1     1    24     2    27    13     0   805    -3.0   22.936  120.308     56.6   1010.4
+     2     1    24     2    27    13     0   805    -3.0   23.275  120.355    142.2   1000.7
+     3     1    24     2    27    13     0   805    -3.0   22.612  120.341    492.9    961.7
+     4     1    24     2    27    13     0   805    -3.0   22.490  120.729    641.9    860.0
+     1     1    24     2    27    12     0   804    -4.0   23.133  120.283     74.9   1010.2
+     2     1    24     2    27    12     0   804    -4.0   23.491  120.413    187.7    994.2
+     3     1    24     2    27    12     0   804    -4.0   22.670  120.382    482.0    960.7
+     4     1    24     2    27    12     0   804    -4.0   22.538  120.774    830.3    801.8
+     1     1    24     2    27    11     0   803    -5.0   23.363  120.276     94.7   1007.0
+     2     1    24     2    27    11     0   803    -5.0   23.702  120.464    202.2    992.8
+     3     1    24     2    27    11     0   803    -5.0   22.708  120.455    327.6    977.5
+     4     1    24     2    27    11     0   803    -5.0   22.576  120.736    629.4    831.5
+     1     1    24     2    27    10     0   802    -6.0   23.611  120.281    123.4   1003.6
+     2     1    24     2    27    10     0   802    -6.0   23.924  120.507    235.0    988.9
+     3     1    24     2    27    10     0   802    -6.0   22.691  120.521     95.4   1004.1
+     4     1    24     2    27    10     0   802    -6.0   22.631  120.727    385.5    834.4
+     1     1    24     2    27     9     0   801    -7.0   23.910  120.338    185.0    996.9
+     2     1    24     2    27     9     0   801    -7.0   24.136  120.570    314.7    974.2
+     3     1    24     2    27     9     0   801    -7.0   22.702  120.585    125.1    999.7
+     4     1    24     2    27     9     0   801    -7.0   22.669  120.749    200.9    817.0
+     1     1    24     2    27     8     0   800    -8.0   24.287  120.521    321.6    982.5
+     2     1    24     2    27     8     0   800    -8.0   24.195  120.654    531.7    942.2
+     3     1    24     2    27     8     0   800    -8.0   22.707  120.608    170.4    993.9
+     4     1    24     2    27     8     0   800    -8.0   22.689  120.756    313.0    822.5
+     1     1    24     2    27     7     0   799    -9.0   24.612  120.810    428.8    969.2
+     2     1    24     2    27     7     0   799    -9.0   24.193  120.736    697.7    917.4
+     3     1    24     2    27     7     0   799    -9.0   22.719  120.610    229.7    986.5
+     4     1    24     2    27     7     0   799    -9.0   22.689  120.751    248.2    833.6
+     1     1    24     2    27     6     0   798   -10.0   24.859  121.146    538.3    932.4
+     2     1    24     2    27     6     0   798   -10.0   24.145  120.757    664.1    919.1
+     3     1    24     2    27     6     0   798   -10.0   22.738  120.605    305.9    978.9
+     4     1    24     2    27     6     0   798   -10.0   22.711  120.761      0.0    888.3
+     1     1    24     2    27     5     0   797   -11.0   24.984  121.430    392.9    975.2
+     2     1    24     2    27     5     0   797   -11.0   24.074  120.749    654.4    923.5
+     3     1    24     2    27     5     0   797   -11.0   22.731  120.603    495.6    958.3
+     4     1    24     2    27     5     0   797   -11.0   22.700  120.785    175.9    829.1
+     1     1    24     2    27     4     0   796   -12.0   25.158  121.638    503.0    954.6
+     2     1    24     2    27     4     0   796   -12.0   23.989  120.737    641.7    927.8
+     3     1    24     2    27     4     0   796   -12.0   22.693  120.592    648.1    943.6
+     4     1    24     2    27     4     0   796   -12.0   22.685  120.810    145.4    853.0
+     1     1    24     2    27     3     0   795   -13.0   25.451  121.825    386.1    978.5
+     2     1    24     2    27     3     0   795   -13.0   23.906  120.739    645.3    929.6
+     3     1    24     2    27     3     0   795   -13.0   22.657  120.583    682.1    940.8
+     4     1    24     2    27     3     0   795   -13.0   22.655  120.810      8.7    903.0
+     1     1    24     2    27     2     0   794   -14.0   25.789  121.989    450.1    971.0
+     2     1    24     2    27     2     0   794   -14.0   23.854  120.722    640.7    924.2
+     3     1    24     2    27     2     0   794   -14.0   22.618  120.552    686.5    941.2
+     4     1    24     2    27     2     0   794   -14.0   22.628  120.817    147.7    890.5
+     1     1    24     2    27     1     0   793   -15.0   26.127  122.136    554.7    959.4
+     2     1    24     2    27     1     0   793   -15.0   23.815  120.684    526.0    936.4
+     3     1    24     2    27     1     0   793   -15.0   22.537  120.538    763.0    933.1
+     4     1    24     2    27     1     0   793   -15.0   22.630  120.825    156.0    898.7
+     1     1    24     2    27     0     0   792   -16.0   26.425  122.289    711.1    940.7
+     2     1    24     2    27     0     0   792   -16.0   23.793  120.665    527.9    942.8
+     3     1    24     2    27     0     0   792   -16.0   22.498  120.642    901.2    897.5
+     4     1    24     2    27     0     0   792   -16.0   22.633  120.822    130.8    898.3
+     1     1    24     2    26    23     0   791   -17.0   26.705  122.454    872.6    921.4
+     2     1    24     2    26    23     0   791   -17.0   23.798  120.671    524.1    940.3
+     3     1    24     2    26    23     0   791   -17.0   22.514  120.689    734.8    863.9
+     4     1    24     2    26    23     0   791   -17.0   22.628  120.825     80.9    904.3
+     2     1    24     2    26    22     0   790   -18.0   23.816  120.688    493.1    939.0
+     3     1    24     2    26    22     0   790   -18.0   22.559  120.685    771.9    865.4
+     4     1    24     2    26    22     0   790   -18.0   22.623  120.829     56.0    907.8
+     2     1    24     2    26    21     0   789   -19.0   23.836  120.711    427.5    946.0
+     3     1    24     2    26    21     0   789   -19.0   22.557  120.666    660.4    899.0
+     4     1    24     2    26    21     0   789   -19.0   22.623  120.833     82.4    907.7
+     2     1    24     2    26    20     0   788   -20.0   23.865  120.751    489.8    941.5
+     3     1    24     2    26    20     0   788   -20.0   22.525  120.624    588.1    942.1
+     4     1    24     2    26    20     0   788   -20.0   22.636  120.848     98.3    914.3
+     2     1    24     2    26    19     0   787   -21.0   23.874  120.791    693.5    909.2
+     3     1    24     2    26    19     0   787   -21.0   22.492  120.608    624.9    943.0
+     4     1    24     2    26    19     0   787   -21.0   22.667  120.865    193.7    891.2
+     2     1    24     2    26    18     0   786   -22.0   23.874  120.793    743.4    902.1
+     3     1    24     2    26    18     0   786   -22.0   22.482  120.601    703.4    938.1
+     4     1    24     2    26    18     0   786   -22.0   22.747  120.984    287.0    906.5
+     2     1    24     2    26    17     0   785   -23.0   23.885  120.793    766.4    902.4
+     3     1    24     2    26    17     0   785   -23.0   22.471  120.582    842.0    922.5
+     4     1    24     2    26    17     0   785   -23.0   22.867  121.070    462.0    920.8
+     2     1    24     2    26    16     0   784   -24.0   23.895  120.778    859.2    898.6
+     3     1    24     2    26    16     0   784   -24.0   22.446  120.618   1090.1    893.0
+     4     1    24     2    26    16     0   784   -24.0   22.973  121.135    790.6    906.9

AeroViz/dataProcess/Chemistry/__init__.py

@@ -0,0 +1,149 @@
+from ..core import Writer, run_process
+
+__all__ = ['Chemistry']
+
+
+class Chemistry(Writer):
+    # Reconstruction
+    @run_process('Chemistry - reconstruction basic', 'reconstrc_basic')
+    def ReConstrc_basic(self, *df_chem, df_ref=None, df_water=None, df_density=None, nam_lst=None):
+        from ._mass_volume import _basic
+
+        if nam_lst is None:
+            nam_lst = ['NH4+', 'SO42-', 'NO3-', 'Fe', 'Na+', 'OC', 'EC']
+
+        out = _basic(df_chem, df_ref, df_water, df_density, nam_lst=nam_lst)
+
+        return self, out
+
+    # Partition (Gas-Particle Partitioning Ratios)
+    @run_process('Chemistry - Partition', 'partition')
+    def Partition(self, df_data):
+        """
+        Calculate gas-particle partitioning ratios (SOR, NOR, NTR, epsilon).
+
+        Parameters
+        ----------
+        df_data : DataFrame
+            Data containing particle and gas concentrations (μg/m³).
+            Required: 'temp' column (temperature in Celsius)
+            Optional species (at least one pair needed):
+            - SO42-, SO2 : for SOR
+            - NO3-, NO2 : for NOR
+            - NH4+, NH3 : for NTR
+            - Cl-, HCl : for chloride partitioning
+
+        Returns
+        -------
+        DataFrame
+            Partitioning ratios: SOR, NOR, NOR_2, NTR, epls_* columns.
+
+        Notes
+        -----
+        - SOR > 0.1 indicates secondary sulfate formation
+        - Values near 1.0: particle phase dominant
+        - Values near 0.0: gas phase dominant
+        """
+        from ._calculate import partition_ratios
+
+        out = partition_ratios(df_data)
+
+        return self, out
+
+    # ISOROPIA
+    @run_process('Chemistry - ISOROPIA', 'isoropia')
+    def ISOROPIA(self, *df_chem, nam_lst=None):
+        from ._isoropia import _basic
+
+        if nam_lst is None:
+            nam_lst = ['Na+', 'SO42-', 'NH4+', 'NO3-', 'Cl-', 'Ca2+',
+                       'K+', 'Mg2+', 'NH3', 'HNO3', 'HCl', 'RH', 'temp']
+
+        if self.path_out is None:
+            raise ValueError('Please Input "path_out" !!')
+
+        out = _basic(df_chem, self.path_out, nam_lst=nam_lst)
+
+        return self, out
+
+    # OCEC
+    @run_process('Chemistry - OC/EC basic', 'ocec_basic')
+    def OCEC_basic(self, df_lcres, df_mass=None, ocec_ratio=None, ocec_ratio_month=1, hr_lim=200,
+                   least_square_range=(0.1, 2.5, 0.1), WISOC_OC_range=(0.2, 0.7, 0.01), ):
+        from ._ocec import _basic
+
+        out = _basic(df_lcres, df_mass, ocec_ratio, ocec_ratio_month, hr_lim, least_square_range, WISOC_OC_range)
+
+        return self, out
+
+    # Volume Average Mixing
+    @run_process('Chemistry - VAM', 'vam')
+    def volume_average_mixing(self, df_volume, df_alwc=None):
+        """
+        Calculate volume-average refractive index and growth factor.
+
+        Parameters
+        ----------
+        df_volume : DataFrame
+            Volume data with AS_volume, AN_volume, etc. and total_dry columns.
+        df_alwc : DataFrame, optional
+            Aerosol liquid water content data.
+
+        Returns
+        -------
+        DataFrame
+            RI data with n_dry, k_dry, n_amb, k_amb, gRH columns.
+        """
+        from ._calculate import volume_average_mixing
+
+        out = volume_average_mixing(df_volume, df_alwc)
+
+        return self, out
+
+    # Kappa calculation
+    @run_process('Chemistry - kappa', 'kappa')
+    def kappa(self, df_data, diameter=0.5):
+        """
+        Calculate hygroscopicity parameter kappa.
+
+        Parameters
+        ----------
+        df_data : DataFrame
+            Data containing gRH, AT (temperature in C), and RH (%).
+        diameter : float, default=0.5
+            Particle dry diameter in micrometers.
+
+        Returns
+        -------
+        DataFrame
+            Kappa values.
+        """
+        from ._calculate import kappa_calculate
+
+        out = kappa_calculate(df_data, diameter)
+
+        return self, out
+
+    # gRH calculation
+    @run_process('Chemistry - gRH', 'gRH')
+    def gRH(self, df_volume, df_alwc):
+        """
+        Calculate hygroscopic growth factor.
+
+        Parameters
+        ----------
+        df_volume : DataFrame
+            Volume data with 'total_dry' column.
+        df_alwc : DataFrame
+            Aerosol liquid water content with 'ALWC' column.
+
+        Returns
+        -------
+        DataFrame
+            Growth factor data.
+        """
+        from ._calculate import gRH_calculate
+
+        out = gRH_calculate(df_volume, df_alwc)
+
+        return self, out