AeroViz 0.1.2__py3-none-any.whl → 0.1.3b0__py3-none-any.whl

AeroViz/rawDataReader/script/AE33.py

@@ -4,28 +4,28 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'AE33'
+    nam = 'AE33'
 
-	def _raw_reader(self, _file):
-		_df = read_table(_file, parse_dates={'time': [0, 1]}, index_col='time',
-						 delimiter=r'\s+', skiprows=5, usecols=range(67))
-		_df.columns = _df.columns.str.strip(';')
+    def _raw_reader(self, _file):
+        _df = read_table(_file, parse_dates={'time': [0, 1]}, index_col='time',
+                         delimiter=r'\s+', skiprows=5, usecols=range(67))
+        _df.columns = _df.columns.str.strip(';')
 
-		# remove data without Status=0, 128 (Not much filter tape), 256 (Not much filter tape)
-		if not self._oth_set.get('ignore_err', False):
-			_df = _df.where((_df['Status'] != 0) | (_df['Status'] != 128) | (_df['Status'] != 256)).copy()
+        # remove data without Status=0, 128 (Not much filter tape), 256 (Not much filter tape)
+        if not self._oth_set.get('ignore_err', False):
+            _df = _df.where((_df['Status'] != 0) | (_df['Status'] != 128) | (_df['Status'] != 256)).copy()
 
-		return _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'Status']]
+        return _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'Status']]
 
-	def _QC(self, _df):
-		# remove negative value
-		_df = _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7']].mask((_df < 0).copy())
+    def _QC(self, _df):
+        # remove negative value
+        _df = _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7']].mask((_df < 0).copy())
 
-		# QC data in 5 min
-		def _QC_func(df):
-			_df_ave, _df_std = df.mean(), df.std()
-			_df_lowb, _df_highb = df < (_df_ave - _df_std * 1.5), df > (_df_ave + _df_std * 1.5)
+        # QC data in 5 min
+        def _QC_func(df):
+            _df_ave, _df_std = df.mean(), df.std()
+            _df_lowb, _df_highb = df < (_df_ave - _df_std * 1.5), df > (_df_ave + _df_std * 1.5)
 
-			return df.mask(_df_lowb | _df_highb).copy()
+            return df.mask(_df_lowb | _df_highb).copy()
 
-		return _df.resample('5min').apply(_QC_func).resample('1h').mean()
+        return _df.resample('5min').apply(_QC_func).resample('1h').mean()

AeroViz/rawDataReader/script/AE43.py

@@ -4,31 +4,31 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'AE43'
+    nam = 'AE43'
 
-	def _raw_reader(self, _file):
-		_df = read_csv(_file, parse_dates={'time': ['StartTime']}, index_col='time')
-		_df_id = _df['SetupID'].iloc[-1]
+    def _raw_reader(self, _file):
+        _df = read_csv(_file, parse_dates={'time': ['StartTime']}, index_col='time')
+        _df_id = _df['SetupID'].iloc[-1]
 
