********************************************************** 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/Tp/Tp2 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) 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 = .0952 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 = 10.504 dir = 270.000 duration = 720.000 Derived Trep = 9.605 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 16.5 percent complete. Average dt 0.276 seconds Time remaining 25 seconds waveparams: Reading from jonswap_2.txt ... XBeach reading from jonswap_2.txt Hm0 = 9.0000 Warning: value > recommended value of 5.0000 fp = .0952 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 = 10.504 dir = 270.000 duration = 720.000 Derived Trep = 9.605 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 28.4 percent complete. Average dt 0.223 seconds Time remaining 30 seconds Simulation 36.0 percent complete. Average dt 0.218 seconds Time remaining 42 seconds waveparams: Reading from jonswap_3.txt ... XBeach reading from jonswap_3.txt Hm0 = 9.0000 Warning: value > recommended value of 5.0000 fp = .0952 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 = 10.504 dir = 270.000 duration = 720.000 Derived Trep = 9.605 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 47.1 percent complete. Average dt 0.194 seconds Time remaining 24 seconds waveparams: Reading from jonswap_4.txt ... XBeach reading from jonswap_4.txt Hm0 = 9.0000 Warning: value > recommended value of 5.0000 fp = .0952 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 = 10.504 dir = 270.000 duration = 720.000 Derived Trep = 9.605 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 59.2 percent complete. Average dt 0.210 seconds Time remaining 17 seconds Simulation 69.5 percent complete. Average dt 0.213 seconds Time remaining 15 seconds waveparams: Reading from jonswap_5.txt ... XBeach reading from jonswap_5.txt Hm0 = 9.0000 Warning: value > recommended value of 5.0000 fp = .0952 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 = 10.504 dir = 270.000 duration = 720.000 Derived Trep = 9.605 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 80.7 percent complete. Average dt 0.193 seconds Time remaining 9 seconds Simulation 87.9 percent complete. Average dt 0.184 seconds Time remaining 8 seconds waveparams: Reading from jonswap_6.txt ... XBeach reading from jonswap_6.txt Hm0 = 9.0000 Warning: value > recommended value of 5.0000 fp = .0952 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 = 10.504 dir = 270.000 duration = 252.200 Derived Trep = 9.605 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 99.4 percent complete. Average dt 0.193 seconds Time remaining 0 seconds Duration : 31.442219000000001 seconds Timesteps : 0 Average dt : +Infinity seconds Unit speed : +Infinity seconds/1 End of program xbeach