# Generated on 06/22/2017 14:24:20 # Deltares, FM-Suite DFlowFM Model Version 1.3.0.37667, DFlow FM Version 1.1.192.51191 [model] Program = D-Flow FM Version = 1.1.192.51191 MDUFormatVersion = 1.02 GuiVersion = 1.3.0.37667 AutoStart = 0 # Autostart simulation after loading MDU or not (0=no, 1=autostart, 2=autostartstop). [geometry] NetFile = t07_net.nc # *_net.nc BathymetryFile = bathymetry0.xyz # *.xyb DryPointsFile = WaterLevIniFile = # Initial water levels sample file *.xyz LandBoundaryFile = # Only for plotting ThinDamFile = # *_thd.pli, Polyline(s) for tracing thin dams. FixedWeirFile = StructureFile = VertplizFile = # *_vlay.pliz), = pliz with x,y, Z, first Z =nr of layers, second Z = laytyp ProflocFile = # *_proflocation.xyz) x,y,z, z = profile refnumber ProfdefFile = # *_profdefinition.def) definition for all profile nrs ProfdefxyzFile = # *_profdefinition.def) definition for all profile nrs Uniformwidth1D = 0.3 # Uniform width for 1D profiles not specified bij profloc ManholeFile = # *... WaterLevIni = 0.4 # Initial water level Bedlevuni = -5 # Uniform bottom level, (only if bedlevtype>=3, used at missing z values in netfile Bedslope = 0 BedlevType = 1 # 1 : Bottom levels at waterlevel cells (=flow nodes), like tiles xz, yz, bl , bob = max(bl left, bl right) Blmeanbelow = -999 Blminabove = -999 # 2 : Bottom levels at velocity points (=flow links), xu, yu, blu, bob = blu, bl = lowest connected link # 3 : Bottom levels at velocity points (=flow links), using mean network levels xk, yk, zk bl = lowest connected link # 4 : Bottom levels at velocity points (=flow links), using min network levels xk, yk, zk bl = lowest connected link # 5 : Bottom levels at velocity points (=flow links), using max network levels xk, yk, zk bl = lowest connected link PartitionFile = # *_part.pol, polyline(s) x,y AngLat = 0 # Angle of latitude S-N (deg), 0=no Coriolis AngLon = 0 # Angle of longitude E-W (deg), 0=Greenwich Conveyance2D = -1 # -1:R=HU,0:R=H, 1:R=A/P, 2:K=analytic-1D conv, 3:K=analytic-2D conv Nonlin2D = 0 # Non-linear 2D volumes, only icm ibedlevtype = 3 Sillheightmin = 0.5 Makeorthocenters = 0 # 1=yes, 0=no switch from circumcentres to orthocentres in geominit Dcenterinside = 1 # limit cell center; 1.0:in cell <-> 0.0:on c/g Bamin = 1E-06 OpenBoundaryTolerance= 3 RenumberFlowNodes = 1 Kmx = 0 Layertype = 1 Numtopsig = 0 SigmaGrowthFactor = 1 ThindykeFile = # *_tdk.pli, Polyline(s) x,y,z, z = thin dyke top levels [numerics] CFLMax = 0.7 # Max. Courant nr. AdvecType = 33 # Adv type, 0=no, 1= Wenneker, qu-udzt, 2=1, q(uio-u), 3=Perot q(uio-u), 4=Perot q(ui-u), 5=Perot q(ui-u) without itself TimeStepType = 2 # 0=only transport, 1=transport + velocity update, 2=full implicit step_reduce, 3=step_jacobi, 4=explicit Limtyphu = 0 Limtypmom = 4 # Limiter type for cell center advection velocity, 0=no, 1=minmod,2=vanLeer,3=Kooren,4=Monotone Central Limtypsa = 0 # Limiter type for salinity transport, 0=no, 1=minmod,2=vanLeer,3=Kooren,4=Monotone Central TransportMethod = 1 Vertadvtypsal = 5 Icgsolver = 4 # Solver type , 1 = sobekGS_OMP, 2 = sobekGS_OMPthreadsafe, 3 = sobekGS, 4 = sobekGS + Saadilud, 5 = parallel/global Saad, 6 = parallel/Petsc, 7 = parallel/GS Maxdegree = 6 FixedWeirScheme = 0 FixedWeirContraction= 1 FixedWeirfrictscheme= 1 Fixedweirtopwidth = 3 Fixedweirtopfrictcoef= -999 Fixedweirtalud = 0.25 Izbndpos = 0 Tlfsmo = 0 # Fourier smoothing time on waterlevel boundaries (s) Slopedrop2D = 0.3 # Apply droplosses only if local bottom slope > Slopedrop2D, <=0 =no droplosses Chkadvd = 0.1 Teta0 = 0.55 Qhrelax = 0.01 Jbasqbnddownwindhs= 0 cstbnd = 0 # Delft-3D type velocity treatment near boundaries for small coastal models (1) or not (0) Maxitverticalforestersal= 100 Maxitverticalforestertem= 0 Jaorgsethu = 1 Turbulencemodel = 3 Turbulenceadvection= 3 AntiCreep = 0 Maxwaterleveldiff = 0 Maxvelocitydiff = 0 Epshu = 0.0001 SobekDFM_umin = 0 [physics] UnifFrictCoef = 45 # Uniform friction coefficient, 0=no friction UnifFrictType = 0 # 0=Chezy, 1=Manning, 2=White Colebrook, 3=z0 etc UnifFrictCoef1D = 45 # Uniform friction coefficient in 1D links, 0=no friction UnifFrictCoefLin = 0 # Uniform linear friction coefficient for ocean models (m/s), 0=no Umodlin = 0 Vicouv = 1E-06 # Uniform horizontal eddy viscosity (m2/s) Dicouv = 1 # Uniform horizontal eddy diffusivity (m2/s) Vicoww = 5E-05 Dicoww = 5E-05 Vicwminb = 0 Smagorinsky = 0 # Add Smagorinsky horizontal turbulence : vicu = vicu + ( (Smagorinsky*dx)**2)*S, e.g. 0.1 Elder = 0 # Add Elder contribution : vicu = vicu + Elder*kappa*ustar*H/6), e.g. 1.0 Irov = 0 wall_ks = 0 # Nikuradse roughness for side walls, wall_z0=wall_ks/30 Rhomean = 1000 # Average water density (kg/m3) Idensform = 1 Ag = 9.81 # Gravitational acceleration TidalForcing = 1 # Tidal forcing (0=no, 1=yes) (only for jsferic == 1) Doodsonstart = 55.565 Doodsonstop = 375.575 Doodsoneps = 0 Salinity = 0 # Include salinity, (0=no, 1=yes) InitialSalinity = 0 Sal0abovezlev = -999 DeltaSalinity = -999 Backgroundsalinity= 0 Temperature = 0 # Include temperature, (0=no, 1=only transport, 5=heat flux model (5) of D3D) InitialTemperature= 6 Secchidepth = 1 Stanton = -1 Dalton = -1 Backgroundwatertemperature= 20 SecondaryFlow = 0 EffectSpiral = 0 BetaSpiral = 0 #BedformFile = trench.bfm [wind] ICdtyp = 2 # ( ),1=const, 2=S&B 2 breakpoints, 3= S&B 3 breakpoints, 4=Charnock constant Cdbreakpoints = 0.0025 0.0025 # ( ), e.g. 0.00063 0.00723 Windspeedbreakpoints= 0 100 # (m/s), e.g. 0.0 100.0 Rhoair = 1.2 PavBnd = 0 PavIni = 0 [waves] Wavemodelnr = 0 WaveNikuradse = 0.01 Rouwav = FR84 Gammax = 1 [time] RefDate = 20000218 # Reference date (yyyymmdd) Tzone = 0 # Data Sources in GMT are interrogated with time in minutes since refdat-Tzone*60 DtUser = 0.25 # User timestep in seconds (interval for external forcing update & his/map output) DtNodal = DtMax = 0.25 # Max timestep in seconds DtInit = 1 # Initial timestep in seconds Tunit = M # Time units in MDU (H, M or S) TStart = 0 # Start time w.r.t. RefDate (in TUnit) TStop = 30 # Stop time w.r.t. RefDate (in TUnit) AutoTimestep = 5 # Use CFL timestep limit or not (1/0) [restart] RestartFile = # Restart file, only from netcdf-file, hence: either *_rst.nc or *_map.nc RestartDateTime = 20000218 # Restart time (YYYYMMDDHHMMSS), only relevant in case of restart from *_map.nc [external forcing] ExtForceFile = t07.ext # *.ext ExtForceFileNew = [trachytopes] TrtRou = N TrtDef = TrtL = TrtDt = 60 [output] OutputDir = dflowfmoutput # Output directory of map-, his-, rst-, dat- and timings-files, default: DFM_OUTPUT_. Set to . for no dir/current dir. FlowGeomFile = # *_flowgeom.nc Flow geometry file in NetCDF format. ObsFile = # *.xyn Coords+name of observation stations. CrsFile = t07_crs.pli # *_crs.pli Polyline(s) definining cross section(s). HisFile = HisInterval = 10 # History output, given as "interval" "start period" "end period" (s) XLSInterval = # Interval (s) between XLS history MapFile = MapInterval = 900 # Map file output, given as "interval" "start period" "end period" (s) RstInterval = 0 # Restart file output, given as "interval" "start period" "end period" (s) S1incinterval = WaqInterval = 0 # Interval (in s) between Delwaq file outputs StatsInterval = # Interval (in s) between simulation statistics output. Writebalancefile = 0 TimingsInterval = # Timings output interval MapFormat = 4 # Map file format, 1: netCDF, 2: Tecplot, 3: netCFD and Tecplot Wrihis_balance = 1 Wrihis_structure_gen= 1 Wrihis_structure_dam= 1 Wrihis_structure_pump= 1 Wrihis_structure_gate= 1 Wrimap_waterlevel_s0= 1 Wrimap_waterlevel_s1= 1 Wrimap_velocity_component_u0= 1 Wrimap_velocity_component_u1= 1 Wrimap_velocity_vector= 1 Wrimap_upward_velocity_component= 0 Wrimap_density_rho= 1 Wrimap_horizontal_viscosity_viu= 1 Wrimap_horizontal_diffusivity_diu= 1 Wrimap_flow_flux_q1= 1 Wrimap_spiral_flow= 1 Wrimap_numlimdt = 1 Wrimap_taucurrent = 1 Wrimap_chezy = 1 Wrimap_turbulence = 1 Wrimap_wind = 1 Wrimap_heat_fluxes= 0 Richardsononoutput= 0 # 1=yes,0=no MapOutputTimeVector= # File (.mpt) containing fixed map output times (s) w.r.t. RefDate FullGridOutput = 0 # 0:compact, 1:full time-varying grid data EulerVelocities = 0 TimeSplitInterval = 0X # Time splitting interval, after which a new output file is started. value+unit, e.g. '1 M', valid units: Y,M,D,h,m,s.