-		# get last SetupID data
-		_df = _df.groupby('SetupID').get_group(_df_id)[
-			['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'Status']].copy()
+        # get last SetupID data
+        _df = _df.groupby('SetupID').get_group(_df_id)[
+            ['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'Status']].copy()
 
-		# remove data without Status=0
-		_df = _df.where(_df['Status'] == 0).copy()
+        # remove data without Status=0
+        _df = _df.where(_df['Status'] == 0).copy()
 
-		return _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7']]
+        return _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7']]
 
-	# QC data
-	def _QC(self, _df):
-		# remove negative value
-		_df = _df.mask((_df < 0).copy())
+    # QC data
+    def _QC(self, _df):
+        # remove negative value
+        _df = _df.mask((_df < 0).copy())
 
-		# QC data in 5 min
-		def _QC_func(df):
-			_df_ave, _df_std = df.mean(), df.std()
-			_df_lowb, _df_highb = df < (_df_ave - _df_std * 1.5), df > (_df_ave + _df_std * 1.5)
+        # QC data in 5 min
+        def _QC_func(df):
+            _df_ave, _df_std = df.mean(), df.std()
+            _df_lowb, _df_highb = df < (_df_ave - _df_std * 1.5), df > (_df_ave + _df_std * 1.5)
 
-			return df.mask(_df_lowb | _df_highb).copy()
+            return df.mask(_df_lowb | _df_highb).copy()
 
-		return _df.resample('5min').apply(_QC_func).resample('1h').mean()
+        return _df.resample('5min').apply(_QC_func).resample('1h').mean()

AeroViz/rawDataReader/script/APS_3321.py

@@ -5,43 +5,43 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'APS_3321'
+    nam = 'APS_3321'
 
-	def _raw_reader(self, _file):
-		with open(_file, 'r', encoding='utf-8', errors='ignore') as f:
-			_df = read_table(f, skiprows=6, parse_dates={'Time': ['Date', 'Start Time']}).set_index('Time')
-			_key = list(_df.keys()[3:54])  ## 542 ~ 1981
+    def _raw_reader(self, _file):
+        with open(_file, 'r', encoding='utf-8', errors='ignore') as f:
+            _df = read_table(f, skiprows=6, parse_dates={'Time': ['Date', 'Start Time']}).set_index('Time')
+            _key = list(_df.keys()[3:54])  ## 542 ~ 1981
 
-			## create new keys
-			_newkey = {}
-			for _k in _key:
-				_newkey[_k] = float(_k).__round__(4)
-			# _newkey['Mode(m)'] = 'mode'
+            ## create new keys
+            _newkey = {}
+            for _k in _key:
+                _newkey[_k] = float(_k).__round__(4)
+            # _newkey['Mode(m)'] = 'mode'
 
-			## get new dataframe
-			_df = _df[_newkey.keys()].rename(_newkey, axis=1)
-			# _df['total'] = _df[list(_newkey.values())[:-1]].sum(axis=1)*(n.diff(n.log(_df.keys()[:-1].to_numpy(float))).mean()).copy()
+            ## get new dataframe
+            _df = _df[_newkey.keys()].rename(_newkey, axis=1)
+            # _df['total'] = _df[list(_newkey.values())[:-1]].sum(axis=1)*(n.diff(n.log(_df.keys()[:-1].to_numpy(float))).mean()).copy()
 
-			_df_idx = to_datetime(_df.index, errors='coerce')
+            _df_idx = to_datetime(_df.index, errors='coerce')
 
-		return _df.set_index(_df_idx).loc[_df_idx.dropna()]
+        return _df.set_index(_df_idx).loc[_df_idx.dropna()]
 
-	## QC data
-	def _QC(self, _df):
-		## mask out the data size lower than 7
-		_df['total'] = _df.sum(axis=1, min_count=1) * (n.diff(n.log(_df.keys().to_numpy(float)))).mean()
-		_df_size = _df['total'].dropna().resample('1h').size().resample(_df.index.freq).ffill()
-		_df = _df.mask(_df_size < 7)
+    ## QC data
+    def _QC(self, _df):
+        ## mask out the data size lower than 7
+        _df['total'] = _df.sum(axis=1, min_count=1) * (n.diff(n.log(_df.keys().to_numpy(float)))).mean()
+        _df_size = _df['total'].dropna().resample('1h').size().resample(_df.index.freq).ffill()
+        _df = _df.mask(_df_size < 7)
 
-		## remove total conc. lower than 700
-		_df = _df.mask(_df['total'] > 700)
+        ## remove total conc. lower than 700
+        _df = _df.mask(_df['total'] > 700)
 
-		# not confirmed
-		"""
-		## remove the bin over 4000 nm which num. conc. larger than 1
-		# _df_remv_ky = _df.keys()[:-2][_df.keys()[:-2]>=4.]
+        # not confirmed
+        """
+        ## remove the bin over 4000 nm which num. conc. larger than 1
+        # _df_remv_ky = _df.keys()[:-2][_df.keys()[:-2]>=4.]
 
-		# _df_1hr[_df_remv_ky] = _df_1hr[_df_remv_ky].copy().mask(_df_1hr[_df_remv_ky]>1.)
-		# """
+        # _df_1hr[_df_remv_ky] = _df_1hr[_df_remv_ky].copy().mask(_df_1hr[_df_remv_ky]>1.)
+        # """
 
-		return _df[_df.keys()[:-1]]
+        return _df[_df.keys()[:-1]]

AeroViz/rawDataReader/script/Aurora.py

@@ -4,35 +4,35 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'Aurora'
+    nam = 'Aurora'
 
-	def _raw_reader(self, _file):
-		with (_file).open('r', encoding='utf-8-sig', errors='ignore') as f:
-			_df = read_csv(f, low_memory=False, index_col=0)
+    def _raw_reader(self, _file):
+        with (_file).open('r', encoding='utf-8-sig', errors='ignore') as f:
+            _df = read_csv(f, low_memory=False, index_col=0)
 
-			_df.index = to_datetime(_df.index, errors='coerce', format=self._oth_set.get('date_format') or 'mixed')
-			_df.index.name = 'time'
+            _df.index = to_datetime(_df.index, errors='coerce', format=self._oth_set.get('date_format') or 'mixed')
+            _df.index.name = 'time'
 
-			_df.columns = _df.keys().str.strip(' ')
+            _df.columns = _df.keys().str.strip(' ')
 
-			_df = _df.loc[
-				_df.index.dropna(), ['0°σspB', '0°σspG', '0°σspR', '90°σspB', '90°σspG', '90°σspR', 'RH']].copy()
-			_df.columns = ['B', 'G', 'R', 'BB', 'BG', 'BR', 'RH']
+            _df = _df.loc[
+                _df.index.dropna(), ['0°σspB', '0°σspG', '0°σspR', '90°σspB', '90°σspG', '90°σspR', 'RH']].copy()
+            _df.columns = ['B', 'G', 'R', 'BB', 'BG', 'BR', 'RH']
 
-		return _df
+        return _df
 
-	## QC data
-	def _QC(self, _df):
-		## remove negative value
-		_df = _df.mask((_df <= 0).copy())
+    ## QC data
+    def _QC(self, _df):
+        ## remove negative value
+        _df = _df.mask((_df <= 0).copy())
 
-		## call by _QC function
-		## QC data in 1 hr
-		def _QC_func(_df_1hr):
-			_df_ave = _df_1hr.mean()
-			_df_std = _df_1hr.std()
-			_df_lowb, _df_highb = _df_1hr < (_df_ave - _df_std * 1.5), _df_1hr > (_df_ave + _df_std * 1.5)
+        ## call by _QC function
+        ## QC data in 1 hr
+        def _QC_func(_df_1hr):
+            _df_ave = _df_1hr.mean()
+            _df_std = _df_1hr.std()
+            _df_lowb, _df_highb = _df_1hr < (_df_ave - _df_std * 1.5), _df_1hr > (_df_ave + _df_std * 1.5)
 
-			return _df_1hr.mask(_df_lowb | _df_highb).copy()
+            return _df_1hr.mask(_df_lowb | _df_highb).copy()
 
-		return _df.resample('1h', group_keys=False).apply(_QC_func)
+        return _df.resample('1h', group_keys=False).apply(_QC_func)

AeroViz/rawDataReader/script/BC1054.py

@@ -4,43 +4,43 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'BC1054'
-
-	def _raw_reader(self, _file):
-		with open(_file, 'r', encoding='utf-8', errors='ignore') as f:
-			_df = read_csv(f, parse_dates=['Time'], index_col='Time')
-
-			_df = _df.rename(columns={
-				'BC1(ng/m3)': 'BC1',
-				'BC2(ng/m3)': 'BC2',
-				'BC3(ng/m3)': 'BC3',
-				'BC4(ng/m3)': 'BC4',
-				'BC5(ng/m3)': 'BC5',
-				'BC6(ng/m3)': 'BC6',
-				'BC7(ng/m3)': 'BC7',
-				'BC8(ng/m3)': 'BC8',
-				'BC9(ng/m3)': 'BC9',
-				'BC10(ng/m3)': 'BC10'
-			})
-
-			# remove data without Status=32 (Automatic Tape Advance), 65536 (Tape Move)
-			# if not self._oth_set.get('ignore_err', False):
-			#     _df = _df.where((_df['Status'] != 32) | (_df['Status'] != 65536)).copy()
-
-			return _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'BC8', 'BC9', 'BC10', 'Status']]
-
-	# QC data
-	def _QC(self, _df):
-		# remove negative value
-		_df = _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'BC8', 'BC9', 'BC10']].mask((_df < 0).copy())
-
-		# call by _QC function
-		# QC data in 1 hr
-		def _QC_func(_df_1hr):
-			_df_ave = _df_1hr.mean()
-			_df_std = _df_1hr.std()
-			_df_lowb, _df_highb = _df_1hr < (_df_ave - _df_std * 1.5), _df_1hr > (_df_ave + _df_std * 1.5)
-
-			return _df_1hr.mask(_df_lowb | _df_highb).copy()
-
-		return _df.resample('1h', group_keys=False).apply(_QC_func).resample('5min').mean()
+    nam = 'BC1054'
+
+    def _raw_reader(self, _file):
+        with open(_file, 'r', encoding='utf-8', errors='ignore') as f:
+            _df = read_csv(f, parse_dates=['Time'], index_col='Time')
+
+            _df = _df.rename(columns={
+                'BC1(ng/m3)': 'BC1',
+                'BC2(ng/m3)': 'BC2',
+                'BC3(ng/m3)': 'BC3',
+                'BC4(ng/m3)': 'BC4',
+                'BC5(ng/m3)': 'BC5',
+                'BC6(ng/m3)': 'BC6',
+                'BC7(ng/m3)': 'BC7',
+                'BC8(ng/m3)': 'BC8',
+                'BC9(ng/m3)': 'BC9',
+                'BC10(ng/m3)': 'BC10'
+            })
+
+            # remove data without Status=32 (Automatic Tape Advance), 65536 (Tape Move)
+            # if not self._oth_set.get('ignore_err', False):
+            #     _df = _df.where((_df['Status'] != 32) | (_df['Status'] != 65536)).copy()
+
+            return _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'BC8', 'BC9', 'BC10', 'Status']]
+
+    # QC data
+    def _QC(self, _df):
+        # remove negative value
+        _df = _df[['BC1', 'BC2', 'BC3', 'BC4', 'BC5', 'BC6', 'BC7', 'BC8', 'BC9', 'BC10']].mask((_df < 0).copy())
+
+        # call by _QC function
+        # QC data in 1 hr
+        def _QC_func(_df_1hr):
+            _df_ave = _df_1hr.mean()
+            _df_std = _df_1hr.std()
+            _df_lowb, _df_highb = _df_1hr < (_df_ave - _df_std * 1.5), _df_1hr > (_df_ave + _df_std * 1.5)
+
+            return _df_1hr.mask(_df_lowb | _df_highb).copy()
+
+        return _df.resample('1h', group_keys=False).apply(_QC_func).resample('5min').mean()

