User manual - forcing¶
Overview¶
The input for SFINCS is supplied using various text and binary files, which are linked through the main input file: sfincs.inp. Within this section of the user manual all different types of forcing settings and files are discussed. The figure below gives an overview of all different types of input files and whether they are required or not. Below an example is given of this file, which uses a keyword/value layout. For more information regarding specific parameters see the pages ‘Input parameters’ or ‘Output parameters’.
NOTE - In the manual below, blocks named ‘Matlab example using OET’ are included, referring to easy setup scripts included in the SFINCS’ Open Earth Tools Matlab set of scripts: https://svn.oss.deltares.nl/repos/openearthtools/trunk/matlab/applications/sfincs
Fig. 13 Overview of input file of SFINCS with indication whther they are required or not¶
Forcing¶
SFINCS has different functionalities regarding different relevant physical processes for compound flooding and what type of model is required. At first nearshore/offshore water levels can be specified at the different locations along the coast to include tides and storm surge levels. Inland drivers of flooding like precipitation and wind can be specified in a number of ways. This varies from simple spatially uniform time-series to spatially varying spiderwebs or grid input types. Furthermore, simple implementations for discharges are included.
Discussed in this user manual are the water-level boundaries, discharge points, wind, pressure, rain and waves.
Water levels¶
To specify water-level time-series to the boundary cells (msk=2), first the input locations have to be specified in ‘sfincs.bnd’. For every boundary point there is interpolated with a weighted average between the two closest input locations.
Fig. 14 Example of the weighted interpolation of boundary input locations to the msk=2 boundary cells.¶
Water level points¶
bndfile = sfincs.bnd
<bnd1 x1> <bnd1 y1>
<bnd2 x2> <bnd2 y2>
e.g.
400000 1200000
480000 1250000
Matlab example using OET
points.x(1) = 400000;
points.y(1) = 1200000;
points.x(2) = 480000;
points.y(2) = 1250000;
sfincs_write_boundary_points(inp.bndfile,points)
Water level time-series¶
Then in the file ‘sfincs.bzs’ the water level time-series are specified per input location. Times are specified in seconds with respect to SFINCS’ internal reference time ‘tref’, as specified in sfincs.inp.
NOTE - The time and length of values you prescibe per boundary input location should be all the same in the bzsfile
bzsfile = sfincs.bzs
<time 1> <zs1 bnd1> <zs1 bnd2>
<time 2> <zs2 bnd1> <zs2 bnd2>
e.g.
0 0.50 0.75
3600 0.60 0.80
7200 0.45 0.85
Matlab example using OET
time = [0, 3600, 7200];
waterlevels = [0.5, 0.75; 0.6, 0.8; 0.45, 0.85];
sfincs_write_boundary_conditions(inp.bzsfile,time,waterlevels)
Waves¶
When forcing waves, besides providing a bzsfile with slowly varying water level time-series, also the same type of file with the quickly varying water level component due to waves can be prescribed. This can contain infragravity and/or short waves. Do note that the forced signal should be the incoming wave component only, not including the reflecting one, since this is computed by SFINCS internally as well. The signal should be around 0.
NOTE - Specified time should be the same in both the bzs and bzi files
bzifile = sfincs.bzi
<time 1> <zi1 bnd1> <zi1 bnd2>
<time 2> <zi2 bnd1> <zi2 bnd2>
e.g.
0 0.05 0.07
2 -0.02 -0.04
4 0.10 0.03
Netcdf format input¶
As alternative, the bnd/bzs/bzi data can also be specified using a single Netcdf file with FEWS input type format ‘netbndbzsbzifile’. Making this format netcdf file can be easily done using the OET Matlab script ‘sfincs_write_netcdf_bndbzsbzifile.m’
Matlab example using OET
inp.netbndbzsbzifile = 'sfincs_netbndbzsbzifile.nc';
x = [0, 100, 200];
y = [50, 150, 250];
EPSGcode = 32631;
UTMname = 'UTM31N';
refdate = '1970-01-01 00:00:00';
% possibly use formatOut = 'yyyy-mm-dd HH:MM:SS'; datestr(tref, formatOut);
time = [0, 60];
rng('default');
bzs = -1 * randi([0 10],length(time),length(x));
bzi = -1 * randi([0 10],length(time),length(x));
sfincs_write_netcdf_bndbzsbzifile(inp.netbndbzsbzifile, x, y, EPSGcode, UTMname, refdate, time, bzs, bzi)
Discharges¶
A simple implementation of discharge points is added to SFINCS, specify values in m^3/s. First specify the locations in ‘sfincs.src’ and then the discharge time-series in ‘sfincs.dis’. Alternatively, you can provide this as netcdf file in the Delft-FEWS format.
Fig. 15 Example of how discharge input points from 2 different sources are snapped to the grid of SFINCS.¶
Discharge points¶
srcfile = sfincs.src
<src1 x1> <src1 y1>
<src2 x2> <src2 y2>
e.g.
300000 1500000
380000 1650000
Matlab example using OET
points.x(1) = 300000;
points.y(1) = 1500000;
points.x(2) = 380000;
points.y(2) = 1650000;
sfincs_write_boundary_points(inp.srcfile,points)
Discharge time-series¶
Then in the file ‘sfincs.dis’ the discharge time-series are specified per input location. Times are specified in seconds with respect to SFINCS’ internal reference time ‘tref’, as specified in sfincs.inp.
disfile = sfincs.dis
<time 1> <dis1 src1> <dis1 src2>
<time 2> <dis2 src1> <dis2 src2>
e.g.
0 100 1000
3600 300 1100
7200 0 1300
Matlab example using OET
time = [0, 3600, 7200];
discharge = [100, 1000; 300, 1100; 0, 1300];
sfincs_write_boundary_conditions(inp.disfile,time,discharge)
Netcdf format input¶
As alternative, the src/dis data can also be specified using a single Netcdf file with FEWS input type format ‘netsrcdisfile’. SFINCS assumes that the input variables ‘x’, ‘y’, ‘time’, ‘discharge’ and ‘stations’ are available in the netcdf file, including a reference time as UNIT in variable ‘time’ of the Fews time format: “minutes since 1970-01-01 00:00:00.0 +0000”
Matlab example using OET
inp.netsrcdisfile = 'sfincs_netsrcdisfile.nc';
x = [0, 100, 200];
y = [50, 150, 250];
EPSGcode = 32631;
UTMname = 'UTM31N';
refdate = '1970-01-01 00:00:00';
% possibly use formatOut = 'yyyy-mm-dd HH:MM:SS'; datestr(tref, formatOut);
time = [0, 60];
rng('default');
dis = -1 * randi([0 10],length(time),length(x));
sfincs_write_netcdf_srcdisfile(inp.netsrcdisfile, x, y, EPSGcode, UTMname, refdate, time, dis)
Meteo¶
There are a few different options to specify wind and rain input:
Use a spatially varying spiderweb input (as in Delft3D/Delft3D FM) for forcing tropical cyclones only the wind and pressure input, or for the wind as well as the rain input.
Use a spatially varying grid input (as in Delft3D) for u- and v- wind velocities and/or the rain and/or pressure input.
Use a spatially varying grid input using a netcdf file based on a FEWS input type format for wind, rain and/or atmospheric pressure input.
Use a spatially uniform input for wind and rain, which is faster but also more simplified.
Make a combination, for instance use a spiderweb for the wind input and a spatially uniform rain-input. When combining, test whether the forcing is as wanted since not all combinations of the above options might be possible and/or changing depending on specific code version.
Fig. 16 Overview of possible meteo input file options and names¶
Spatially varying spiderweb¶
The option of forcing spiderweb files is only relevant for tropical cyclones, best is to put grid units in the same projected coordinate reference system (UTM zone) as the SFINCS grid. For generation of these spiderweb files use Deltares’ Wind Enhancement Scheme tool (WES, see https://content.oss.deltares.nl/delft3d/manuals/Delft3D-WES_User_Manual.pdf or OET Matlab equivalent) or get in touch.
Spiderweb-input:
spwfile = tropical_cyclone.spw
Spatially varying gridded¶
Spatially varying meteo input on constant grid can be forced using the native Delft3D type meteo input files, using the same file conventions. For wind this is wind in x-&y-direction (amu, amv), precipitation (ampr) and atmospheric pressure (amp). The grid has a constant resolution dx&dy, which can be in the native (usually coarser than your SFINCS grid) resolution of the meteo data. Within SFINCS this is interpolated onto the actual SFINCS grid.
Wind:
**amufile = sfincs.amu**
within amufile:
quantity1 = x_wind
unit1 = m s-1
**amvfile = sfincs.amv**
within amvfile:
quantity1 = y_wind
unit1 = m s-1
Rain:
**amprfile = sfincs.ampr**
within amprfile:
quantity1 = precipitation
unit1 = mm/hr
Atmospheric pressure:
**ampfile = sfincs.amp**
within ampfile:
quantity1 = air_pressure
unit1 = Pa
Delft3D-meteo ascii type input:
These files have this general header of 13 lines which SFINCS expects (Check this after creating your input files!), after which the TIME and data blocks are given per time frame. Only use 1 quantity per file:
FileVersion = 1.03
filetype = meteo_on_equidistant_grid
n_cols = 2
n_rows = 4
grid_unit = m
x_llcorner = 417328
y_llcorner = 3495537
dx = 5000
dy = 5000
n_quantity = 1
quantity1 = x_wind
unit1 = m s-1
NODATA_value = -999
TIME = 90831.0 hours since 1970-01-01 00:00:00 +00:00 # 1980-05-12 15:00:00
0 0 0 0
0 0 0 0
TIME = 90831.0 hours since 1970-01-01 00:00:00 +00:00 # 1980-05-12 15:00:00
0 0 0 0
0 0 0 0
Matlab example using OET
data.parameter.time = datenum(2018,01,01):3/24:datenum(2018,01,02);
data.parameter.x = 0:5000:25000;
data.parameter.y = 10000:5000:40000;
data.parameter.val = ones(length(data.parameter.time), length(data.parameter.y), length(data.parameter.x));
write_meteo_file_delft3d(inp.amufile, data, 'x_wind', 'm s-1', datenum(1970,01,01), varargin);
see 'write_meteo_file_delft3d.m' for more information.
Spatially varying gridded netcdf¶
The same spatially varying gridded input as using Delft3d’ ascii input files can be specified using FEWS compatible Netcdf input files. Here for the wind the amu&amv files are combined into 1 Netcdf file (netamuamvfile), the precipitation is in a separate input file (netamprfile) as well as the atmospheric pressure (netampfile).
Making this format netcdf file can be easily done using the OET Matlab scripts ‘sfincs_write_netcdf_amuamvfile.m’, ‘sfincs_write_netcdf_amprfile.m’and ‘sfincs_write_netcdf_ampfile.m’. See those files for more information.
Matlab example using OET - netamuamvfile
inp.netamuamvfile = 'sfincs_netamuamvfile.nc';
x = [0, 100, 200];
y = [50, 150, 250];
EPSGcode = 32631;
UTMname = 'UTM31N';
refdate = '1970-01-01 00:00:00';
% possibly use formatOut = 'yyyy-mm-dd HH:MM:SS'; datestr(tref, formatOut);
time = [0, 60];
rng('default');
amu = -1 * randi([0 10],length(time),length(y),length(x));
amv = 1 * randi([0 10],length(time),length(y),length(x));
sfincs_write_netcdf_amuamvfile(inp.netamuamvfile, x, y, EPSGcode, UTMname, refdate, time, amu, amv)
Matlab example using OET - netamprfile
inp.netamprfile = 'sfincs_netamprfiles.nc';
x = [0, 100, 200];
y = [50, 150, 250];
EPSGcode = 32631;
UTMname = 'UTM31N';
refdate = '1970-01-01 00:00:00';
% possibly use formatOut = 'yyyy-mm-dd HH:MM:SS'; datestr(tref, formatOut);
time = [0, 60];
rng('default');
ampr = -1 * randi([0 10],length(time),length(y),length(x));
sfincs_write_netcdf_amprfile(inp.netamprfile, x, y, EPSGcode, UTMname, refdate, time, ampr)
Matlab example using OET - netampfile
inp.netampfile = 'sfincs_netampfiles.nc';
x = [0, 100, 200];
y = [50, 150, 250];
EPSGcode = 32631;
UTMname = 'UTM31N';
refdate = '1970-01-01 00:00:00';
% possibly use formatOut = 'yyyy-mm-dd HH:MM:SS'; datestr(tref, formatOut);
time = [0, 60];
rng('default');
amp = -1 * randi([0 10],length(time),length(y),length(x));
sfincs_write_netcdf_ampfile(inp.netampfile, x, y, EPSGcode, UTMname, refdate, time, amp)
Spatially uniform¶
Spatially uniform wind:
‘vmag’ is the wind speed in m/s, ‘vdir’ is the wind direction in nautical from where the wind is coming. The file can be make using OET Matlab script ‘sfincs_write_boundary_conditions.m’. Times are specified in seconds with respect to SFINCS’ internal reference time ‘tref’, as specified in sfincs.inp.
wndfile = sfincs.wnd
<time 1> <vmag1> <vdir1>
<time 2> <vmag2> <vdir2>
e.g.
0 5 120
3600 15 180
7200 10 165
Spatially uniform rain:
Rain input in mm/hr, times are specified in seconds with respect to SFINCS’ internal reference time ‘tref’, as specified in sfincs.inp. The file can be make using OET Matlab script ‘sfincs_write_boundary_conditions.m’.
precipfile = sfincs.prcp
<time 1> <prcp0>
<time 2> <prcp1>
e.g.
0 0
3600 15
7200 10