check_communes.py
View Options

			import xarray as xr
			import sys
			import os

			fileDir = sys.argv[1]
			batches = int(sys.argv[2])

			year = int(fileDir.split('_')[-1])

			tmpDir = 'tmp_%d' %year

			for batch in range(batches):
			print(batch)
			comm_fp = os.path.join(fileDir, '%d_commune_hourly' %year, '%d_commune_hourly_%d.nc' %(year, batch))
			irrad_fp = os.path.join(tmpDir, 'tilted_irrad_%d.nc' %batch)
			pvPot_fp = os.path.join(tmpDir, 'pv_potential_%d.nc' %batch)

			try:
			commPot = xr.open_dataset(comm_fp)

			if commPot.pv_potential.sum() > 0:
			continue

			print('ERROR in %s' %pvPot_fp)
			commPot.close()

			except:
			print('Cannot open %s' %comm_fp)

			os.system('python merge_commune_batch.py %s %s %s' %(pvPot_fp, irrad_fp, comm_fp))

solarPOT/Potential_hourly/check_communes.py29f993e4425aaw_solar