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

AeroViz/dataProcess/Chemistry/_mass_volume.py

@@ -1,175 +1,175 @@
-from pandas import date_range, concat, DataFrame, to_numeric
+from pandas import concat, DataFrame
 
 
 def _basic(df_che, df_ref, df_water, df_density, nam_lst):
-	df_all = concat(df_che, axis=1)
-	index = df_all.index.copy()
-	df_all.columns = nam_lst
-
-	## parameter
-	mol_A, mol_S, mol_N = df_all['NH4+'] / 18, df_all['SO42-'] / 96, df_all['NO3-'] / 62
-	df_all['status'] = (mol_A) / (2 * mol_S + mol_N)
-
-	convert_nam = {'AS': 'SO42-',
-				   'AN': 'NO3-',
-				   'OM': 'OC',
-				   'Soil': 'Fe',
-				   'SS': 'Na+',
-				   'EC': 'EC',
-				   }
-
-	mass_coe = {'AS': 1.375,
-				'AN': 1.29,
-				'OM': 1.8,
-				'Soil': 28.57,
-				'SS': 2.54,
-				'EC': 1,
-				}
-
-	vol_coe = {'AS': 1.76,
-			   'AN': 1.73,
-			   'OM': 1.4,
-			   'Soil': 2.6,
-			   'SS': 2.16,
-			   'EC': 1.5,
-			   }
-
-	RI_coe = {'550': {'ALWC': 1.333 + 0j,
-					  'AS': 1.53 + 0j,
-					  'AN': 1.55 + 0j,
-					  'OM': 1.55 + 0.0163j,
-					  'Soil': 1.56 + 0.006j,
-					  'SS': 1.54 + 0j,
-					  'EC': 1.80 + 0.72j,
-					  },
-
-			  ## m + kj -> m value is same as 550 current
-			  '450': {'ALWC': 1.333 + 0j,
-					  'AS': 1.57 + 0j,
-					  'AN': 1.57 + 0j,
-					  'OM': 1.58 + 0.056,
-					  'Soil': 1.56 + 0.009j,
-					  'SS': 1.54 + 0j,
-					  'EC': 1.80 + 0.79j,
-					  },
-			  }
-
-	## mass
-	## NH4 Enough
-	df_mass = DataFrame()
-	df_enough = df_all.where(df_all['status'] >= 1).dropna().copy()
-
-	for _mass_nam, _coe in mass_coe.items():
-		df_mass[_mass_nam] = df_all[convert_nam[_mass_nam]] * _coe
-
-	## NH4 Deficiency
-	defic_idx = df_all['status'] < 1
-
-	if defic_idx.any():
-		residual = mol_A - 2 * mol_S
-
-		## residual > 0
-		_status = residual > 0
-		if _status.any():
-			_cond = _status & (residual <= mol_N)
-			df_mass.loc[_cond, 'AN'] = residual.loc[_cond] * 80
-
-			_cond = _status & (residual > mol_N)
-			df_mass.loc[_cond, 'AN'] = mol_N.loc[_cond] * 80
-
-		## residual < 0
-		_status = residual <= 0
-		if _status.any():
-			df_mass.loc[_status, 'AN'] = 0
-
-			_cond = _status & (mol_A <= 2 * mol_S)
-			df_mass.loc[_cond, 'AS'] = mol_A.loc[_cond] / 2 * 132
-
-			_cond = _status & (mol_A > 2 * mol_S)
-			df_mass.loc[_cond, 'AS'] = mol_S.loc[_cond] * 132
-
-	df_mass_cal = df_mass.dropna().copy()
-	df_mass['total'] = df_mass.sum(axis=1, min_count=6)
-
-	qc_ratio = df_mass['total'] / df_ref
-	qc_cond = (qc_ratio >= 0.7) & (qc_ratio <= 1.3)
-
-	## volume
-	df_vol = DataFrame()
-	for _vol_nam, _coe in vol_coe.items():
-		df_vol[_vol_nam] = df_mass_cal[_vol_nam] / _coe
-
-	if df_water is not None:
-		df_vol['ALWC'] = df_water
-		df_vol = df_vol.dropna()
-		df_vol['total_wet'] = df_vol.sum(axis=1, min_count=6)
-
-	df_vol['total_dry'] = df_vol[vol_coe.keys()].sum(axis=1, min_count=6)
-
-	## density
-	df_vol_cal = DataFrame()
-	df_den_rec = df_mass['total'] / df_vol['total_dry']
-	if df_density is not None:
-		df_den_all = concat([df_all[['SO42-', 'NO3-', 'NH4+', 'EC']], df_density, df_mass['OM']], axis=1).dropna()
-
-		df_vol_cal = (df_den_all[['SO42-', 'NO3-', 'NH4+']].sum(axis=1) / 1.75) + \
-					 df_den_all['Cl-'] / 1.52 + \
-					 df_den_all['OM'] / 1.4 + df_den_all['EC'] / 1.77
-
-		df_den = df_den_all.sum(axis=1, min_count=6) / df_vol_cal
-	# df_den = df_den_all.sum(axis=1) / df_vol_cal
-	# df_den = df_mass['total'].loc[df_den_all.index] / df_vol_cal
-
-	## refractive index
-	ri_dic = {}
-	for _lambda, _coe in RI_coe.items():
-
-		df_RI = DataFrame()
-
-		for _ky, _df in df_vol.items():
-			if 'total' in _ky: continue
-			df_RI[_ky] = (_df * _coe[_ky])
-
-		df_RI['RI_wet'] = None
-		if df_water is not None:
-			df_RI['RI_wet'] = (df_RI / df_vol['total_wet'].to_frame().values).sum(axis=1)
+    df_all = concat(df_che, axis=1)
+    index = df_all.index.copy()
+    df_all.columns = nam_lst
+
+    ## parameter
+    mol_A, mol_S, mol_N = df_all['NH4+'] / 18, df_all['SO42-'] / 96, df_all['NO3-'] / 62
+    df_all['status'] = (mol_A) / (2 * mol_S + mol_N)
+
+    convert_nam = {'AS': 'SO42-',
+                   'AN': 'NO3-',
+                   'OM': 'OC',
+                   'Soil': 'Fe',
+                   'SS': 'Na+',
+                   'EC': 'EC',
+                   }
+
+    mass_coe = {'AS': 1.375,
+                'AN': 1.29,
+                'OM': 1.8,
+                'Soil': 28.57,
+                'SS': 2.54,
+                'EC': 1,
+                }
+
+    vol_coe = {'AS': 1.76,
+               'AN': 1.73,
+               'OM': 1.4,
+               'Soil': 2.6,
+               'SS': 2.16,
+               'EC': 1.5,
+               }
+
+    RI_coe = {'550': {'ALWC': 1.333 + 0j,
+                      'AS': 1.53 + 0j,
+                      'AN': 1.55 + 0j,
+                      'OM': 1.55 + 0.0163j,
+                      'Soil': 1.56 + 0.006j,
+                      'SS': 1.54 + 0j,
+                      'EC': 1.80 + 0.72j,
+                      },
+
+              ## m + kj -> m value is same as 550 current
+              '450': {'ALWC': 1.333 + 0j,
+                      'AS': 1.57 + 0j,
+                      'AN': 1.57 + 0j,
+                      'OM': 1.58 + 0.056,
+                      'Soil': 1.56 + 0.009j,
+                      'SS': 1.54 + 0j,
+                      'EC': 1.80 + 0.79j,
+                      },
+              }
+
+    ## mass
+    ## NH4 Enough
+    df_mass = DataFrame()
+    df_enough = df_all.where(df_all['status'] >= 1).dropna().copy()
+
+    for _mass_nam, _coe in mass_coe.items():
+        df_mass[_mass_nam] = df_all[convert_nam[_mass_nam]] * _coe
+
+    ## NH4 Deficiency
+    defic_idx = df_all['status'] < 1
+
+    if defic_idx.any():
+        residual = mol_A - 2 * mol_S
+
+        ## residual > 0
+        _status = residual > 0
+        if _status.any():
+            _cond = _status & (residual <= mol_N)
+            df_mass.loc[_cond, 'AN'] = residual.loc[_cond] * 80
+
+            _cond = _status & (residual > mol_N)
+            df_mass.loc[_cond, 'AN'] = mol_N.loc[_cond] * 80
+
+        ## residual < 0
+        _status = residual <= 0
+        if _status.any():
+            df_mass.loc[_status, 'AN'] = 0
+
+            _cond = _status & (mol_A <= 2 * mol_S)
+            df_mass.loc[_cond, 'AS'] = mol_A.loc[_cond] / 2 * 132
+
+            _cond = _status & (mol_A > 2 * mol_S)
+            df_mass.loc[_cond, 'AS'] = mol_S.loc[_cond] * 132
+
+    df_mass_cal = df_mass.dropna().copy()
+    df_mass['total'] = df_mass.sum(axis=1, min_count=6)
+
+    qc_ratio = df_mass['total'] / df_ref
+    qc_cond = (qc_ratio >= 0.7) & (qc_ratio <= 1.3)
+
+    ## volume
+    df_vol = DataFrame()
+    for _vol_nam, _coe in vol_coe.items():
+        df_vol[_vol_nam] = df_mass_cal[_vol_nam] / _coe
+
+    if df_water is not None:
+        df_vol['ALWC'] = df_water
+        df_vol = df_vol.dropna()
+        df_vol['total_wet'] = df_vol.sum(axis=1, min_count=6)
+
+    df_vol['total_dry'] = df_vol[vol_coe.keys()].sum(axis=1, min_count=6)
+
+    ## density
+    df_vol_cal = DataFrame()
+    df_den_rec = df_mass['total'] / df_vol['total_dry']
+    if df_density is not None:
+        df_den_all = concat([df_all[['SO42-', 'NO3-', 'NH4+', 'EC']], df_density, df_mass['OM']], axis=1).dropna()
+
+        df_vol_cal = (df_den_all[['SO42-', 'NO3-', 'NH4+']].sum(axis=1) / 1.75) + \
+                     df_den_all['Cl-'] / 1.52 + \
+                     df_den_all['OM'] / 1.4 + df_den_all['EC'] / 1.77
+
+        df_den = df_den_all.sum(axis=1, min_count=6) / df_vol_cal
+    # df_den = df_den_all.sum(axis=1) / df_vol_cal
+    # df_den = df_mass['total'].loc[df_den_all.index] / df_vol_cal
+
+    ## refractive index
+    ri_dic = {}
+    for _lambda, _coe in RI_coe.items():
+
+        df_RI = DataFrame()
+
+        for _ky, _df in df_vol.items():
+            if 'total' in _ky: continue
+            df_RI[_ky] = (_df * _coe[_ky])
+
+        df_RI['RI_wet'] = None
+        if df_water is not None:
+            df_RI['RI_wet'] = (df_RI / df_vol['total_wet'].to_frame().values).sum(axis=1)
 
