********************************************************** Welcome to XBeach revisionRange 2714 Not mixed revisions and with local modificatio ns date2012/02/20 16:45:55 URL:https://svn.oss.deltares.nl/repos/xbeach/trunk/src/xbeachlibrary ********************************************************** Simulation started: YYYYMMDD hh:mm:ss time zone (UTC) 20120411 11:59:23 +0200 General Input Module Reading input parameters: -------------------------------- Physical processes: XBeach reading fromparams.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) bchwiz =0 (no record found, default value used) -------------------------------- Grid parameters: gridform =xbeach (no record found, default value used) xori =-8308.5365 yori =0.0000 alfa =0.0000 nx =619 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 of180.0000 thetamax =315.0000 Warning: value > recommended value of180.0000 dtheta =90.0000 Warning: value > recommended value of20.0000 thetanaut =1 -------------------------------- Model time parameters: CFL =0.7000 (no record found, default value used) tstop =22000.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 hotstartflo =0 (no record found, default value used) -------------------------------- Wave boundary condition parameters: instat =jons bcfile =filelist.txt taper =100.0000 (no record found, default value used) nmax =0.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: nonhspectru =0 (no record found, default value used) random =0 fcutoff =0.0000 (no record found, default value used) nspr =0 (no record found, default value used) trepfac =0.0100 (no record found, default value used) sprdthr =0.0800 (no record found, default value used) oldwbc =0 (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) nc =1 (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_f = None specified disch_times = 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 =0.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) breakerdela =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 =0.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 =0.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 =0.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 =0.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 =0.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 =0.4000 (no record found, default value used) D50 =0.0002 D90 =0.0003 rhos =2650.0000 (no record found, default value used) dzg =0.1000 (no record found, default value used) dzg1 =0.1000 (no record found, default value used) dzg2 =0.1000 (no record found, default value used) dzg3 =0.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 =0.0060 (no record found, default value used) smax =-1.0000 (no record found, default value used) tsfac =0.1000 (no record found, default value used) facua =0.1000 (no record found, default value used) facSk =0.1000 (no record found, default value used) facAs =0.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 =0.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 =4000.0000 morstop =22000.0000 (no record found, default value used) wetslp =0.3000 (no record found, default value used) dryslp =1.0000 (no record found, default value used) hswitch =0.1000 (no record found, default value used) dzmax =0.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 rotate =1 (no record found, default value used) tstart =4000.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 =-1 (no record found, default value used) 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 =3 nmeanvar: Will generate mean, min, max and variance output for variable:H nmeanvar: Will generate mean, min, max and variance output for variable:zs nmeanvar: Will generate mean, min, max and variance output for variable:zb ncross =0 (no record found, default value used) outputforma =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 =0.0050 (no record found, default value used) umin =0.0000 (no record found, default value used) hmin =0.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 =0.1000 (no record found, default value used) -------------------------------- Bed update numerics parameters: frac_dz =0.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 =0.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 (0.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 parameterrt inparams.txt Unknown, unused or multiple statements of parameterdtbc inparams.txt ------------------------------------ Building Grid and Bathymetry ------------------------------------ Initializing ..... -------------------------------- Initializing spectral wave boundary conditions -------------------------------- NetCDF outputformat Creating netcdf variable:H Creating netcdf variable:zs Creating netcdf variable:zs0 Creating netcdf variable:zb Creating netcdf variable:hh Creating netcdf variable:u Creating netcdf variable:v Creating netcdf variable:ue Creating netcdf variable:ve Creating netcdf variable:urms Creating netcdf variable:Fx Creating netcdf variable:Fy Creating netcdf variable:ccg Creating netcdf variable:ceqsg Creating netcdf variable:ceqbg Creating netcdf variable:Susg Creating netcdf variable:Svsg Creating netcdf variable:E Creating netcdf variable:R Creating netcdf variable:D Creating netcdf variable:DR Creating netcdf variable:H_mean Creating netcdf variable:H_var Creating netcdf variable:H_min Creating netcdf variable:H_max Creating netcdf variable:zs_mean Creating netcdf variable:zs_var Creating netcdf variable:zs_min Creating netcdf variable:zs_max Creating netcdf variable:zb_mean Creating netcdf variable:zb_var Creating netcdf variable:zb_min Creating netcdf variable:zb_max Writing file definition. Setting up boundary conditions waveparams: Reading fromjonswap_1.txt ... XBeach reading fromjonswap_1.txt Hm0 =9.0000 Warning: value > recommended value of5.0000 fp =0.0833 fnyq =1.0000 dfj =0.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 =0.5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin1 of1 writing wave energy toE_jonswap_1.txt ... Apply symmetric energy distribution w.r.t mean wave direction file done Calculating flux at boundary writing long wave mass flux toq_jonswap_1.txt ... file done Boundary conditions complete, starting computation Simulation 1.4 percent complete. Average dt 0.136 seconds Time remaining 5 minutes and 49 seconds Simulation 2.8 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 44 seconds Simulation 4.3 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 35 seconds Simulation 5.7 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 34 seconds Simulation 7.1 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 24 seconds Simulation 8.5 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 26 seconds Simulation 9.9 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 16 seconds Simulation 11.4 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 13 seconds Simulation 12.8 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 6 seconds Simulation 14.2 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 4 seconds Simulation 15.6 percent complete. Average dt 0.137 seconds Time remaining 4 minutes and 56 seconds waveparams: Reading fromjonswap_2.txt ... XBeach reading fromjonswap_2.txt Hm0 =9.0000 Warning: value > recommended value of5.0000 fp =0.0833 fnyq =1.0000 dfj =0.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 =0.5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin1 of1 writing wave energy toE_jonswap_2.txt ... Apply symmetric energy distribution w.r.t mean wave direction file done Calculating flux at boundary writing long wave mass flux toq_jonswap_2.txt ... file done Simulation 17.0 percent complete. Average dt 0.137 seconds Time remaining 5 minutes and 3 seconds Simulation 18.4 percent complete. Average dt 0.137 seconds Time remaining 4 minutes and 53 seconds Simulation 19.6 percent complete. Average dt 0.136 seconds Time remaining 5 minutes and 21 seconds Simulation 20.8 percent complete. Average dt 0.136 seconds Time remaining 5 minutes and 35 seconds Simulation 21.9 percent complete. Average dt 0.136 seconds Time remaining 5 minutes and 49 seconds Simulation 23.0 percent complete. Average dt 0.136 seconds Time remaining 6 minutes and 0 seconds Simulation 24.0 percent complete. Average dt 0.136 seconds Time remaining 6 minutes and 9 seconds Simulation 25.1 percent complete. Average dt 0.136 seconds Time remaining 6 minutes and 10 seconds Simulation 25.9 percent complete. Average dt 0.135 seconds Time remaining 7 minutes and 30 seconds Simulation 26.7 percent complete. Average dt 0.136 seconds Time remaining 7 minutes and 26 seconds Simulation 28.0 percent complete. Average dt 0.136 seconds Time remaining 4 minutes and 49 seconds Simulation 29.3 percent complete. Average dt 0.136 seconds Time remaining 4 minutes and 25 seconds Simulation 30.6 percent complete. Average dt 0.136 seconds Time remaining 4 minutes and 20 seconds Simulation 32.0 percent complete. Average dt 0.136 seconds Time remaining 4 minutes and 17 seconds waveparams: Reading fromjonswap_3.txt ... XBeach reading fromjonswap_3.txt Hm0 =9.0000 Warning: value > recommended value of5.0000 fp =0.0833 fnyq =1.0000 dfj =0.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 =0.5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin1 of1 writing wave energy toE_jonswap_3.txt ... Apply symmetric energy distribution w.r.t mean wave direction file done Calculating flux at boundary writing long wave mass flux toq_jonswap_3.txt ... file done Simulation 33.3 percent complete. Average dt 0.136 seconds Time remaining 4 minutes and 15 seconds Simulation 34.6 percent complete. Average dt 0.132 seconds Time remaining 4 minutes and 12 seconds Simulation 35.9 percent complete. Average dt 0.134 seconds Time remaining 4 minutes and 10 seconds Simulation 37.0 percent complete. Average dt 0.119 seconds Time remaining 4 minutes and 28 seconds Simulation 38.4 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 52 seconds Simulation 39.7 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 43 seconds Simulation 41.0 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 44 seconds Simulation 42.4 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 35 seconds Simulation 43.7 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 33 seconds Simulation 45.0 percent complete. Average dt 0.133 seconds Time remaining 3 minutes and 30 seconds Simulation 46.2 percent complete. Average dt 0.124 seconds Time remaining 3 minutes and 42 seconds Simulation 47.5 percent complete. Average dt 0.130 seconds Time remaining 3 minutes and 24 seconds Simulation 48.9 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 11 seconds waveparams: Reading fromjonswap_4.txt ... XBeach reading fromjonswap_4.txt Hm0 =9.0000 Warning: value > recommended value of5.0000 fp =0.0833 fnyq =1.0000 dfj =0.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 =0.5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin1 of1 writing wave energy toE_jonswap_4.txt ... Apply symmetric energy distribution w.r.t mean wave direction file done Calculating flux at boundary writing long wave mass flux toq_jonswap_4.txt ... file done Simulation 50.2 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 6 seconds Simulation 51.5 percent complete. Average dt 0.136 seconds Time remaining 3 minutes and 3 seconds Simulation 52.8 percent complete. Average dt 0.134 seconds Time remaining 2 minutes and 60 seconds Simulation 54.2 percent complete. Average dt 0.136 seconds Time remaining 2 minutes and 54 seconds Simulation 55.4 percent complete. Average dt 0.125 seconds Time remaining 3 minutes and 0 seconds Simulation 56.7 percent complete. Average dt 0.136 seconds Time remaining 2 minutes and 44 seconds Simulation 58.1 percent complete. Average dt 0.136 seconds Time remaining 2 minutes and 36 seconds Simulation 59.4 percent complete. Average dt 0.136 seconds Time remaining 2 minutes and 31 seconds Simulation 60.8 percent complete. Average dt 0.136 seconds Time remaining 2 minutes and 25 seconds Simulation 62.1 percent complete. Average dt 0.136 seconds Time remaining 2 minutes and 21 seconds Simulation 63.4 percent complete. Average dt 0.136 seconds Time remaining 2 minutes and 16 seconds Simulation 64.7 percent complete. Average dt 0.125 seconds Time remaining 2 minutes and 24 seconds waveparams: Reading fromjonswap_5.txt ... XBeach reading fromjonswap_5.txt Hm0 =9.0000 Warning: value > recommended value of5.0000 fp =0.0833 fnyq =1.0000 dfj =0.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 =0.5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin1 of1 writing wave energy toE_jonswap_5.txt ... Apply symmetric energy distribution w.r.t mean wave direction file done Calculating flux at boundary writing long wave mass flux toq_jonswap_5.txt ... file done Simulation 66.0 percent complete. Average dt 0.131 seconds Time remaining 2 minutes and 12 seconds Simulation 67.3 percent complete. Average dt 0.134 seconds Time remaining 2 minutes and 5 seconds Simulation 68.6 percent complete. Average dt 0.136 seconds Time remaining 1 minutes and 56 seconds Simulation 70.0 percent complete. Average dt 0.136 seconds Time remaining 1 minutes and 51 seconds Simulation 71.3 percent complete. Average dt 0.136 seconds Time remaining 1 minutes and 47 seconds Simulation 72.5 percent complete. Average dt 0.123 seconds Time remaining 1 minutes and 56 seconds Simulation 73.9 percent complete. Average dt 0.136 seconds Time remaining 1 minutes and 37 seconds Simulation 75.2 percent complete. Average dt 0.136 seconds Time remaining 1 minutes and 33 seconds Simulation 76.7 percent complete. Average dt 0.134 seconds Time remaining 1 minutes and 19 seconds Simulation 78.0 percent complete. Average dt 0.123 seconds Time remaining 1 minutes and 21 seconds Simulation 79.5 percent complete. Average dt 0.134 seconds Time remaining 1 minutes and 9 seconds Simulation 81.0 percent complete. Average dt 0.135 seconds Time remaining 1 minutes and 4 seconds waveparams: Reading fromjonswap_6.txt ... XBeach reading fromjonswap_6.txt Hm0 =9.0000 Warning: value > recommended value of5.0000 fp =0.0833 fnyq =1.0000 dfj =0.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 =0.5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin1 of1 writing wave energy toE_jonswap_6.txt ... Apply symmetric energy distribution w.r.t mean wave direction file done Calculating flux at boundary writing long wave mass flux toq_jonswap_6.txt ... file done Simulation 82.5 percent complete. Average dt 0.136 seconds Time remaining 59 seconds Simulation 83.9 percent complete. Average dt 0.128 seconds Time remaining 57 seconds Simulation 85.4 percent complete. Average dt 0.133 seconds Time remaining 49 seconds Simulation 86.8 percent complete. Average dt 0.129 seconds Time remaining 46 seconds Simulation 88.3 percent complete. Average dt 0.136 seconds Time remaining 39 seconds Simulation 89.8 percent complete. Average dt 0.134 seconds Time remaining 34 seconds Simulation 91.2 percent complete. Average dt 0.126 seconds Time remaining 32 seconds Simulation 92.6 percent complete. Average dt 0.129 seconds Time remaining 26 seconds Simulation 94.1 percent complete. Average dt 0.136 seconds Time remaining 19 seconds Simulation 95.6 percent complete. Average dt 0.130 seconds Time remaining 15 seconds Simulation 97.1 percent complete. Average dt 0.133 seconds Time remaining 10 seconds waveparams: Reading fromjonswap_7.txt ... XBeach reading fromjonswap_7.txt Hm0 =9.0000 Warning: value > recommended value of5.0000 fp =0.0833 fnyq =1.0000 dfj =0.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 =80.000 Derived Trep =10.980 Changing dtbc in wave boundary conditions to satisfy Nyquist condition: New dtbc =0.5025 s. Calculating wave energy at boundary Calculating wave energy for theta bin1 of1 writing wave energy toE_jonswap_7.txt ... Apply symmetric energy distribution w.r.t mean wave direction file done Calculating flux at boundary writing long wave mass flux toq_jonswap_7.txt ... file done Simulation 98.5 percent complete. Average dt 0.132 seconds Time remaining 5 seconds Simulation 99.9 percent complete. Average dt 0.124 seconds Time remaining 0 seconds Duration : 371.079578700000 seconds Timesteps : 0 Average dt : Infinity seconds Unit speed : Infinity seconds/1 End of program xbeach