AeroViz/rawDataReader/script/EPA_vertical.py

@@ -4,15 +4,15 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'EPA_vertical'
+    nam = 'EPA_vertical'
 
-	def _raw_reader(self, _file):
-		with _file.open('r', encoding='big5', errors='ignore') as f:
-			_df = read_csv(f, names=['time', 'station', 'comp', 'data', None], skiprows=1, na_values=['-'],
-						   parse_dates=['time'], index_col='time')
-			_df['data'] = to_numeric(_df['data'], errors='coerce')
+    def _raw_reader(self, _file):
+        with _file.open('r', encoding='big5', errors='ignore') as f:
+            _df = read_csv(f, names=['time', 'station', 'comp', 'data', None], skiprows=1, na_values=['-'],
+                           parse_dates=['time'], index_col='time')
+            _df['data'] = to_numeric(_df['data'], errors='coerce')
 
-			_df_piv = _df.pivot_table(values='data', columns='comp', index='time')
-			_df_piv.index.name = 'time'
+            _df_piv = _df.pivot_table(values='data', columns='comp', index='time')
+            _df_piv.index.name = 'time'
 
-		return _df_piv
+        return _df_piv

AeroViz/rawDataReader/script/GRIMM.py

@@ -4,32 +4,32 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'GRIMM'
+    nam = 'GRIMM'
 