-		df_RI['RI_dry'] = (df_RI[vol_coe.keys()] / df_vol['total_dry'].to_frame().values).sum(axis=1)
+        df_RI['RI_dry'] = (df_RI[vol_coe.keys()] / df_vol['total_dry'].to_frame().values).sum(axis=1)
 
-		ri_dic[f'RI_{_lambda}'] = df_RI[['RI_dry', 'RI_wet']]
+        ri_dic[f'RI_{_lambda}'] = df_RI[['RI_dry', 'RI_wet']]
 
-	## mole and equivalent
-	df_eq = concat((mol_A, mol_S, mol_N, mol_A * 1, mol_S * 2, mol_N * 1), axis=1)
-	df_eq.columns = ['mol_NH4', 'mol_SO4', 'mol_NO3', 'eq_NH4', 'eq_SO4', 'eq_NO3', ]
+    ## mole and equivalent
+    df_eq = concat((mol_A, mol_S, mol_N, mol_A * 1, mol_S * 2, mol_N * 1), axis=1)
+    df_eq.columns = ['mol_NH4', 'mol_SO4', 'mol_NO3', 'eq_NH4', 'eq_SO4', 'eq_NO3', ]
 
-	## out
-	out = {'mass': df_mass,
-		   'volume': df_vol,
-		   'vol_cal': df_vol_cal,
-		   'eq': df_eq,
-		   'density_mat': df_den,
-		   'density_rec': df_den_rec,
-		   }
-	out.update(ri_dic)
+    ## out
+    out = {'mass': df_mass,
+           'volume': df_vol,
+           'vol_cal': df_vol_cal,
+           'eq': df_eq,
+           'density_mat': df_den,
+           'density_rec': df_den_rec,
+           }
+    out.update(ri_dic)
 
