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))