import numpy
execfile('helperTwoD.py')
trialsPer = 3 # Number of times to run each experiment

materials = getMaterialList() # Load in material Data
matrix = "ESBR" # Choose the first experiment from TCNanoFillers
fillers = ["Alumina", "ZincOxide"] # Two fillers associated with first experiment

fileName = matrix+fillers[0] + fillers[1] # Name file after matrix, filler materials
f = open(fileName, 'w')
f.write('Matrix\tFiller\tPortion\tRadius\tNumber\tSide\tSeed\tNodes\tq\tdT\tk\tWarn\n')

for i in range(len(materials[matrix]['fillers'])): # "For each filler material"
	modelObject, modelName = createModel(1) # Create model database "Model-1"
	side, radius, portions, dP, dM = defExperiment(matrix, fillers[i]) # Define material attributes for specified matrix, fillers
	
	for j in range(len(portions)): # "For each PHR/volume portion specified"
		for k in range(trialsPer): # "Run trialsPer times for each material/filler/PHR combination"
			seed = numpy.random.randint(1000) # Random seed for coordinate generation
			radius, number = invPHR(portions[j], dP, radius, dM, side) # Returns specific radius size and number of inclusions for closest PHR value.
			xVals, yVals, warningPoints = getPoints(seed, side, radius, number) # returns coordinates for inclusion locations. 
			part = createMatrix(side) # Create the matrix
			edges1, vertices1, face1 = assignGeomSequence(part) # Create references to important sets in our part
			matrixSet, fillerSet, allSet = createCircleInclusion(radius, number, xVals, yVals) # Draw inclusions in the matrix
			createSection(matrix, matrixSet) # Create section for matrix material
			createSection(fillers[i], fillerSet) # Create section for filler material
			assemblyTop, assemblyBottom = makeAssembly() # Create assembly and return references to assembly sets
			temp1, temp2 = 328.15, 298.15 # Assign heat temperature to be used in experiment
			heatStep(temp1, temp2) # apply heat BC
			elements, nodes = makeMesh() # Draw mesh and return number of nodes and elements
			warningString = submitJob() # Submit job and take note of any warnings
			avgHF, TC = getThermalProperties() # Extract relevant information about thermal properties
			f.write(dataString(matrix, fillers[i], portions[j])) # Write the data to file
	
f.close()