-	for _ky, _df in out.items():
-		out[_ky] = _df.reindex(index).where(qc_cond)
+    for _ky, _df in out.items():
+        out[_ky] = _df.reindex(index).where(qc_cond)
 
-	return out
+    return out
 
 
 # '''
 
 
 def mass_ratio(_df):
-	if _df['PM25'] >= _df['total_mass']:
-		_df['others'] = _df['PM25'] - _df['total_mass']
-		for _val, _species in zip(_df.values, _df.index):
-			_df[f'{_species}_ratio'] = _val / _df['PM25'].__round__(3)
+    if _df['PM25'] >= _df['total_mass']:
+        _df['others'] = _df['PM25'] - _df['total_mass']
+        for _val, _species in zip(_df.values, _df.index):
+            _df[f'{_species}_ratio'] = _val / _df['PM25'].__round__(3)
 
-	if _df['PM25'] < _df['total_mass']:
-		_df['others'] = 0
-		for _val, _species in zip(_df.values, _df.index):
-			_df[f'{_species}_ratio'] = _val / _df['PM25'].__round__(3)
+    if _df['PM25'] < _df['total_mass']:
+        _df['others'] = 0
+        for _val, _species in zip(_df.values, _df.index):
+            _df[f'{_species}_ratio'] = _val / _df['PM25'].__round__(3)
 
