import project, os from pyvolution.data_manager import sww2dem from pyvolution.ermapper_grids import read_ermapper_grid #name = project.outputname + '_0.0tide_notsunami' #name = project.outputname + '_0.0tide' name = project.outputname #sww2dem(name, basename_out = name + '_speed', sww2dem(name, basename_out = name + '_momentum', #quantity = 'stage-elevation', #Depth quantity = '(xmomentum**2 + ymomentum**2)**0.5', #Absolute momentum #quantity = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5', #Speed cellsize = 50, #cellsize = 500, #cellsize = 10, #easting_min = 465646, #easting_max = 480286, #northing_min = 7712865, #northing_max = 7724233, reduction = max, verbose = True, format = 'asc') #Check #data = read_ermapper_grid(name) #print 'Values from %s are in [%f, %f]' %(name, min(data.flat), max(data.flat))