********************************************************** Welcome to XBeach version 1.23.5852 date01-09-2021 16:16:24 URL:http://svn.oss.deltares.nl/repos/xbeach/trunk/ ********************************************************** Simulation started: YYYYMMDD hh:mm:ss time zone (UTC) 20220302 11:08:07 +0100 General Input Module Reading input parameters: XBeach reading fromparams.txt -------------------------------- Backward compatibility: instat =jons -------------------------------- Physical processes: wavemodel =surfbeat cyclic =0 (no record found, default value used) swave =1 (no record found, default value used) single_dir =0 (no record found, default value used) lwave =1 (no record found, default value used) flow =1 (no record found, default value used) sedtrans =1 (no record found, default value used) morphology =1 (no record found, default value used) avalanching =1 (no record found, default value used) gwflow =0 (no record found, default value used) ships =0 (no record found, default value used) vegetation =0 (no record found, default value used) setbathy =0 (no record found, default value used) viscosity =1 (no record found, default value used) advection =1 (no record found, default value used) wind =0 (no record found, default value used) -------------------------------- Grid parameters: gridform =xbeach (no record found, default value used) xori =.0000 (no record found, default value used) yori =.0000 (no record found, default value used) alfa =.0000 nx =63 ny =0 posdwn =-1.0000 depfile =bed.dep vardx =1 dx =-1.0000 (no record found, default value used) dy =-1.0000 (no record found, default value used) xfile =x.grd yfile =y.grd nz =1 (no record found, default value used) thetamin =-90.0000 thetamax =90.0000 thetanaut =0 dtheta =180.0000 -------------------------------- Model time parameters: CFL =.9000 dtset =.0000 (no record found, default value used) tstop =162000.0000 maxdtfac =50.0000 (no record found, default value used) -------------------------------- Physical constants: rho =1025.0000 g =9.8100 (no record found, default value used) depthscale =1.0000 (no record found, default value used) -------------------------------- Initial conditions: zsinitfile = None specified -------------------------------- Wave boundary condition parameters: wbctype =parametric bcfile =waves.lst taper =100.0000 (no record found, default value used) nmax =.8000 (no record found, default value used) lateralwave =neumann (no record found, default value used) -------------------------------- Wave-spectrum boundary condition parameters: random =1 (no record found, default value used) fcutoff =.0000 (no record found, default value used) trepfac =.0100 (no record found, default value used) sprdthr =.0800 (no record found, default value used) Tm01switch =0 (no record found, default value used) nspectrumloc =1 (no record found, default value used) -------------------------------- Flow boundary condition parameters: front =abs_1d left =neumann right =neumann back =abs_1d ARC =1 (no record found, default value used) order =2.0000 (no record found, default value used) highcomp =0 (no record found, default value used) freewave =0 (no record found, default value used) epsi =-1.0000 tidetype =velocity (no record found, default value used) -------------------------------- Tide boundary conditions: tideloc =2 paulrevere =land (no record found, default value used) zs0file =tide.txt -------------------------------- Discharge boundary conditions: disch_loc_file = None specified disch_timeseries_file = None specified ndischarge =0 (no record found, default value used) ntdischarge =0 (no record found, default value used) beta =.1380 -------------------------------- Wave breaking parameters: break =roelvink2 (no record found, default value used) gamma =.5410 gammax =2.3640 alpha =1.2620 n =10.0000 (no record found, default value used) delta =.0000 (no record found, default value used) deltahmin =.1000 (no record found, default value used) fw =.0000 fwfile = None specified fwcutoff =1000.0000 (no record found, default value used) breakerdelay =1.0000 (no record found, default value used) -------------------------------- Roller parameters: roller =1 (no record found, default value used) rfb =0 (no record found, default value used) -------------------------------- Wave-current interaction parameters: wci =0 (no record found, default value used) hwci =.1000 (no record found, default value used) hwcimax =100.0000 (no record found, default value used) cats =4.0000 (no record found, default value used) -------------------------------- Flow parameters: bedfriction =chezy (no record found, default value used) bedfricfile = None specified bedfriccoef =55.0000 (no record found, default value used) droot =.5000 (no record found, default value used) dstem =.5000 (no record found, default value used) maxcf =.0400 (no record found, default value used) nuh =.1000 (no record found, default value used) nuhfac =1.0000 (no record found, default value used) smag =1 (no record found, default value used) -------------------------------- Coriolis force parameters: wearth =.0417 (no record found, default value used) lat =.0000 (no record found, default value used) -------------------------------- Sediment transport parameters: form =vanthiel_vanrijn (no record found, default value use d) waveform =vanthiel (no record found, default value used) sws =1 (no record found, default value used) lws =1 (no record found, default value used) BRfac =1.0000 (no record found, default value used) facua =.1000 (no record found, default value used) facSk =.3750 facAs =.1230 Tbfac =1.0000 (no record found, default value used) turb =bore_averaged (no record found, default value used) turbadv =none sus =1 (no record found, default value used) bed =1 (no record found, default value used) bulk =0 (no record found, default value used) facsl =.1500 (no record found, default value used) z0 =.0060 (no record found, default value used) smax =-1.0000 (no record found, default value used) bdslpeffmag =roelvink_total (no record found, default value used) bdslpeffini =none (no record found, default value used) bdslpeffdir =none (no record found, default value used) reposeangle =30.0000 (no record found, default value used) tsfac =.1000 (no record found, default value used) Tsmin =.5000 (no record found, default value used) facDc =1.0000 (no record found, default value used) jetfac =.1000 lwt =0 (no record found, default value used) betad =1.0000 (no record found, default value used) fallvelred =0 (no record found, default value used) dilatancy =0 (no record found, default value used) -------------------------------- Bed composition parameters: ngd =1 (no record found, default value used) nd =3 (no record found, default value used) por =.4000 (no record found, default value used) D50 =.0002 D90 =.0003 rhos =2650.0000 dzg =.1000 (no record found, default value used) dzg1 =.1000 (no record found, default value used) dzg2 =.1000 (no record found, default value used) dzg3 =.1000 (no record found, default value used) sedcal =1.0000 (no record found, default value used) ucrcal =1.0000 (no record found, default value used) -------------------------------- Morphology parameters: morfac =1.0000 (no record found, default value used) morfacopt =1 (no record found, default value used) morstart =.0000 morstop =162000.0000 (no record found, default value used) wetslp =.2600 dryslp =1.0000 (no record found, default value used) hswitch =.0000 (no record found, default value used) dzmax =.0500 (no record found, default value used) struct =1 ne_layer =nelayer.dep -------------------------------- Output variables: timings =1 (no record found, default value used) tunits = None specified tstart =.0000 tint =1.0000 (no record found, default value used) tsglobal = None specified tintg =1.0000 tspoints = None specified tintp =1.0000 tsmean = None specified tintm =162000.0000 nglobalvar =4 nglobalvar: Will generate global output for variable:zb nglobalvar: Will generate global output for variable:zs nglobalvar: Will generate global output for variable:sedero nglobalvar: Will generate global output for variable:H npoints =0 (no record found, default value used) nrugauge =1 npointvar =0 (no record found, default value used) Output runup gaugerugau001 xpoint:-1190.00 ypoint:.00 nrugdepth =1 (no record found, default value used) rugdepth =.0100 nmeanvar =4 nmeanvar: Will generate mean, min, max and variance output for variable:zb nmeanvar: Will generate mean, min, max and variance output for variable:zs nmeanvar: Will generate mean, min, max and variance output for variable:seder o nmeanvar: Will generate mean, min, max and variance output for variable:H outputformat =netcdf outputprecision =double (no record found, default value used) ncfilename = None specified netcdf output to:xboutput.nc remdryoutput =1 (no record found, default value used) -------------------------------- Output projection: projection = None specified rotate =1 (no record found, default value used) -------------------------------- Wave numerics parameters: scheme =warmbeam (no record found, default value used) snells =1 (no record found, default value used) -------------------------------- Flow numerics parameters: eps =.0050 (no record found, default value used) eps_sd =.5000 (no record found, default value used) umin =.0000 (no record found, default value used) hmin =.2000 (no record found, default value used) secorder =0 (no record found, default value used) -------------------------------- Sediment transport numerics parameters: thetanum =1.0000 (no record found, default value used) cmax =.1000 (no record found, default value used) dtlimTs =5.0000 (no record found, default value used) -------------------------------- Bed update numerics parameters: frac_dz =.7000 (no record found, default value used) nd_var =2 (no record found, default value used) split =1.0100 (no record found, default value used) merge =.0100 (no record found, default value used) nTrepavaltime =1.0000 (no record found, default value used) -------------------------------- Finished reading input parameters -------------------------------- Unknown, unused or multiple statements of parameter%MORPHOLOGY inparams.txt Unknown, unused or multiple statements of parameter%SEDTRANS inparams.txt Unknown, unused or multiple statements of parameter%SNELLS inparams.txt Unknown, unused or multiple statements of parameter%SWRUNUP inparams.txt Unknown, unused or multiple statements of parameterFACRUN inparams.txt Unknown, unused or multiple statements of parameterCF inparams.txt ------------------------------------ Building Grid and Bathymetry ------------------------------------ Initializing ..... readtide: reading tide time series fromtide.txt ... -------------------------------- Initializing spectral wave boundary conditions -------------------------------- Runup gauge at grid line iy=1 NetCDF outputformat Creating netcdf variable:zb Creating netcdf variable:zs Creating netcdf variable:sedero Creating netcdf variable:H Creating netcdf variable:point_xz Creating netcdf variable:point_yz Creating netcdf variable:point_zs Creating netcdf variable:zb_mean Creating netcdf variable:zb_var Creating netcdf variable:zb_min Creating netcdf variable:zb_max Creating netcdf variable:zs_mean Creating netcdf variable:zs_var Creating netcdf variable:zs_min Creating netcdf variable:zs_max Creating netcdf variable:sedero_mean Creating netcdf variable:sedero_var Creating netcdf variable:sedero_min Creating netcdf variable:sedero_max Creating netcdf variable:H_mean Creating netcdf variable:H_var Creating netcdf variable:H_min Creating netcdf variable:H_max Writing file definition. Writing point vars. Setting up boundary conditions -------------------------------- Calculating spectral wave boundary conditions -------------------------------- Reading spectrum at location1 waveparams: Reading from:waves1.bnd ... XBeach reading fromwaves1.bnd nmodal =1 (no record found, default value used) Hm0 =9.8900 Warning: value > recommended value of5.0000 fp =.0597 Warning: value < recommended value of.0625 gammajsp =3.3000 s =10000.0000 Warning: value > recommended value of1000 .0000 tma =0 (no record found, default value used) mainang =270.0000 fnyq =.3000 (no record found, default value used) dfj =.0015 (no record found, default value used) Interpreting spectrum at location1 Values calculated from interpolated spectrum: Hm0 =9.89 m Trep =15.31 s Mean dir =270.28 degN Overall Trep from all spectra calculated:15.31 s Calculating Fourier components 5.0% done 10.0% done 15.0% done 20.0% done 25.0% done 30.0% done 35.0% done 40.0% done 45.0% done 50.0% done 55.0% done 60.0% done 65.0% done 70.0% done 75.0% done 80.0% done 85.0% done 90.0% done 95.0% done 100.0% done Wave energy outside computational domain at offshore boundary:.0% Calculating short wave time series for theta bin1 of1 Calculating wave energy envelope at boundary. 100.0% done Wave bin1 contains100.00% of total energy Writing wave energy toE_series00001.bcf ... file done Calculating primary wave interaction 5.0% done 10.1% done 15.1% done 20.1% done 25.1% done 30.2% done 35.2% done 40.2% done 45.2% done 50.3% done 55.3% done 60.3% done 65.3% done 70.4% done 75.4% done 80.4% done 85.4% done 90.5% done 95.5% done 100.0% done Warning: shallow water so long wave variance is reduced using par%nmax Calculating flux at boundary 100.0% done Writing long wave mass flux toq_series00001.bcf ... file done -------------------------------- Spectral wave boundary conditions complete -------------------------------- Boundary conditions complete, starting computation Simulation 2.1 percent complete. Average dt 1.000 seconds Time remaining 3 minutes and 57 seconds -------------------------------- Calculating spectral wave boundary conditions -------------------------------- Reading spectrum at location1 waveparams: Reading from:waves1.bnd ... nmodal =1 (no record found, default value used) Hm0 =9.8900 Warning: value > recommended value of5.0000 fp =.0597 Warning: value < recommended value of.0625 gammajsp =3.3000 s =10000.0000 Warning: value > recommended value of1000 .0000 tma =0 (no record found, default value used) mainang =270.0000 fnyq =.3000 (no record found, default value used) dfj =.0015 (no record found, default value used) Interpreting spectrum at location1 Values calculated from interpolated spectrum: Hm0 =9.89 m Trep =15.31 s Mean dir =270.28 degN Overall Trep from all spectra calculated:15.31 s Calculating Fourier components 5.0% done 10.0% done 15.0% done 20.0% done 25.0% done 30.0% done 35.0% done 40.0% done 45.0% done 50.0% done 55.0% done 60.0% done 65.0% done 70.0% done 75.0% done 80.0% done 85.0% done 90.0% done 95.0% done 100.0% done Wave energy outside computational domain at offshore boundary:.0% Calculating short wave time series for theta bin1 of1 Calculating wave energy envelope at boundary. 100.0% done Wave bin1 contains100.00% of total energy Writing wave energy toE_series00002.bcf ... file done Calculating primary wave interaction 5.0% done 10.1% done 15.1% done 20.1% done 25.1% done 30.2% done 35.2% done 40.2% done 45.2% done 50.3% done 55.3% done 60.3% done 65.3% done 70.4% done 75.4% done 80.4% done 85.4% done 90.5% done 95.5% done 100.0% done Warning: shallow water so long wave variance is reduced using par%nmax Calculating flux at boundary 100.0% done Writing long wave mass flux toq_series00002.bcf ... file done -------------------------------- Spectral wave boundary conditions complete --------------------------------