sums=[]
resFile=open(resFilePath,"a")
resFile.write(str(id)+";"+str(x)+";"+str(y))
pdb.set_trace()
template.modifyParam("GRID :",str(int(grid_llx))+" "+str(int(grid_lly))+" 1000 1000 1 1")
#The results are stored as a nested list. List edbResults contains one post for each edb. This post contains the sums for all radiuses
sumVelTA=0
edbResults=[]
for edb in edbs:
template.modifyParam("#1015 0 0",cf.user+" "+edb)
template.modifyParam("edb.edb:",edb)
template.modifyParam("EDB :",edb)
report(log,"Executing xrepedb for edb: "+edb)
rastList=execute("xrepedb -i "+cf.template)
radResults=[]
ncellsNonZero=[]
for i in range(len(cf.rads)):
radResults.append(0)
ncellsNonZero.append(0)
report(log,"Processing grid")
for row in rastList[4:-1]:
cols=row.split()
cell_x=float(cols[0])
cell_y=float(cols[1])
cell_val=float(cols[2])
sums=[]
#Summation presumes that the radiuses are in a ordered list with the smallest radius first
for i in range(len(cf.rads)):
sys.exit("Cannot find parameter old_edb in control file")
if parameters==None:
sys.exit("No parameters specified for update in control file")
if kommuner==None:
sys.exit("No kommun specified for update")
report(log,"\n"+50*"=",reportToPrompt)
report(log,"Kommuner to update:",reportToPrompt)
for kn in kommuner:
report(log,kn+", edb: "+knDict[kn],reportToPrompt)
report(log,20*"="+"UPDATING"+20*"=",reportToPrompt)
for kn in kommuner:
report(log,"Processing kommun: "+kn,reportToPrompt)
report(log,"edb: "+knDict[kn],reportToPrompt)
oldFileName=path.join(wd,"old_"+kn+"_"+knDict[kn]+".asc")
old_cmd="roaddb -o"+oldFileName+" -u "+kn+" -e "+knDict[kn]
logs=execute(old_cmd)
report(log,"edb: "+new_edb,reportToPrompt)
newFileName=path.join(wd,"new_"+kn+"_"+new_edb+".asc")
new_cmd="roaddb -o"+newFileName+" -u "+kn+" -e "+new_edb
logs=execute(new_cmd)
outFileName=path.join(od,"updated_"+kn+"_"+new_edb+".asc")
roaddb.updateRoadFile(newFileName,oldFileName,outFileName,parameters)
report(log,"Road ascii db updated",reportToPrompt)
report(log,25*"-",reportToPrompt)
report(log,"Update procedure finished",reportToPrompt)
