# HYD 000

The control file HYD_000.fvc should include the syntax commands listed on the lefthand side of the following table.

! TUFLOW FV TUTORIAL | The first lines are a description of the model simulation. |

! Estuary Model – 2D ONLY | |

GIS Format == SHP | These lines contain the commands to define the GIS format and model projection. For more information on configuring a model this way please refer to Configure model using QGIS . |

SHP Projection == ..\model\gis\projection.prj | |

Tutorial Model == ON | Enables licence free modelling. |

! GENERAL PARAMETERS | The ks friction approach is the preferred bottom drag model for 3D simulations. Temperature and salinity are added here for example purposes, in 2D these essentially act as tracers. Spatial order (horizontal, vertical) specifies the order of accuracy of the solution schemes used in the simulation. Second order spatial accuracy will typically be required in the vertical direction when trying to resolve stratification. Atmospheric heat calculations are not to be included. |

Bottom Drag Model == ks | |

Include Salinity == 1,0 | |

Include Temperature == 1,0 | |

Include Heat == 0 | |

! TIME COMMANDS | The time commands include the time format (ISO = date and time), start and end times. The CFL limit is 0.60 (less than the default – TUFLOW FV then assigns a timestep at each computational step according to the CFL limit and within the ranges specified in the timestep limits. Display dt sets the simulation time interval (in seconds) for displaying timestep information. |

Time Format == ISODATE | |

Start Time == 01/05/2011 00:00:00 | |

End Time == 07/05/2011 00:00:00 | |

CFL == 0.6 | |

Timestep Limits == 0.1,30. | |

Display dt == 900 | |

Cell Dry/Wet Depths == 0.001, 0.02 | |

! MODEL PARAMETERS | The stability limits command specifies a maximum water level and maximum velocity that indicate an unstable model. The simulation will stop if these limits are exceeded. A Smagorinsky eddy viscosity approach has been specified, with a Smagorinsky factor of 0.2. |

Stability Limits == 10.0,10.0 | |

! Horizontal Turbulent Mixing | |

Momentum Mixing Model == Smagorinsky | |

Global Horizontal Eddy Viscosity == 0.2 | |

Global Horizontal Eddy Viscosity Limits == 0.05, 99999. | |

! Horizontal Scalar Mixing | |

Scalar Mixing Model == Smagorinsky | |

Global Horizontal Scalar Diffusivity == 0.2 | |

Global Horizontal Scalar Diffusivity Limits == 0.05, 99999. | |

! GEOMETRY | A spherical coordinate system is specified, whereby geometry inputs and computational coordinates are in degrees.
The model geometry is the 2dm mesh provided for this tutorial module. This 2dm file includes mesh/element configuration and topography/bathymetry data. External boundary condition locations are included in the SHP file that is read in by the Read GIS Nodestring command. For more information regarding GIS nodestring definition please refer to Configure model using QGIS . The echo function ensures the model writes a *_geo.nc (netCDF format geometry) check file. This will be required for curtain plots. |

Spherical == 1 | |

Geometry 2D == ..\model\geo\hydraul_006.2dm | |

Read GIS Nodestring == ..\model\gis\2d_ns_OpenBoundaries_001_L.shp | |

Echo Geometry == 1 | |

! MATERIAL PROPERTIES | The Set Mat command will set the cell material ID to 1 globally. The second command, Read GIS Mat will read in an SHP file to override the cell material ID for the polygons defined with materials 2 and 3. In this way, three material types have been specified, each with a corresponding ks value. These include:
1 = Main river area (Default) |

Set Mat == 1 | |

Read GIS Mat == ..\model\gis\2d_mat_Materials_002_R.shp | |

Material == 1 | |

Bottom Roughness == 0.02 | |

End Material | |

Material == 2 | |

Bottom Roughness == 0.75 | |

End Material | |

Material == 3 | |

Bottom Roughness == 1.3 | |

End Material | |

! INITIAL CONDITIONS | The initial internal model boundary conditions. |

Initial Temperature == 20. | |

Initial Water Level == 0.052 | |

Initial Salinity ==10. | |

! BOUNDARY CONDITIONS | The boundary conditions are defined in csv files containing water level and flow information along the external boundaries of the model. The upstream flow (Q) and downstream water level (WL) boundaries are specified and linked to nodestrings 1 and 2 as defined in the 2dm mesh file. |

BC == Q, 1, ..\bc_dbase\Upstream_Q_Temp_Sal_Sed_001.csv | |

BC Header == time_hr,flow_m3s-1,sal_ppt,temp_degC | |

End BC | |

BC == WL, 2, ..\bc_dbase\Downstream_H_Temp_Sal_Sed_001.csv | |

End BC | |

! OUTPUT COMMANDS | In this instance, datv and netcdf format files are specified. While the datv format is easily read into SMS for viewing, the netcdf is ideal for processing in matlab. The h, d, v, temp, and sal mean that output files containing water level, depth, velocity, temperature, and salinity will be created. Output interval is a command used to specify the desired output interval in seconds. |

Logdir == log | |

Output Dir == ..\results\ | |

Output == netcdf | |

Output Parameters == h,v,d, temp, sal | |

Output Interval == 3600. | |

End Output | |

Output == flux | |

Output Interval == 900. | |

End Output | |

Output == mass | |

Output Interval == 900. | |

End Output |

To return to Tutorial Module 5 please refer to Tutorial Module 5.