-	def _raw_reader(self, _file):
+    def _raw_reader(self, _file):
 
-		_df = read_csv(_file, header=233, delimiter='\t', index_col=0, parse_dates=[0], encoding='ISO-8859-1',
-					   dayfirst=True).rename_axis("Time")
-		_df.index = to_datetime(_df.index, format="%d/%m/%Y %H:%M:%S", dayfirst=True)
+        _df = read_csv(_file, header=233, delimiter='\t', index_col=0, parse_dates=[0], encoding='ISO-8859-1',
+                       dayfirst=True).rename_axis("Time")
+        _df.index = to_datetime(_df.index, format="%d/%m/%Y %H:%M:%S", dayfirst=True)
 
-		if _file.name.startswith("A407ST"):
-			_df.drop(_df.columns[0:11].tolist() + _df.columns[128:].tolist(), axis=1, inplace=True)
-		else:
-			_df.drop(_df.columns[0:11].tolist() + _df.columns[-5:].tolist(), axis=1, inplace=True)
+        if _file.name.startswith("A407ST"):
+            _df.drop(_df.columns[0:11].tolist() + _df.columns[128:].tolist(), axis=1, inplace=True)
+        else:
+            _df.drop(_df.columns[0:11].tolist() + _df.columns[-5:].tolist(), axis=1, inplace=True)
 
