********************************************************** 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 11:23:17 +0100 running in: /u/heijer/sources/run02/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 = 0 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 = 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.6 percent complete. Average dt 0.259 seconds Time remaining 27 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 26.3 percent complete. Average dt 0.178 seconds Time remaining 34 seconds Simulation 36.1 percent complete. Average dt 0.168 seconds Time remaining 33 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 46.0 percent complete. Average dt 0.164 seconds Time remaining 27 seconds Simulation 55.9 percent complete. Average dt 0.162 seconds Time remaining 22 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 65.8 percent complete. Average dt 0.163 seconds Time remaining 17 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.6 percent complete. Average dt 0.167 seconds Time remaining 14 seconds Simulation 82.2 percent complete. Average dt 0.156 seconds Time remaining 13 seconds Simulation 91.0 percent complete. Average dt 0.143 seconds Time remaining 5 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 Duration : 38.716113999999997 seconds Timesteps : 0 Average dt : +Infinity seconds Unit speed : +Infinity seconds/1 End of program xbeach