# 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 | |

Spatial Order == 1,2 | |

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.0001,100 | |

Display dt == 900 | |

! 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. An external GOTM model is used to model turbulence. The timestep for updating the vertical eddy-viscosity and scalar-diffusivity terms is set to 1800 seconds. The limits for these terms are also specified. |

stability limits == 10.0,10.0 | |

Momentum mixing model == Smagorinsky | |

Global horizontal eddy viscosity == 0.2 | |

vertical mixing model == External | |

External turbulence model directory == ..\GOTM\ | |

turbulence update dt == 1800 | |

global vertical eddy viscosity limits == 1.0e-4, 1.0 | |

global vertical scalar diffusivity limits == 0., 1.0 | |

! 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 == ..\geo\hydraul_006.2dm | |

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

echo geometry netcdf == 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 salinity == 10. | |

Initial Water Level == 0.052 | |

! 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\Upstream_Q_Temp_Sal_Sed_001.csv | |

bc header == time_hr,flow_m3s-1,sal_ppt,temp_degC | |

end bc | |

bc == WL, 2, ..\bc\Downstream_H_Temp_Sal_Sed_001.csv | |

bc header == time_hr,tidelevel_m,sal_ppt,temp_degC | |

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, sal and Rhow mean that output files containing water level, depth, velocity, temperature, salinity and density will be created. Output interval is a command used to specify the desired output interval in seconds. |

Logdir == log | |

output dir == ..\output\ | |

output == datv | |

output parameters == h,v,d, temp, sal | |

output interval == 3600. | |

end output | |

output == netcdf | |

output parameters == h,v,d, temp, sal, Rhow | |

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.