Tutorial Module04 Draft

From TUFLOW FV Wiki
Jump to: navigation, search

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.

Tute4 026.jpg

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).

NETCDF 000.png