source: anuga_work/production/new_south_wales/batemans_bay/run_up.py @ 7586

""" This code implements run up solutions to the shallow water wave
equations.

References:
Madsen and Fuhrman 2008. Run-up of tsunamis and long waves
in terms of surf-similarity. Coastal Engineering, 55, p209.

Creator: Jonathan Griffin, Geoscience Australia
Created: 10 November 2009
"""

def Madsen(depth, slope, height, period, gravity = 9.81):

    freq = (2*pi)/period
    a = (depth*power(freq,2.0))/(gravity*power(slope,2.0))
    #print 'a',a
    b = power(a,0.25)
    #print 'b',b
    run_up_madsen =  2*power(np.pi,0.5)*b*height

    R_break = gravity*(power(slope,2.0))/(power(freq,2.0))
    run_up = min(run_up_madsen, R_break)
    return run_up

def surf_sim(depth, slope, height, period, gravity = 9.81):
    wavelength = period*sqrt(gravity*depth)
    print 'wavelength', wavelength
    surf_sim_param = slope/(sqrt(height/wavelength))
    print 'surf_sim_param', surf_sim_param
    #print 'height',height
##    print 'depth', depth
    a = 2*(power(pi,(0.75)))
    HA = height/depth
    #print HA
    b = power(HA,-0.25)
    c =power(surf_sim_param,(-0.5))
##    print 'a', a
##    print 'b' ,b
##    print 'c', c
    R = height*a*b*c
    #R = height*2*(power(pi,(3/4)))*(power((height/depth),(-1/4)))*(power(surf_sim_param,(-1/2)))
    print 'R',R
    R_break = height*(1.0/pi)*power(surf_sim_param,2.0)
    print 'R_break',R_break
    run_up = min(R, R_break)
##    print 'run_up',run_up
    return run_up
    




def energy_conserv(crest_integrated_energy,theta = 3.,alpha = 2.,gravity = 9.81):
    E = crest_integrated_energy
    print 'crest_integrated_energy', crest_integrated_energy
    a = 1/tan(alpha*pi/180.)
    b = 1/tan(theta*pi/180.)                                
 
    F = a-b
#    print 'F',F
    R3 = 6.*E/(1000.*gravity*F)
#    print 'R3', R3
    R = power(R3,(1./3.))

    return R

import numpy as np
from numpy import sin,tan,cos,power,pi, sqrt
if __name__ == '__main__':
    sys.exit(main())