np.random.seed(opts.seed)
if opts.RandomWaves == True:
phi = 2 * np.pi * np.random.rand(opts.N)
Tend = opts.Ntotalwaves * opts.Tp / 1.1
wave = wt.RandomWavesFast(Tstart=opts.Tstart,
Tend=Tend,
x0=opts.x0,
Tp=opts.Tp,
Hs=opts.Hs,
mwl=opts.mwl,
depth=opts.depth,
waveDir=opts.waveDir,
g=opts.g,
N=opts.N,
bandFactor=opts.bandFactor,
spectName=opts.spectName,
spectral_params=opts.spectral_params,
phi=phi,
Lgen=opts.Lgen,
Nwaves=opts.Nwaves,
Nfreq=opts.Nfreq,
checkAcc=True,
fast=True)
Duration = Tend
wave_length = wave.wavelength
else:
wave = wt.NewWave(Tp=opts.Tp,
Tend = 2.*opts.T
x0 = np.array([0., 0. ,0. ])
waveDir = np.array(opts.wave_dir) #[1.,0.,0.]
g = np.array(opts.g)
N = 2000
phi = np.loadtxt("phases.txt")
Lgen = np.array([opts.tank_sponge[0], 0., 0.])
wave = wt.RandomWavesFast(0.,
Tend,
x0,
Tp,
Hs,
mwl,#m significant wave height
depth , #m depth
waveDir,
g, #peak frequency
N,
bandFactor, #accelerationof gravity
spectName ,# random words will result in error and return the available spectra
spectral_params, #JONPARAMS = {"gamma": 3.3, "TMA":True,"depth": depth}
phi,
Lgen = Lgen,
Nfreq = 32,
Nwaves = 15,
checkAcc = True)
# tank options
waterLevel = opts.water_level
tank_dim = opts.tank_dim
tank_sponge = opts.tank_sponge
# ----- DOMAIN ----- #