Tutorial Module04 Draft
Cyclone/Hurricane Hydrodynamic Model (Tide + External Wind/Pressure and Wave Model)
As an update to the previous Holland parametric wind model example, this example will use NETCDF format files as wave and wind boundary condition inputs. These boundary condition inputs have been derived from data which was collected from Cyclone Paul, which crossed the Gulf in 2010.
The syntax additions/modifications to the holland parametric wind/pressure model are highlighted below. After making the syntax updates, save the TUFLOW FV control file with the name, 'GCarp_Holland_002.fvc'. The control file should be saved to the input folder within the model directory.
! TUFLOW FV TUTORIAL | Model description update |
! Gulf of Carpentaria – NETCDF Boundary Condition Inputs
| |
! GENERAL PARAMETERS | No syntax changes |
bottom drag model == ks
| |
! TIME COMMANDS | updated simulation end time:
end time == 03/04/2010 23:00:00 |
time format == ISODATE | |
start time == 26/03/2010 16:30:00 | |
end time == 03/04/2010 23:00:00 | |
cfl == 0.95 | |
timestep limits == 0. 1, 30.
| |
! MODEL PARAMETERS | No syntax changes |
cell wet/dry depths == 5.0e-3, 5.0e-2 | |
stability limits == 50. ,10. | |
scalar mixing model == Elder | |
momentum mixing model == Smagorinsky | |
global horizontal eddy viscosity == 0.5
| |
!GEOMETRY | No syntax changes |
geometry 2d == ..\geo\Gulf_of_Carpentaria_001.2dm | |
cell elevation file == ..\geo\Gulf_of_Carpentaria_001_2dm.csv, coordinate | |
spherical == 1
| |
! MATERIAL PROPERTIES | No syntax changes |
material == 1 | |
- bottom roughness == 0.05 | |
-end material
| |
! INITIAL CONDITIONS | No syntax changes |
initial waterlevel ==0
| |
! BOUNDARY CONDITIONS | The parametric Holland wind/pressure model input from the previous example (bc == CYC_HOLLAND) have been replaced with NETCDF inputs.
grid definition file == ..\bc\BOM_HOLLAND_WIND.nc grid definition variables == lon,lat end grid bc == W10_GRID, 1, ..\bc\BOM_HOLLAND_WIND.nc bc header == time,u,v bc update dt == 450.0 end bc Wave boundary condition inputs have also been added, from which TUFOW FV calculates the hydrodynamics. grid definition file == ..\bc\WAVE.nc grid definition variables == longitude,latitude end grid bc == Wave, 2, ..\bc\WAVE.nc bc header == time,hs,tps,thetap bc reference time == 01/01/1970 00:00 bc time units == seconds bc update dt == 3600. end bc
|
bc == WLS, 1, ..\bc\TIDE.csv | |
bc header == TIME,NS1,NS2 | |
end bc | |
bc == WLS, 2, ..\bc\TIDE.csv | |
bc header == TIME,NS2,NS3 | |
end bc | |
bc == WLS, 3, ..\bc\TIDE.csv | |
bc header == TIME,NS3,NS4 | |
end bc | |
bc == WLS, 4, ..\bc\TIDE.csv | |
bc header == TIME,NS5,NS6 | |
end bc | |
bc == WLS, 5, ..\bc\TIDE.csv | |
bc header == TIME,NS6,NS7 | |
end bc | |
grid definition file == ..\bc\BOM_HOLLAND_WIND.nc | |
grid definition variables == lon,lat | |
end grid | |
bc == W10_GRID, 1, ..\bc\BOM_HOLLAND_WIND.nc | |
bc header == time,u,v | |
bc update dt == 450.0 | |
end bc | |
grid definition file == ..\bc\WAVE.nc | |
grid definition variables == longitude,latitude | |
end grid | |
bc == Wave, 2, ..\bc\WAVE.nc | |
bc header == time,hs,tps,thetap | |
bc reference time == 01/01/1970 00:00 | |
bc time units == seconds | |
bc update dt == 3600. | |
end bc
| |
! OUTPUT COMMANDS | Additional wave height, period and direction outputs have been added.
output parameters == h,v,W10,wvht,wvper,wvdir In addition to datv output, NECDF output has also been specified. output == netcdf output parameters == H,V,W10,wvht,wvper,wvdir output interval == 900. end output
|
output dir == ..\Output\ | |
output == datv | |
output parameters == h,v,W10,wvht,wvper,wvdir | |
output interval == 900 | |
end output | |
output == netcdf | |
output parameters == H,V,W10,wvht,wvper,wvdir | |
output interval == 900. | |
end output |
NETCDF Boundary Condition Inputs
As an update to the previous Holland parametric wind model example, this example will use NETCDF format files as wave and wind boundary condition inputs. These boundary condition inputs have been derived from data which was collected from Cyclone Paul, which crossed the Gulf in 2010.
SWAN is a third-generation wave model developed at the Delft University of Technology (Delft University of Technology, 2006). SWAN version 40.91AB can be configured to output NETCDF format results which can be used as wave inputs to TUFLOW-FV models. The outputs can be produced for the SWAN computational grid or for a defined output grid using FRAME or GROUP commands (SWAN User Manual). The format of the SWAN output command is:
BLOCK ‘grid_name’ NOHEADER ‘out_file_name’ & LAY-OUT 3 HSIGN TPS PDIR OUTPUT yyyymmdd.HHMMSS dt dtunit
For the purposes of this tutorial, SWAN model wave forcing output from a simulation of Cyclone Paul (which crossed the Gulf in 2010) has been provided.
As an update to the previous Holland wind model example, in addition to the wave boundary condition input, wind has be applied as a boundary condition using NETCDF format. The NETCDF file specifies 10m wind speed time series on an array of grid points. This file has also been provided.
Run TUFLOW FV
Once you’re happy with the fvc file contents, run TUFLOW FV.
Useful links:
Tide + External Wind/Pressure and Wave Model Results
'DAT' and 'NETCDF' output types have been specified for this simulation. Review the DAT results either in SMS or one of the following GIS packages; MapInfo, ArcGIS, QGIS and SAGA (see DOS Utilities).
NETCDF results can be read using a range of freeware software. We however recommend using Matlab. A range of matlab functions have been created which automated many common TUFLOW FV result processing tasks (see Matlab Utilities).