-		if _df.empty:
-			print(_file, "is empty")
-			return None
+        if _df.empty:
+            print(_file, "is empty")
+            return None
 
-		return _df / 0.035
+        return _df / 0.035
 
-	def _QC(self, _df):
-		# QC data in 1 hr
-		def _QC_func(_df_1hr):
-			_df_ave = _df_1hr.mean()
-			_df_std = _df_1hr.std()
-			_df_lowb, _df_highb = _df_1hr < (_df_ave - _df_std * 1.5), _df_1hr > (_df_ave + _df_std * 1.5)
+    def _QC(self, _df):
+        # QC data in 1 hr
+        def _QC_func(_df_1hr):
+            _df_ave = _df_1hr.mean()
+            _df_std = _df_1hr.std()
+            _df_lowb, _df_highb = _df_1hr < (_df_ave - _df_std * 1.5), _df_1hr > (_df_ave + _df_std * 1.5)
 
-			return _df_1hr.mask(_df_lowb | _df_highb).copy()
+            return _df_1hr.mask(_df_lowb | _df_highb).copy()
 
-		return _df.resample('5min').apply(_QC_func).resample('1h').mean()
+        return _df.resample('5min').apply(_QC_func).resample('1h').mean()

AeroViz/rawDataReader/script/IGAC_TH.py

@@ -8,97 +8,97 @@ from AeroViz.rawDataReader.core import AbstractReader
 
 
 class Reader(AbstractReader):
-	nam = 'IGAC_TH'
+    nam = 'IGAC_TH'
 
-	def _raw_reader(self, _file):
+    def _raw_reader(self, _file):
 
-		self.meta['freq'] = self._oth_set.get('data_freq') or self.meta['freq']
+        self.meta['freq'] = self._oth_set.get('data_freq') or self.meta['freq']
 
