********************************************************** Welcome to XBeach revision 2429M date Tue Nov 29 13:07:22 CET 2011 URL: URL: https://svn.oss.deltares.nl/repos/xbeach/trunk ********************************************************** Simulation started: YYYYMMDD hh:mm:ss time zone (UTC) 20111202 10:06:47 +0100 running in: /u/heijer/sources/run00/Hsig/Hsig1 General Input Module Reading input parameters: -------------------------------- Physical processes: XBeach reading from params.txt swave = 1 (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) nonh = 0 (no record found, default value used) gwflow = 0 (no record found, default value used) q3d = 0 (no record found, default value used) swrunup = 0 (no record found, default value used) ships = 0 (no record found, default value used) -------------------------------- Grid parameters: gridform = xbeach (no record found, default value used) xori = -2930.6250 yori = 0.0000 alfa = 0.0000 nx = 233 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 thetamin = 225.0000 Warning: value > recommended value of 180.0000 thetamax = 315.0000 Warning: value > recommended value of 180.0000 dtheta = 90.0000 Warning: value > recommended value of 20.0000 thetanaut = 1 -------------------------------- Model time parameters: CFL = .7000 (no record found, default value used) tstop = 19260.0000 -------------------------------- Physical constants: rho = 1025.0000 (no record found, default value used) 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: instat = jons bcfile = filelist.txt taper = 100.0000 (no record found, default value used) nmax = .8000 (no record found, default value used) leftwave = neumann (no record found, default value used) rightwave = neumann (no record found, default value used) -------------------------------- Wave-spectrum boundary condition parameters: random = 0 fcutoff = 0.0000 (no record found, default value used) nspr = 0 (no record found, default value used) trepfac = .0100 (no record found, default value used) sprdthr = .0800 (no record found, default value used) oldwbc = 0 (no record found, default value used) newstatbc = 1 (no record found, default value used) correctHm0 = 1 (no record found, default value used) oldnyq = 0 (no record found, default value used) Tm01switch = 0 (no record found, default value used) wbcversion = 2 (no record found, default value used) -------------------------------- Flow boundary condition parameters: front = abs_2d (no record found, default value used) left = neumann (no record found, default value used) right = neumann (no record found, default value used) back = abs_2d (no record found, default value used) ARC = 1 (no record found, default value used) order = 2.0000 (no record found, default value used) carspan = 0 (no record found, default value used) freewave = 0 (no record found, default value used) epsi = 0.0000 (no record found, default value used) tidetype = velocity (no record found, default value used) -------------------------------- Tide boundary conditions: tideloc = 0 (no record found, default value used) zs0 = 5.0000 -------------------------------- 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) -------------------------------- Wave breaking parameters: break = roelvink2 (no record found, default value used) gamma = .5500 (no record found, default value used) alpha = 1.0000 (no record found, default value used) n = 10.0000 (no record found, default value used) gammax = 2.0000 (no record found, default value used) delta = 0.0000 (no record found, default value used) fw = 0.0000 (no record found, default value used) fwcutoff = 1000.0000 (no record found, default value used) breakerdelay = 1 (no record found, default value used) shoaldelay = 0 (no record found, default value used) facsd = 1.0000 (no record found, default value used) facrun = 1.0000 (no record found, default value used) -------------------------------- Roller parameters: roller = 1 (no record found, default value used) beta = .1000 (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) cats = 4.0000 (no record found, default value used) -------------------------------- Flow parameters: bedfriction = chezy (no record found, default value used) bedfricfile = None specified C = 55.0000 (no record found, default value used) nuh = .1000 (no record found, default value used) nuhfac = 1.0000 (no record found, default value used) nuhv = 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 = 0.0000 (no record found, default value used) -------------------------------- Wind parameters: rhoa = 1.2500 (no record found, default value used) Cd = .0020 (no record found, default value used) windfile = None specified windv = 0.0000 (no record found, default value used) windth = 270.0000 (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 (no record found, default value used) 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) -------------------------------- Sediment transport parameters: form = vanthiel_vanrijn (no record found, default value used) 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) facsl = 1.6000 (no record found, default value used) z0 = .0060 (no record found, default value used) smax = -1.0000 (no record found, default value used) tsfac = .1000 (no record found, default value used) facua = .1000 (no record found, default value used) facSk = .1000 (no record found, default value used) facAs = .1000 (no record found, default value used) turb = bore_averaged (no record found, default value used) Tbfac = 1.0000 (no record found, default value used) Tsmin = .5000 (no record found, default value used) lwt = 0 (no record found, default value used) betad = 1.0000 (no record found, default value used) sus = 1 (no record found, default value used) bed = 1 (no record found, default value used) bulk = 1 (no record found, default value used) -------------------------------- Morphology parameters: morfac = 5.0000 morfacopt = 1 (no record found, default value used) morstart = 1260.0000 morstop = 19260.0000 (no record found, default value used) wetslp = .3000 (no record found, default value used) dryslp = 1.0000 (no record found, default value used) hswitch = .1000 (no record found, default value used) dzmax = .0500 (no record found, default value used) struct = 0 (no record found, default value used) -------------------------------- Output variables: timings = 1 (no record found, default value used) tunits = None specified projection = None specified tstart = 1260.0000 tint = 1.0000 (no record found, default value used) tsglobal = None specified tintg = 100.0000 tspoints = None specified tintp = 1.0000 (no record found, default value used) tscross = None specified tintc = 1.0000 (no record found, default value used) tsmean = None specified tintm = 18000.0000 (no record found, default value used) nglobalvar = 3 nglobalvar: Will generate global output for variable:H nglobalvar: Will generate global output for variable:zs nglobalvar: Will generate global output for variable:zb npoints = 0 (no record found, default value used) nrugauge = 0 (no record found, default value used) npointvar = 0 (no record found, default value used) rugdepth = 0.0000 (no record found, default value used) nmeanvar = 0 (no record found, default value used) ncross = 0 (no record found, default value used) outputformat = netcdf ncfilename = None specified netcdf output to:xboutput.nc -------------------------------- Drifters parameters: drifterfile = None specified ndrifter = 0 (no record found, default value used) -------------------------------- Shipwaves parameters: shipfile = None specified -------------------------------- Wave numerics parameters: scheme = upwind_2 (no record found, default value used) -------------------------------- Flow numerics parameters: eps = .0050 (no record found, default value used) umin = 0.0000 (no record found, default value used) hmin = .2000 (no record found, default value used) secorder = 0 (no record found, default value used) oldhu = 0 (no record found, default value used) -------------------------------- Sediment transport numerics parameters: thetanum = 1.0000 (no record found, default value used) sourcesink = 0 (no record found, default value used) cmax = .1000 (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) -------------------------------- MPI parameters: mpiboundary = auto (no record found, default value used) -------------------------------- Finished reading input parameters -------------------------------- Warning: Setting rugdepth to minimum value greater than eps (.00500) Warning: 2D absorbing boundary condition [front=abs_2d] reduces to a 1D absorbing boundary condition [front=abs_1d] in 1D mode [ny=0] Warning: 2D absorbing boundary condition [back=abs_2d] reduces to a 1D absorbing boundary condition [back=abs_1d] in 1D mode [ny=0] Unknown, unused or multiple statements of parameter rt in params.txt Unknown, unused or multiple statements of parameter dtbc in params.txt ------------------------------------ Building Grid and Bathymetry ------------------------------------ Initializing ..... -------------------------------- Initializing spectral wave boundary conditions -------------------------------- NetCDF outputformat Creating netcdf variable: H Creating netcdf variable: zs Creating netcdf variable: zb Writing file definition. Setting up boundary conditions waveparams: Reading from jonswap_1.txt ... XBeach reading from jonswap_1.txt Hm0 = 6.7500 Warning: value > recommended value of 5.0000 fp = .0833 fnyq = 1.0000 dfj = .0050 (no record found, default value used) gammajsp = 3.3000 s = 20.0000 mainang = 270.0000 Input checked: Hm0 = 6.750 Tp = 12.005 dir = 270.000 duration = 720.000 Derived Trep = 10.980 Changing dtbc in wave boundary conditions to satisfy Nyquist condition: New dtbc = .5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin 1 of 1 writing wave energy to E_jonswap_1.txt ... file done Calculating flux at boundary Flux qx 1 of 1 Flux qy 1 of 1 Flux qtot 1 of 1 Flux eta 1 of 1 writing long wave mass flux to q_jonswap_1.txt ... file done Boundary conditions complete, starting computation Simulation 15.2 percent complete. Average dt 0.274 seconds Time remaining 28 seconds waveparams: Reading from jonswap_2.txt ... XBeach reading from jonswap_2.txt Hm0 = 6.7500 Warning: value > recommended value of 5.0000 fp = .0833 fnyq = 1.0000 dfj = .0050 (no record found, default value used) gammajsp = 3.3000 s = 20.0000 mainang = 270.0000 Input checked: Hm0 = 6.750 Tp = 12.005 dir = 270.000 duration = 720.000 Derived Trep = 10.980 Changing dtbc in wave boundary conditions to satisfy Nyquist condition: New dtbc = .5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin 1 of 1 writing wave energy to E_jonswap_2.txt ... file done Calculating flux at boundary Flux qx 1 of 1 Flux qy 1 of 1 Flux qtot 1 of 1 Flux eta 1 of 1 writing long wave mass flux to q_jonswap_2.txt ... file done Simulation 27.2 percent complete. Average dt 0.208 seconds Time remaining 30 seconds waveparams: Reading from jonswap_3.txt ... XBeach reading from jonswap_3.txt Hm0 = 6.7500 Warning: value > recommended value of 5.0000 fp = .0833 fnyq = 1.0000 dfj = .0050 (no record found, default value used) gammajsp = 3.3000 s = 20.0000 mainang = 270.0000 Input checked: Hm0 = 6.750 Tp = 12.005 dir = 270.000 duration = 720.000 Derived Trep = 10.980 Changing dtbc in wave boundary conditions to satisfy Nyquist condition: New dtbc = .5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin 1 of 1 writing wave energy to E_jonswap_3.txt ... file done Calculating flux at boundary Flux qx 1 of 1 Flux qy 1 of 1 Flux qtot 1 of 1 Flux eta 1 of 1 writing long wave mass flux to q_jonswap_3.txt ... file done Simulation 39.4 percent complete. Average dt 0.211 seconds Time remaining 25 seconds Simulation 48.7 percent complete. Average dt 0.173 seconds Time remaining 27 seconds Simulation 55.0 percent complete. Average dt 0.194 seconds Time remaining 36 seconds waveparams: Reading from jonswap_4.txt ... XBeach reading from jonswap_4.txt Hm0 = 6.7500 Warning: value > recommended value of 5.0000 fp = .0833 fnyq = 1.0000 dfj = .0050 (no record found, default value used) gammajsp = 3.3000 s = 20.0000 mainang = 270.0000 Input checked: Hm0 = 6.750 Tp = 12.005 dir = 270.000 duration = 720.000 Derived Trep = 10.980 Changing dtbc in wave boundary conditions to satisfy Nyquist condition: New dtbc = .5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin 1 of 1 writing wave energy to E_jonswap_4.txt ... file done Calculating flux at boundary Flux qx 1 of 1 Flux qy 1 of 1 Flux qtot 1 of 1 Flux eta 1 of 1 writing long wave mass flux to q_jonswap_4.txt ... file done Simulation 64.9 percent complete. Average dt 0.178 seconds Time remaining 18 seconds waveparams: Reading from jonswap_5.txt ... XBeach reading from jonswap_5.txt Hm0 = 6.7500 Warning: value > recommended value of 5.0000 fp = .0833 fnyq = 1.0000 dfj = .0050 (no record found, default value used) gammajsp = 3.3000 s = 20.0000 mainang = 270.0000 Input checked: Hm0 = 6.750 Tp = 12.005 dir = 270.000 duration = 720.000 Derived Trep = 10.980 Changing dtbc in wave boundary conditions to satisfy Nyquist condition: New dtbc = .5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin 1 of 1 writing wave energy to E_jonswap_5.txt ... file done Calculating flux at boundary Flux qx 1 of 1 Flux qy 1 of 1 Flux qtot 1 of 1 Flux eta 1 of 1 writing long wave mass flux to q_jonswap_5.txt ... file done Simulation 75.0 percent complete. Average dt 0.187 seconds Time remaining 13 seconds Simulation 84.1 percent complete. Average dt 0.162 seconds Time remaining 9 seconds waveparams: Reading from jonswap_6.txt ... XBeach reading from jonswap_6.txt Hm0 = 6.7500 Warning: value > recommended value of 5.0000 fp = .0833 fnyq = 1.0000 dfj = .0050 (no record found, default value used) gammajsp = 3.3000 s = 20.0000 mainang = 270.0000 Input checked: Hm0 = 6.750 Tp = 12.005 dir = 270.000 duration = 252.200 Derived Trep = 10.980 Changing dtbc in wave boundary conditions to satisfy Nyquist condition: New dtbc = .5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin 1 of 1 writing wave energy to E_jonswap_6.txt ... file done Calculating flux at boundary Flux qx 1 of 1 Flux qy 1 of 1 Flux qtot 1 of 1 Flux eta 1 of 1 writing long wave mass flux to q_jonswap_6.txt ... file done Simulation 94.6 percent complete. Average dt 0.179 seconds Time remaining 3 seconds Duration : 34.407769999999999 seconds Timesteps : 0 Average dt : +Infinity seconds Unit speed : +Infinity seconds/1 End of program xbeach