-	return _df['others':].drop(labels=['PM25_ratio', 'total_mass_ratio'])
+    return _df['others':].drop(labels=['PM25_ratio', 'total_mass_ratio'])

AeroViz/dataProcess/Chemistry/_ocec.py

@@ -1,161 +1,161 @@
-from AeroViz.dataProcess.core import _union_index
-
-from pandas import date_range, concat, DataFrame, to_numeric
-from scipy.optimize import curve_fit
 import numpy as np
+from pandas import concat, DataFrame
+from scipy.optimize import curve_fit
+
+from AeroViz.dataProcess.core import _union_index
 
 __all__ = [
-	'_basic',
-	# '_ocec_ratio_cal',
+    '_basic',
+    # '_ocec_ratio_cal',
 ]
 
 
 def _min_Rsq(_oc, _ec, _rng):
-	_val_mesh, _oc_mesh = np.meshgrid(_rng, _oc)
-	_val_mesh, _ec_mesh = np.meshgrid(_rng, _ec)
+    _val_mesh, _oc_mesh = np.meshgrid(_rng, _oc)
+    _val_mesh, _ec_mesh = np.meshgrid(_rng, _ec)
 
-	_out_table = DataFrame(_oc_mesh - _val_mesh * _ec_mesh, index=_oc.index, columns=_rng)
+    _out_table = DataFrame(_oc_mesh - _val_mesh * _ec_mesh, index=_oc.index, columns=_rng)
 
-	## calculate R2
-	_r2_dic = {}
-	_func = lambda _x, _sl, _inte: _sl * _x + _inte
-	for _ocec, _out in _out_table.items():
-		_df = DataFrame([_out.values, _ec.values]).T.dropna()
+    ## calculate R2
+    _r2_dic = {}
+    _func = lambda _x, _sl, _inte: _sl * _x + _inte
+    for _ocec, _out in _out_table.items():
+        _df = DataFrame([_out.values, _ec.values]).T.dropna()
 
-		_x, _y = _df[0], _df[1]
-		_opt, _ = curve_fit(_func, _x, _y)
+        _x, _y = _df[0], _df[1]
+        _opt, _ = curve_fit(_func, _x, _y)
 
-		_tss = np.sum((_y - _y.mean()) ** 2.)
-		_rss = np.sum((_y - _func(_x, *_opt)) ** 2.)
+        _tss = np.sum((_y - _y.mean()) ** 2.)
+        _rss = np.sum((_y - _func(_x, *_opt)) ** 2.)
 
-		_r2_dic[round(_ocec, 3)] = 1. - _rss / _tss
+        _r2_dic[round(_ocec, 3)] = 1. - _rss / _tss
 
-	## get the min R2
-	_ratio = DataFrame(_r2_dic, index=[0]).idxmin(axis=1).values[0]
+    ## get the min R2
+    _ratio = DataFrame(_r2_dic, index=[0]).idxmin(axis=1).values[0]
 
