********************************************************** 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/Hsig/Hsig5 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 = 11.2500 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 = 11.250 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 14.8 percent complete. Average dt 0.254 seconds Time remaining 29 seconds waveparams: Reading from jonswap_2.txt ... XBeach reading from jonswap_2.txt Hm0 = 11.2500 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 = 11.250 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 25.0 percent complete. Average dt 0.187 seconds Time remaining 37 seconds Simulation 30.7 percent complete. Average dt 0.168 seconds Time remaining 1 minutes and 1 seconds waveparams: Reading from jonswap_3.txt ... XBeach reading from jonswap_3.txt Hm0 = 11.2500 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 = 11.250 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 41.3 percent complete. Average dt 0.187 seconds Time remaining 28 seconds Simulation 49.8 percent complete. Average dt 0.148 seconds Time remaining 29 seconds waveparams: Reading from jonswap_4.txt ... XBeach reading from jonswap_4.txt Hm0 = 11.2500 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 = 11.250 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 58.5 percent complete. Average dt 0.181 seconds Time remaining 24 seconds Simulation 67.6 percent complete. Average dt 0.157 seconds Time remaining 18 seconds Simulation 73.6 percent complete. Average dt 0.149 seconds Time remaining 22 seconds waveparams: Reading from jonswap_5.txt ... XBeach reading from jonswap_5.txt Hm0 = 11.2500 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 = 11.250 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.3 percent complete. Average dt 0.151 seconds Time remaining 10 seconds Simulation 90.8 percent complete. Average dt 0.144 seconds Time remaining 5 seconds waveparams: Reading from jonswap_6.txt ... XBeach reading from jonswap_6.txt Hm0 = 11.2500 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 = 11.250 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 : 39.754956000000000 seconds Timesteps : 0 Average dt : +Infinity seconds Unit speed : +Infinity seconds/1 End of program xbeach