source: trunk/anuga_work/development/2010-projects/anuga_1d/sww-sudi/parameter.py @ 8427

from scipy.special import jn
from scipy import pi, sqrt, linspace, pi, sin, cos
from config import g
from Numeric import zeros, Float
from rootsearch import *
from bisect_function import *

def j0(x):
    return jn(0.0, x)

def j1(x):
    return jn(1.0, x)

def j2(x):
    return jn(2.0, x)

def j3(x):
    return jn(3.0, x)

def jm1(x):
    return jn(-1.0, x)

def jm2(x):
    return jn(-2.0, x)


#DIMENSIONAL PARAMETERS
L = 5e4                 # Length of channel (m)
h_0 = 5e2               # Height at origin when the water is still
numb = 550#650              # initial Number of computational cells
cell_len = 100.0        # is actually 1.1*L/numb
Tp = 900.0#1020.0#3600.0#1020.0#900.0# Period of oscillation
a = 1.0#5.0#1.0#                                  # Amplitude at origin
               
points = zeros(numb+2,Float)
for i in range(numb+2):
        points[i] = i*cell_len - 0.5*cell_len
N = len(points)-1       # Number of computational cells


#DIMENSIONLESS PARAMETERS
eps = a/h_0
T = Tp*sqrt(g*h_0)/L
A = eps/j0(4.0*pi/T)