-		with (_file).open('r', encoding='utf-8-sig', errors='ignore') as f:
-			_df = read_csv(f, low_memory=False, index_col=0)
+        with (_file).open('r', encoding='utf-8-sig', errors='ignore') as f:
+            _df = read_csv(f, low_memory=False, index_col=0)
 
-			_df.index = to_datetime(_df.index, errors='coerce', format=self._oth_set.get('date_format') or 'mixed')
-			_df.index.name = 'time'
+            _df.index = to_datetime(_df.index, errors='coerce', format=self._oth_set.get('date_format') or 'mixed')
+            _df.index.name = 'time'
 
-			_df.columns = _df.keys().str.strip(' ')
+            _df.columns = _df.keys().str.strip(' ')
 
-			_df = _df.loc[_df.index.dropna()].copy()
+            _df = _df.loc[_df.index.dropna()].copy()
 
-		return _df.loc[~_df.index.duplicated()]
+        return _df.loc[~_df.index.duplicated()]
 
-	## QC data
-	def _QC(self, _df):
+    ## QC data
+    def _QC(self, _df):
 
-		## QC parameter, function (MDL SE LE)
-		_mdl = {
-			'Na+': 0.05,
-			'NH4+': 0.05,
-			'K+': 0.05,
-			'Mg2+': 0.05,
-			'Ca2+': 0.05,
-			'Cl-': 0.05,
-			'NO2-': 0.05,
-			'NO3-': 0.05,
-			'SO42-': 0.05,
-		}
+        ## QC parameter, function (MDL SE LE)
+        _mdl = {
+            'Na+': 0.05,
+            'NH4+': 0.05,
+            'K+': 0.05,
+            'Mg2+': 0.05,
+            'Ca2+': 0.05,
+            'Cl-': 0.05,
+            'NO2-': 0.05,
+            'NO3-': 0.05,
+            'SO42-': 0.05,
+        }
 
-		def _se_le(_df_, _log=False):
-			_df_ = np.log10(_df_) if _log else _df_
+        def _se_le(_df_, _log=False):
+            _df_ = np.log10(_df_) if _log else _df_
 
-			_df_qua = _df_.quantile([.25, .75])
-			_df_q1, _df_q3 = _df_qua.loc[.25].copy(), _df_qua.loc[.75].copy()
-			_df_iqr = _df_q3 - _df_q1
+            _df_qua = _df_.quantile([.25, .75])
+            _df_q1, _df_q3 = _df_qua.loc[.25].copy(), _df_qua.loc[.75].copy()
+            _df_iqr = _df_q3 - _df_q1
 
-			_se = concat([_df_q1 - 1.5 * _df_iqr] * len(_df_), axis=1).T.set_index(_df_.index)
-			_le = concat([_df_q3 + 1.5 * _df_iqr] * len(_df_), axis=1).T.set_index(_df_.index)
+            _se = concat([_df_q1 - 1.5 * _df_iqr] * len(_df_), axis=1).T.set_index(_df_.index)
+            _le = concat([_df_q3 + 1.5 * _df_iqr] * len(_df_), axis=1).T.set_index(_df_.index)
 
-			if _log:
-				return 10 ** _se, 10 ** _le
-			return _se, _le
+            if _log:
+                return 10 ** _se, 10 ** _le
+            return _se, _le
 
-		_cation, _anion, _main = ['Na+', 'NH4+', 'K+', 'Mg2+', 'Ca2+'], ['Cl-', 'NO2-', 'NO3-', 'SO42-', ], ['SO42-',
-																											 'NO3-',
-																											 'NH4+']
+        _cation, _anion, _main = ['Na+', 'NH4+', 'K+', 'Mg2+', 'Ca2+'], ['Cl-', 'NO2-', 'NO3-', 'SO42-', ], ['SO42-',
+                                                                                                             'NO3-',
+                                                                                                             'NH4+']
 