-	return _ratio, _out_table[_ratio]
+    return _ratio, _out_table[_ratio]
 
 
 def _ocec_ratio_cal(_nam, _lcres_splt, _hr_lim, _range_, _wisoc_range_):
-	## parameter
-	_out = DataFrame(index=_lcres_splt.index)
-	(_, _oc), (_, _ec) = _lcres_splt.items()
-	# _oc, _ec = _lcres_splt['Thermal_OC'], _lcres_splt['Thermal_EC']
+    ## parameter
+    _out = DataFrame(index=_lcres_splt.index)
+    (_, _oc), (_, _ec) = _lcres_splt.items()
+    # _oc, _ec = _lcres_splt['Thermal_OC'], _lcres_splt['Thermal_EC']
 
-	## real data OC/EC
-	_ocec_ratio_real = (_oc / _ec).quantile(.5)
+    ## real data OC/EC
+    _ocec_ratio_real = (_oc / _ec).quantile(.5)
 
-	_out[f'OC/EC_real_{_nam}'] = _ocec_ratio_real
-	_out[f'POC_real_{_nam}'] = _ocec_ratio_real * _ec
-	_out[f'SOC_real_{_nam}'] = _oc - _out[f'POC_real_{_nam}']
+    _out[f'OC/EC_real_{_nam}'] = _ocec_ratio_real
+    _out[f'POC_real_{_nam}'] = _ocec_ratio_real * _ec
+    _out[f'SOC_real_{_nam}'] = _oc - _out[f'POC_real_{_nam}']
 
-	## the least R2 method
-	## estimated OC/EC
-	if (len(_lcres_splt) <= _hr_lim):
-		print(f"\t\t{_lcres_splt.index[0].strftime('%Y-%m-%d %X')} to {_lcres_splt.index[-1].strftime('%Y-%m-%d %X')}")
-		print('\t\tPlease Modify the Values of "hour_limit" or Input Sufficient Amount of Data !!')
+    ## the least R2 method
+    ## estimated OC/EC
+    if (len(_lcres_splt) <= _hr_lim):
+        print(f"\t\t{_lcres_splt.index[0].strftime('%Y-%m-%d %X')} to {_lcres_splt.index[-1].strftime('%Y-%m-%d %X')}")
+        print('\t\tPlease Modify the Values of "hour_limit" or Input Sufficient Amount of Data !!')
 
-		_out[[f'OC/EC_{_nam}', f'POC_{_nam}', f'SOC_{_nam}', f'WISOC/OC_{_nam}', f'WSOC_{_nam}',
-			  f'WISOC_{_nam}']] = np.nan
+        _out[[f'OC/EC_{_nam}', f'POC_{_nam}', f'SOC_{_nam}', f'WISOC/OC_{_nam}', f'WSOC_{_nam}',
+              f'WISOC_{_nam}']] = np.nan
 
-		return _out
+        return _out
 
-	if (len(_lcres_splt.dropna()) == 0):
-		_out[[f'OC/EC_{_nam}', f'POC_{_nam}', f'SOC_{_nam}', f'WISOC/OC_{_nam}', f'WSOC_{_nam}',
-			  f'WISOC_{_nam}']] = np.nan
+    if (len(_lcres_splt.dropna()) == 0):
+        _out[[f'OC/EC_{_nam}', f'POC_{_nam}', f'SOC_{_nam}', f'WISOC/OC_{_nam}', f'WSOC_{_nam}',
+              f'WISOC_{_nam}']] = np.nan
 
-		return _out
+        return _out
 
-	## OC/EC
-	_ocec_ratio = False
-	_st, _ed, _stp = _range_
+    ## OC/EC
+    _ocec_ratio = False
+    _st, _ed, _stp = _range_
 
-	for _ in range(2):
-		if _ocec_ratio:
-			_ocec_rng = np.arange(_ocec_ratio - _stp / 2, _ocec_ratio + _stp / 2, .01).round(3)
-		else:
-			_ocec_rng = np.arange(_st, _ed + _stp, _stp).round(3)
+    for _ in range(2):
+        if _ocec_ratio:
+            _ocec_rng = np.arange(_ocec_ratio - _stp / 2, _ocec_ratio + _stp / 2, .01).round(3)
+        else:
+            _ocec_rng = np.arange(_st, _ed + _stp, _stp).round(3)
 
