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)