# Generated on 12/10/2015 11:43:06 # Deltares, FM-Suite DFlowFM Model Version 1.0.0.31438, DFlow FM Version 1.1.148.41699 [model] Program = D-Flow FM Version = 1.1.148.41699 MDUFormatVersion = 1.02 GuiVersion = 1.0.0.31438 AutoStart = 0 # Autostart simulation after loading MDU or not (0=no, 1=autostart, 2=autostartstop). [geometry] NetFile = scenarioA_04_wbthy_net.nc# *_net.nc BathymetryFile = # *.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 = ProflocFile = # *_proflocation.xyz) x,y,z, z = profile refnumber ProfdefFile = # *_profdefinition.def) definition for all profile nrs ProfdefxyzFile = Uniformwidth1D = 2 ManholeFile = # *... WaterLevIni = -0.22 # Initial water level Bedlevuni = -5 # Uniform bottom level, (only if bedlevtype>=3, used at missing z values in netfile Bedslope = 0 # bedslopeination, sets zk = bedlevuni + x*bedslope ans sets zbndz = xbndz*bedslope BedlevType = 3 Blmeanbelow = -999 Blminabove = -999 PartitionFile = # *_part.pol, polyline(s) x,y AngLat = 52 # Angle of latitude (deg), 0=no Coriolis AngLon = 0 Conveyance2D = 3 # -1:R=HU,0:R=H, 1:R=A/P, 2:K=analytic-1D conv, 3:K=analytic-2D conv Nonlin2D = 0 Sillheightmin = 0.5 # Weir treatment only if both sills larger than this value (m) 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 = 0 Kmx = 0 # Max nr of vertical layers 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 = 3 # 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 Limtyphu = 0 # Limiter type for waterdepth in continuity eq., 0=no, 1=minmod,2=vanLeer,3=Kooren,4=Monotone Central Limtypmom = 4 # Limiter type for cell center advection velocity, 0=no, 1=minmod,2=vanLeer,3=Kooren,4=Monotone Central Limtypsa = 4 # 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 # Maximum degree in Gauss elimination FixedWeirScheme = 0 FixedWeirContraction= 1 FixedWeirfrictscheme= 1 Fixedweirtopwidth = 3 Fixedweirtopfrictcoef= -999 Fixedweirtalud = 0.25 Izbndpos = 0 # Position of z boundary, 0=D3Dflow, 1=on net boundary, 2 = on specified polyline Tlfsmo = 3600 # Fourier smoothing time on waterlevel boundaries (s) Slopedrop2D = 0 # Apply droplosses only if local bottom slope > Slopedrop2D Chkadvd = 0.1 # Check advection terms if depth < chkadvdp, => less setbacks Teta0 = 0.55 # Teta of time integration, 0.5 < Teta < 1d0 Qhrelax = 0.01 Jbasqbnddownwindhs= 0 # 0 : original hu on qbnd, 1 = downwind hs on qbnd cstbnd = 0 Maxitverticalforestersal= 100 Maxitverticalforestertem= 0 Jaorgsethu = 1 Turbulencemodel = 1 Turbulenceadvection= 0 AntiCreep = 0 Maxwaterleveldiff = 0 Maxvelocitydiff = 0 Epshu = 0.0001 CFLWaveFrac = 0.1 # Wave velocity fraction, total courant vel = u + cflw*wavevelocity Icoriolistype = # Coriolis weigthing, icorio 0, 4 or 5 (4 : v = tang. comp of total velocity at u-point, continuity weighted ucnxq, ucnyq = default so far, 5 : 4, scaling down coriolis force below trshcorio) Thindykescheme = # 0 = no, 1 = compact stencil, 2 = whole tile lifted, full subgrid weir + factor Thindykecontraction= # flow width = flow width*thindykecontraction Lincontin = # 0 = no, 1 = yes linear continuity Jaembed1D = # 1 : use embedded 1d channels, first run: Operations: Embed 1D channels Maxitverticalforester= # 0 : no vertical filter, > 0 = Max nr of iterations [physics] UnifFrictCoef = 0.023 # Uniform friction coefficient, 0=no friction UnifFrictType = 1 # 0=Chezy, 1=Manning, 2=White Colebrook, 3=z0 etc UnifFrictCoef1D = 0.023 UnifFrictCoefLin = 0 Umodlin = 0 # linear friction umod, ifrctyp 4,5,6 Vicouv = 1 # Uniform horizontal eddy viscosity (m2/s) Dicouv = 1 # Uniform horizontal eddy diffusivity (m2/s) Vicoww = 0 # Uniform vertical eddy viscosity (m2/s) Dicoww = 0 # Uniform vertical eddy diffusivity (m2/s) 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 # 0=free slip, 1 = partial slip using wall_ks wall_ks = 0 # Nikuradse roughness for side walls, wall_z0=wall_ks/30 Rhomean = 1023 # Average water density (kg/m3) Idensform = 1 Ag = 9.81 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 # Initial salinity concentration (ppt) Sal0abovezlev = -999 DeltaSalinity = -999 Backgroundsalinity= 30 Temperature = 0 InitialTemperature= 6 Secchidepth = 1 Stanton = -1 Dalton = -1 Backgroundwatertemperature= 6 SecondaryFlow = 0 EffectSpiral = 0 BetaSpiral = 0 [wind] ICdtyp = 3 # ( ), Cd = const, 2=S&B 2 breakpoints, 3= S&B 3 breakpoints Cdbreakpoints = 0.00063 0.00723 0.00723# ( ), e.g. 0.00063 0.00723 Windspeedbreakpoints= 0 100 100 # (m/s), e.g. 0.0 100.0 Rhoair = 1.25 # Air density (kg/m3) PavBnd = 101325 # Average air Pressure on open boundaries, (N/m2), only applied if value > 0 PavIni = 0 Gapres = 101325.0000000 # Global Atmospheric Pressure, (N/m2) [Zsediment] Sedimentmodelnr = 2 # Sediment model nr, (0=no, 1=Krone, 2=SvR2007) SedFile = # Sediment characteristics file (*.sed) MorFile = # Morphology settings file (*.mor) Nr_of_sedfractions = 1 # Nr of sediment fractions, (specify the next parameters for each fraction) MxgrKrone = 0 # Highest fraction index treated by Krone D50 = .0005 # Mean Sandgrain diameter (m), e.g. 0.0001 Rhosed = 2650. # Mean Sandgrain rho (kg/m3) , e.g. 2650 InitialSedimentConcentration = 0. # Initial sediment concentration (kg /m3) Uniformerodablethickness = 20. # Uniform erodable layer thickness (m) Numintverticaleinstein = 7 # Number of vertical intervals in Einstein integrals ( ) Jaceneqtr = 2 # 1=equilibriumtransport at cell centre, 2= at netnode (default) ( ) Morfac = 100. # Morphological acceleration factor (), bottom updates active for morfac > 0, 1d0=realtime, etc TMorfspinup = 0. # Spin up time for morphological adaptations (s) [waves] Wavemodelnr = 0 WaveNikuradse = 0.01 Rouwav = FR84 Gammax = 1 [time] RefDate = 20110801 # Reference date (yyyymmdd) Tzone = 0 DtUser = 6 # User timestep in seconds (interval for external forcing update & his/map output) DtNodal = DtMax = 6 # Max timestep in seconds DtInit = 3 # Initial timestep in seconds Tunit = M # Time units in MDU (H, M or S) TStart = 479520 # Start time w.r.t. RefDate (in TUnit) TStop = 485280 # Stop time w.r.t. RefDate (in TUnit) AutoTimestep = 1 # Use CFL timestep limit or not (1/0) [restart] RestartFile = # Restart file, only map, hence: *_map.nc RestartDateTime = 20110801 # Restart time (yyyymmddhhmmss) [external forcing] ExtForceFile = zm_tide.ext # *.ext ExtForceFileNew = [trachytopes] TrtRou = N TrtDef = TrtL = TrtDt = 60 [output] OutputDir = # 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 = zm_tide_obs.xyn # *.xyn Coords+name of observation stations. CrsFile = zm_tide_crs.pli # *_crs.pli Polyline(s) definining cross section(s). HisFile = zm_tide_his.nc # *_his.nc History file in NetCDF format. HisInterval = 300 # History output, given as "interval" "start period" "end period" (s) XLSInterval = # Interval (s) between XLS history MapFile = zm_tide_map.nc # *_map.nc Map file in NetCDF format. MapInterval = 900 # Map file output, given as "interval" "start period" "end period" (s) RstInterval = 86400 # Interval (s) between map file outputs S1incinterval = # Interval (m) in incremental file for waterlevels S1 WaqInterval = # Interval (in s) between Delwaq file outputs StatsInterval = 3600 # Interval (in s) between simulation statistics output. Writebalancefile = 0 TimingsInterval = # Timings output interval MapFormat = 1 # 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= 1 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 WaqFileBase = # Basename (without extension) for all Delwaq files to be written. SnapshotDir = # Directory where snapshots/screendumps are saved. 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.