-		_df_salt = _df[_mdl.keys()].copy()
-		_df_pm = _df['PM2.5'].copy()
+        _df_salt = _df[_mdl.keys()].copy()
+        _df_pm = _df['PM2.5'].copy()
 
-		## lower than PM2.5
-		## conc. of main salt should be present at the same time (NH4+, SO42-, NO3-)
-		_df_salt = _df_salt.mask(_df_salt.sum(axis=1, min_count=1) > _df_pm).dropna(subset=_main).copy()
+        ## lower than PM2.5
+        ## conc. of main salt should be present at the same time (NH4+, SO42-, NO3-)
+        _df_salt = _df_salt.mask(_df_salt.sum(axis=1, min_count=1) > _df_pm).dropna(subset=_main).copy()
 
-		## mdl
-		for (_key, _df_col), _mdl_val in zip(_df_salt.items(), _mdl.values()):
-			_df_salt[_key] = _df_col.mask(_df_col < _mdl_val, _mdl_val / 2)
+        ## mdl
+        for (_key, _df_col), _mdl_val in zip(_df_salt.items(), _mdl.values()):
+            _df_salt[_key] = _df_col.mask(_df_col < _mdl_val, _mdl_val / 2)
 
-		## group by time (per month)
-		_df_salt['tm'] = _df_salt.index.strftime('%Y-%m')
+        ## group by time (per month)
+        _df_salt['tm'] = _df_salt.index.strftime('%Y-%m')
 
-		_df_lst = []
-		for _ky, _df_grp in _df_salt.groupby('tm'):
-			_df_grp = _df_grp[_mdl.keys()].copy()
+        _df_lst = []
+        for _ky, _df_grp in _df_salt.groupby('tm'):
+            _df_grp = _df_grp[_mdl.keys()].copy()
 
-			## calculate SE LE
-			## salt < LE
-			_se, _le = _se_le(_df_grp, _log=True)
-			_df_grp = _df_grp.mask(_df_grp > _le).copy()
+            ## calculate SE LE
+            ## salt < LE
+            _se, _le = _se_le(_df_grp, _log=True)
+            _df_grp = _df_grp.mask(_df_grp > _le).copy()
 
-			## C/A, A/C
-			_rat_CA = (_df_grp[_cation].sum(axis=1) / _df_grp[_anion].sum(axis=1)).to_frame()
-			_rat_AC = (1 / _rat_CA).copy()
+            ## C/A, A/C
+            _rat_CA = (_df_grp[_cation].sum(axis=1) / _df_grp[_anion].sum(axis=1)).to_frame()
+            _rat_AC = (1 / _rat_CA).copy()
 
-			_se, _le = _se_le(_rat_CA, )
-			_cond_CA = (_rat_CA < _le) & (_rat_CA > 0)
+            _se, _le = _se_le(_rat_CA, )
+            _cond_CA = (_rat_CA < _le) & (_rat_CA > 0)
 
-			_se, _le = _se_le(_rat_AC, )
-			_cond_AC = (_rat_AC < _le) & (_rat_AC > 0)
+            _se, _le = _se_le(_rat_AC, )
+            _cond_AC = (_rat_AC < _le) & (_rat_AC > 0)
 
-			_df_grp = _df_grp.where((_cond_CA * _cond_AC)[0]).copy()
+            _df_grp = _df_grp.where((_cond_CA * _cond_AC)[0]).copy()
 
-			## conc. of main salt > SE
-			_se, _le = _se_le(_df_grp[_main], _log=True)
-			_df_grp[_main] = _df_grp[_main].mask(_df_grp[_main] < _se).copy()
+            ## conc. of main salt > SE
+            _se, _le = _se_le(_df_grp[_main], _log=True)
+            _df_grp[_main] = _df_grp[_main].mask(_df_grp[_main] < _se).copy()
 
-			_df_lst.append(_df_grp)
+            _df_lst.append(_df_grp)
 
-		return concat(_df_lst).reindex(_df.index)
+        return concat(_df_lst).reindex(_df.index)