for i in range(len(vRange)):
for j in range(len(alphaRange)):
valueList.append({'velocity':vRange[i],'alpha':alphaRange[j]})
valueList.append({'velocity':10,'alpha':4})
valueList.append({'velocity':5,'alpha':10})
valueList.append({'velocity':7.5,'alpha':10})
print valueList
for k in range(6,len(valueList)):
vel = valueList[k]['velocity']
alpha = valueList[k]['alpha']
caseDir.clear(functionObjectData=True)
caseDir.clearOther()
caseDir.clearPattern("postProcessing")
caseDir.clearResults(functionObjectData=True)
velocity = [0,math.sin(alpha*math.pi/180)*vel,-math.cos(alpha*math.pi/180)*vel]
print "Simulation {:d} from {:d}: velocity = {:f} alpha = {:f}".format(k,len(valueList),vel,alpha)
parameters['flowVelocity'] = "({:f} {:f} {:f})".format(velocity[0],velocity[1],velocity[2])
parameters['magUInf'] = "{:f}".format(vel)
parameters.writeFile()
decompose = BasicRunner(argv=["decomposePar"],silent=True)
print " Decompose Dictionary"
start = time.clock()
decompose.start()
print " Run SimpleFoam"
run = AnalyzedRunner(StandardLogAnalyzer(),argv=["mpirun -np 4","simpleFoam","-parallel"],silent=True)
run.start()
run.picklePlots()
print " Reconstruct Par"
