clean_tilted_irrad.py
b9297dda7962aw_solar

clean_tilted_irrad.py
View Options

			import os
			import pandas as pd
			import xarray as xr
			import numpy as np

			import time
			import util

			IRRAD_FP = sys.argv[1]
			OUT_FP = sys.argv[2]

			data = xr.open_mfdataset( IRRAD_FP )

			data['poa_global'] = data.poa_direct + data.poa_sky_diffuse + data.poa_ground_diffuse
			outliers = data.where(data.poa_global > 2000, drop = True)
			good_timestamps = data.timestamp[ np.logical_not( data.timestamp.isin(outliers.timestamp) )]

			poa_irrad_clean = data.sel( timestamp = good_timestamps).drop(['poa_global'] )
			poa_irrad_clean.to_netcdf( OUT_FP )