profilegenerator.py

#!/usr/bin/python3

    #Copyright (C) 2023 University of Twente

    #This program is free software: you can redistribute it and/or modify
    #it under the terms of the GNU General Public License as published by
    #the Free Software Foundation, either version 3 of the License, or
    #(at your option) any later version.

    #This program is distributed in the hope that it will be useful,
    #but WITHOUT ANY WARRANTY; without even the implied warranty of
    #MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    #GNU General Public License for more details.

    #You should have received a copy of the GNU General Public License
    #along with this program.  If not, see <http://www.gnu.org/licenses/>.

    
#from configLoader import *
#config = importlib.import_module(cfgFile)
import os, sys, getopt

print("Profilegenerator 1.3.2\n", flush=True)
print("Copyright (C) 2023 University of Twente", flush=True)
print("This program comes with ABSOLUTELY NO WARRANTY.", flush=True)
print("This is free software, and you are welcome to redistribute it under certain conditions.", flush=True)
print("See the acompanying license for more information.\n", flush=True)

if(len(sys.argv) > 1):
	#Get arguments:
	try:
		opts, args = getopt.getopt(sys.argv[1:],"c:o:f",["config=","output=", "force"])
	except getopt.GetoptError:
		print("Usage:", flush=True)
		print('profilegenerator.py -c <config> [-o <output subfolder> -f]', flush=True)
		sys.exit(2)

	# Default write into a new folder
	cfgOutputDir = 'output/output/'
	forceDeletion = False

	#Parse arguments
	for opt, arg in opts:
		if opt in ("-c", "--config"):
			cfgFile = arg
		elif opt in ("-o", "--output"):
			cfgOutputDir = 'output/'+arg+'/'
		elif opt in ("-f", "--force"):
			forceDeletion = True



	#Check if the output dir exists, otherwise make it
	os.makedirs(os.path.dirname(cfgOutputDir), exist_ok=True)

	if os.listdir(cfgOutputDir):
		# Empty the directory
		if forceDeletion:
			for tf in os.listdir(cfgOutputDir):
				fp = os.path.join(cfgOutputDir, tf)
				try:
					if os.path.isfile(fp):
						os.unlink(fp)
				except Exception as e:
					print(e, flush=True)
		else:
			print("Config directory is not empty! Provide the --force flag to delete the contents", flush=True)
			exit()
			
			


else:
	print("Usage:", flush=True)
	print('profilegenerator.py -c <config> [-o <output subfolder> -f]', flush=True)
	exit()



import profilegentools
from configLoader import *
config = importlib.import_module(cfgFile)

print('Loading config: '+cfgFile, flush=True)
print("The current config will create and simulate "+str(len(config.householdList))+" households", flush=True)
print("Results will be written into: "+cfgOutputDir+"\n", flush=True)
print("NOTE: Simulation may take a (long) while...\n", flush=True)

# Check the config:
if config.penetrationEV + config.penetrationPHEV > 100:
	print("Error, the combined penetration of EV and PHEV exceed 100!", flush=True)
	exit()
if config.penetrationPV < config.penetrationBattery:
	print("Error, the penetration of PV must be equal or higher than PV!", flush=True)
	exit()
if config.penetrationHeatPump + config.penetrationCHP > 100:
	print("Error, the combined penetration of heatpumps and CHPs exceed 100!", flush=True)
	exit()



# Randomize using the seed
random.seed(config.seed)

# Create empty files
config.writer.createEmptyFiles()

#import neighbourhood
import neighbourhood
neighbourhood.neighbourhood()

configFile = []

hnum = 0

householdList = config.householdList
numOfHouseholds = len(householdList)

while len(householdList) > 0:
	print("Household "+str(hnum+1)+" of "+str(numOfHouseholds), flush=True)
	householdList[0].simulate()
	
	#Warning: On my PC the random number is still the same at this point, but after calling scaleProfile() it isn't!!!
	householdList[0].scaleProfile()
	householdList[0].reactivePowerProfile()
	householdList[0].thermalGainProfile()

	config.writer.writeHousehold(householdList[0], hnum)
	config.writer.writeNeighbourhood(hnum)
	
	householdList[0].Consumption = None
	householdList[0].Occupancy = None
	for p in householdList[0].Persons:
		del(p)
	del(householdList[0])
	
	hnum = hnum + 1