TAU0 is an input in the fort.15 file that influences the degree of numerical diffusion in ADCIRC's governing equations. Specifically, it influences the weighting factor that determines the relative contribution of the primitive and wave portions of the the Generalized Wave-Continuity Equation (GWCE). The weighting factor, , is affected by values in both the fort.15 file and the fort.13 file, if the primitive weighting in continuity equation or min and max primitive weighting in continuity equation nodal attributes are specified. This page addresses both the
TAU0 value in the fort.15 file and the parameter more generally.
TAU0 value specified in the fort.15 file can be either a flag (indicating how ADCIRC should operate) or the actual value used in solving the GWCE, it is important to distinguish between the two. All negative
TAU0 are flags, all positive
TAU0 are . The following table is a summary of possible
TAU0 values and their meaning. Note that for
TAU0 = -x.1 where
x is an integer, behavior is the same as
-x, but the values are written to the fort.90 file. More on this below in Outputting.
|Varies in space||no||yes||yes||yes||yes||yes||yes|
|Varies in time||no||no||no||yes||yes||yes||yes|
|Requires nodal attribute||no||no||no||yes||no||yes||no|
For all positive values, the value specified in the fort.15 file is spatially and temporally constant and applied directly, i.e.
TAU0 = , so:
TAU0 = 0, the GWCE is a pure wave equation
TAU0 = 1, the GWCE behaves like a pure primitive continuity equation
0 < TAU0 < 1, the GWCE is a weighted sum of the primitive and wave portions of the equation
Spatially varying, temporally constant
TAU0 = -1, is computed as follows:
If depth >= 10; then = 0.005
If depth < 10; then = 0.020
TAU0 = -2, is is computed according to depth as follows:
If depth >= 200; then = 0.005
If 1 < depth < 200; then = 1/depth
If depth < 1; then
Spatially and temporally varying
TAU0 = -3, is computed from
TAU0Base(read in from the nodal attribute file) as follows:
If TAU0Base < 0.025; then = TAU0Base(constant in time)
If TAU0Base >= 0.025; then = TAU0Base+1.5*TKwhere
TAU0 = -5, is computed similar to
TAU0 = -3as follows:
- It is limited by
Tau0FullDomainMin <= <= Tau0FullDomainMax, which are specified on the following line of the fort.15 file (only when
TAU0 = -5)
- If the min and max primitive weighting in continuity equation nodal attribute is used, then the nodal minimum and maximum values replace the full-domain values in the above calculation.
Spatial and Temporal Updating
TAU0=-3, -5, -6, or -7, is updated in space and in time.
For TAU0 formulations that vary spatially or temporally, ADCIRC can output the internally-calculated nodal values. They are written to the fort.90 file, which has the same format and output frequency as the water surface elevation output file (fort.63). fort.90 output is activated by placing a 1 in the tenths place of TAU0 in the fort.15 file. For example, if TAU0=-3.1, the calculation of is still carried out according to the description of TAU0=-3 above, and the fort.90 output file will also be produced.
For positive values (
TAU0=), a good rule of thumb for setting
TAU0 is to set it equal to the largest value of an equivalent linear friction factor (e.g, for linear friction
TAU0 = TAU; for quadratic friction
TAU0 = maximum (speed*Cd/depth). Typical values for
TAU0 are in the range of 0.005 – 0.1.
TAU0 = -3, paired with the primitive weighting in continuity equation nodal attribute is generally the most popular formulation. In this case,
TAU0Base nodal attribute values can be generated with the ADCIRC utility program tau0_gen.f. The program bases generation on the following logic:
- If the avg. distance between a node and its neighbors < 1750 m TAU0Base = 0.03
- If the avg. distance between a node and its neighbors > 1750 m AND depth < 10m; TAU0Base = 0.02 (TAU0 is constant in time)
- If the avg. distance between a node and its neighbors > 1750 m AND depth > 10m; TAU0Base = 0.005 (TAU0 is constant in time)