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