-		_ocec_ratio, _soc = _min_Rsq(_oc, _ec, _ocec_rng)
+        _ocec_ratio, _soc = _min_Rsq(_oc, _ec, _ocec_rng)
 
-	## WISOC
-	_st, _ed, _stp = _wisoc_range_
-	_wisoc_rng = (np.arange(_st, _ed + _stp, _stp) * _ocec_ratio).round(5)
-	_wisoc_ratio, _wsoc = _min_Rsq(_oc, _ec, _wisoc_rng)
+    ## WISOC
+    _st, _ed, _stp = _wisoc_range_
+    _wisoc_rng = (np.arange(_st, _ed + _stp, _stp) * _ocec_ratio).round(5)
+    _wisoc_ratio, _wsoc = _min_Rsq(_oc, _ec, _wisoc_rng)
 
-	## out
-	_out[f'OC/EC_{_nam}'] = _ocec_ratio
-	_out[f'SOC_{_nam}'] = _soc
-	_out[f'POC_{_nam}'] = _oc - _out[f'SOC_{_nam}']
-	_out[f'WISOC/OC_{_nam}'] = _wisoc_ratio
-	_out[f'WSOC_{_nam}'] = _wsoc
-	_out[f'WISOC_{_nam}'] = _oc - _out[f'WSOC_{_nam}']
+    ## out
+    _out[f'OC/EC_{_nam}'] = _ocec_ratio
+    _out[f'SOC_{_nam}'] = _soc
+    _out[f'POC_{_nam}'] = _oc - _out[f'SOC_{_nam}']
+    _out[f'WISOC/OC_{_nam}'] = _wisoc_ratio
+    _out[f'WSOC_{_nam}'] = _wsoc
+    _out[f'WISOC_{_nam}'] = _oc - _out[f'WSOC_{_nam}']
 
-	return _out[[f'OC/EC_{_nam}', f'POC_{_nam}', f'SOC_{_nam}', f'WISOC/OC_{_nam}', f'WSOC_{_nam}', f'WISOC_{_nam}',
-				 f'OC/EC_real_{_nam}', f'POC_real_{_nam}', f'SOC_real_{_nam}']]
+    return _out[[f'OC/EC_{_nam}', f'POC_{_nam}', f'SOC_{_nam}', f'WISOC/OC_{_nam}', f'WSOC_{_nam}', f'WISOC_{_nam}',
+                 f'OC/EC_real_{_nam}', f'POC_real_{_nam}', f'SOC_real_{_nam}']]
 
 
 def _basic(_lcres, _res, _mass, _ocec_ratio, _ocec_ratio_month, _hr_lim, _range, _wisoc_range):
-	_lcres, _res, _mass = _union_index(_lcres, _res, _mass)
+    _lcres, _res, _mass = _union_index(_lcres, _res, _mass)
 
-	_out = {}
+    _out = {}
 
-	## OC1, OC2, OC3, OC4, PC
-	_df_bsc = _res / _lcres['Sample_Volume'].to_frame().values.copy()
+    ## OC1, OC2, OC3, OC4, PC
+    _df_bsc = _res / _lcres['Sample_Volume'].to_frame().values.copy()
 
-	## SOC, POC, OC/EC
-	if _ocec_ratio is not None:
-		try:
-			iter(_ocec_ratio)
-		except TypeError:
-			raise TypeError('"ocec_ratio" Only Accept a Single Value !!')
+    ## SOC, POC, OC/EC
+    if _ocec_ratio is not None:
+        try:
+            iter(_ocec_ratio)
+        except TypeError:
+            raise TypeError('"ocec_ratio" Only Accept a Single Value !!')
 
-		_prcs_df = DataFrame(index=_df_bsc.index)
-		_prcs_df['OC/EC'] = _ocec_ratio
-		_prcs_df['POC'] = _ocec_ratio * _lcres['Thermal_EC']
-		_prcs_df['SOC'] = _lcres['Thermal_OC'] - _prcs_df['POC']
+        _prcs_df = DataFrame(index=_df_bsc.index)
+        _prcs_df['OC/EC'] = _ocec_ratio
+        _prcs_df['POC'] = _ocec_ratio * _lcres['Thermal_EC']
+        _prcs_df['SOC'] = _lcres['Thermal_OC'] - _prcs_df['POC']
 
