********************************************************** Welcome to XBeach revision 2467M date Fri Dec 9 16:38:40 CET 2011 URL: URL: https://svn.oss.deltares.nl/repos/xbeach/trunk ********************************************************** Simulation started: YYYYMMDD hh:mm:ss time zone (UTC) 20111212 10:13:03 +0100 running in: /u/heijer/sources/run01/D50/D504 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 = .0003 D90 = .0004 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) facDc = 1.0000 (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 = 4 nglobalvar: Will generate global output for variable:H nglobalvar: Will generate global output for variable:zs nglobalvar: Will generate global output for variable:zb nglobalvar: Will generate global output for variable:dzav 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 Creating netcdf variable: dzav Writing file definition. Setting up boundary conditions waveparams: Reading from jonswap_1.txt ... XBeach reading from jonswap_1.txt Hm0 = 9.0000 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 = 9.000 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 12.9 percent complete. Average dt 0.265 seconds Time remaining 34 seconds waveparams: Reading from jonswap_2.txt ... XBeach reading from jonswap_2.txt Hm0 = 9.0000 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 = 9.000 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 23.7 percent complete. Average dt 0.215 seconds Time remaining 35 seconds Simulation 28.2 percent complete. Average dt 0.176 seconds Time remaining 1 minutes and 20 seconds waveparams: Reading from jonswap_3.txt ... XBeach reading from jonswap_3.txt Hm0 = 9.0000 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 = 9.000 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 40.1 percent complete. Average dt 0.207 seconds Time remaining 25 seconds Simulation 49.3 percent complete. Average dt 0.165 seconds Time remaining 27 seconds waveparams: Reading from jonswap_4.txt ... XBeach reading from jonswap_4.txt Hm0 = 9.0000 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 = 9.000 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 57.4 percent complete. Average dt 0.188 seconds Time remaining 27 seconds Simulation 66.5 percent complete. Average dt 0.174 seconds Time remaining 19 seconds Simulation 73.3 percent complete. Average dt 0.166 seconds Time remaining 20 seconds waveparams: Reading from jonswap_5.txt ... XBeach reading from jonswap_5.txt Hm0 = 9.0000 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 = 9.000 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 82.4 percent complete. Average dt 0.163 seconds Time remaining 10 seconds Simulation 90.9 percent complete. Average dt 0.154 seconds Time remaining 5 seconds waveparams: Reading from jonswap_6.txt ... XBeach reading from jonswap_6.txt Hm0 = 9.0000 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 = 9.000 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 Duration : 35.504601999999998 seconds Timesteps : 0 Average dt : +Infinity seconds Unit speed : +Infinity seconds/1 End of program xbeach