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.0e3, 5.0e2  
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 thirdgeneration 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 TUFLOWFV 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’ & LAYOUT 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).