-	else:
-		_df_lst = []
-		for _, _df in _lcres.resample(f'{_ocec_ratio_month}MS', closed='left'):
-			_thm_cal = _ocec_ratio_cal('thm', _df[['Thermal_OC', 'Thermal_EC']], _hr_lim, _range, _wisoc_range)
-			_opt_cal = _ocec_ratio_cal('opt', _df[['Optical_OC', 'Optical_EC']], _hr_lim, _range, _wisoc_range)
-			_df_lst.append(concat([_thm_cal, _opt_cal], axis=1))
+    else:
+        _df_lst = []
+        for _, _df in _lcres.resample(f'{_ocec_ratio_month}MS', closed='left'):
+            _thm_cal = _ocec_ratio_cal('thm', _df[['Thermal_OC', 'Thermal_EC']], _hr_lim, _range, _wisoc_range)
+            _opt_cal = _ocec_ratio_cal('opt', _df[['Optical_OC', 'Optical_EC']], _hr_lim, _range, _wisoc_range)
+            _df_lst.append(concat([_thm_cal, _opt_cal], axis=1))
 
-		_prcs_df = concat(_df_lst)
+        _prcs_df = concat(_df_lst)
 
-	_df_bsc = concat((_df_bsc.copy(), _prcs_df), axis=1)
+    _df_bsc = concat((_df_bsc.copy(), _prcs_df), axis=1)
 
-	## ratio
-	_df_ratio = DataFrame(index=_df_bsc.index)
+    ## ratio
+    _df_ratio = DataFrame(index=_df_bsc.index)
 
-	for _ky, _val in _df_bsc.items():
-		if 'OC/EC' in _ky: continue
-		_df_ratio[f'{_ky}/Thermal_OC'] = _val / _lcres['Thermal_OC']
-		_df_ratio[f'{_ky}/Optical_OC'] = _val / _lcres['Optical_OC']
+    for _ky, _val in _df_bsc.items():
+        if 'OC/EC' in _ky: continue
+        _df_ratio[f'{_ky}/Thermal_OC'] = _val / _lcres['Thermal_OC']
+        _df_ratio[f'{_ky}/Optical_OC'] = _val / _lcres['Optical_OC']
 
-	if _mass is not None:
-		for _ky, _val in _df_bsc.items():
-			_df_ratio[f'{_ky}/PM'] = _val / _mass
+    if _mass is not None:
+        for _ky, _val in _df_bsc.items():
+            _df_ratio[f'{_ky}/PM'] = _val / _mass
 
-		_df_ratio[f'Thermal_OC/PM'] = _lcres['Thermal_OC'] / _mass
-		_df_ratio[f'Thermal_EC/PM'] = _lcres['Thermal_EC'] / _mass
+        _df_ratio[f'Thermal_OC/PM'] = _lcres['Thermal_OC'] / _mass
+        _df_ratio[f'Thermal_EC/PM'] = _lcres['Thermal_EC'] / _mass
 
-		_df_ratio[f'Optical_OC/PM'] = _lcres['Optical_OC'] / _mass
-		_df_ratio[f'Optical_EC/PM'] = _lcres['Optical_EC'] / _mass
+        _df_ratio[f'Optical_OC/PM'] = _lcres['Optical_OC'] / _mass
+        _df_ratio[f'Optical_EC/PM'] = _lcres['Optical_EC'] / _mass
 
-	## ratio status
-	_df_bsc = concat((_lcres, _df_bsc.copy()), axis=1)
+    ## ratio status
+    _df_bsc = concat((_lcres, _df_bsc.copy()), axis=1)
 
-	for _ky, _df in _df_ratio.items():
-		_df_bsc[f'{_ky}_status'] = 'Normal'
-		_df_bsc[f'{_ky}_status'] = _df_bsc[f'{_ky}_status'].mask(_df > 1, 'Warning')
+    for _ky, _df in _df_ratio.items():
+        _df_bsc[f'{_ky}_status'] = 'Normal'
+        _df_bsc[f'{_ky}_status'] = _df_bsc[f'{_ky}_status'].mask(_df > 1, 'Warning')
 
-	## out
-	_out['ratio'] = _df_ratio
-	_out['basic'] = _df_bsc
+    ## out
+    _out['ratio'] = _df_ratio
+    _out['basic'] = _df_bsc
 
-	return _out
+    return _out
 
 
 '''