View source for TUFLOW HR Output

<font size="18">Page Under Construction</font> 


=Introduction=
High Resolution (HR) Grid/Raster Map Outputs was first introduced in TUFLOW 2020-10-AB release. When an SGS model uses “<font color="blue"><tt>SGS Approach </tt></font> <font color="red"><tt>== </tt></font> <font color="black"><tt>Method C</tt></font>”, the sampled elevations are retained at the end of the geometry processing. These sub-grid elevations include topography modifiers such as breaklines, and they allow a high-resolution elevation check file to be written and used for high resolution depth map outputs. Currently, ASC, FLT and TIF raster formats are supported.

For the HR output the water level at each HR output location is interpolated from the computed 2D water levels. The depth is the difference between the interpolated water level and the sub-grid elevation. This differs from the standard depth output which calculates the depths at cell centres and corners first, then interpolate the depths to the standard output grid locations(default is half the cell size). The difference in the interpolation process is illustrated below.

[[File:Sgs_std_output.png|360px]] [[File:Sgs_hr_output.png|360px]]

The advantage of the HR output is that it can retain the sub-grid detail of the terrain information even at coarse cell size. As compared in the example below, the HR depth output shows clear flow path even at 100m grid. However, nicer depth output does not mean reliable hydraulic result. In fact, the 100m cell size is too coarse to produce reliable/converged hydraulic results for any real-world flood model. It's is strongly recommended that the model convergence/benchmarking tests must be conducted based on the standard map output, and the HR output should be used for presentation purpose only.

[[File:Sgs_std_vs_hr_output.png|720px]]

=Enable High Resolution Output=
To produce high-resolution output, add “HRASC”, “HRFLT” or “HRTIF” to the “<font color="blue"><tt>Map Output Format </tt></font> <font color="red"><tt>==</tt></font>” command, for example: 

<font color="blue"><tt>Map Output Format </tt></font> <font color="red"><tt>== </tt></font> <font color="black"><tt>XMDF TIF HRTIF</tt></font>

At the moment, only water level and depth are supported as the High Resolution Output formats, i.e.:

<font color="blue"><tt>HRTIF Map Output Data Types </tt></font> <font color="red"><tt>== </tt></font> <font color="black"><tt>h d</tt></font>

The map output data types and output interval can be defined separately for the High Resolution Output formats. For example: 

<font color="blue"><tt>Map Output Interval </tt></font> <font color="red"><tt>== </tt></font> <font color="black"><tt>600</tt></font> <font color="green"><tt>! Sets all map output interval to 600 seconds</tt></font>

<font color="blue"><tt>HRTIF Map Output Interval </tt></font> <font color="red"><tt>== </tt></font> <font color="black"><tt>3600</tt></font> <font color="green"><tt>! Sets HR map output interval to 3600 seconds</tt></font>
 
The output resolution for high resolution grid output is defined using the following .tcf command:

<font color="blue"><tt>HR Grid Output Cell Size </tt></font> <font color="red"><tt>== </tt></font> <font color="black"><tt>0.5</tt></font> <font color="green"><tt>! m or ft</tt></font>

If this command is omitted, the default output resolution is set as the SGS sampling distance.


HR outputs are supported by HPC and Quadtree Output Zones. Below is an example of tcf commands to set up HR output in a output zone:
 <font color="blue"><tt>Model Output Zones</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>ZoneA</tt></font><br>
 <font color="blue"><tt>Define Output Zone</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>ZoneA</tt></font>
    <font color="blue"><tt>Read GIS Output Zone</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>..\model\gis\2d_oz_S05_002_R.shp</tt></font>
    <font color="blue"><tt>Write Check Files INCLUDE</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>dem_Z</tt></font> <font color="green"><tt>! includes only dem_Z check file for this output zone</tt></font>
    <font color="blue"><tt>Map Output Format</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>HRTIF</tt></font>
    <font color="blue"><tt>HRFLT Map Output Data Types</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>h d</tt></font>
    <font color="blue"><tt>HRFLT Map Output Interval</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>3600</tt></font>
    <font color="blue"><tt>HR Grid Output Cell Size</tt></font> <font color="red"><tt>==</tt></font> <font color="black"><tt>0.5</tt></font>
 <font color="blue"><tt>End Define</tt></font>

If HR output is defined in the Map Output Format, a high-resolution DEM_Z check file (DEM_Z_HR) is produced instead of the DEM_Zmin check file.

=Water Level Interpolation Method=
Similar to the standard output, the HR output also need to interpolate cell centre water level to cell corners. However, the interpolation methods for the standard output (<font color="blue"><tt>Map Output Corner Interpolation </tt></font> <font color="red"><tt>== </tt></font> <font color="black"><tt>Method C</tt></font>) can produce "bumpy" HR water level output in direct rainfall model along narrow streams. This often happens between fully wet cells and sheet flow cells, as the water level is linearly interpolated 
[[]]
applies a minimum corner water level of “corner elevation + cell wet/dry depth” if all the surrounding 4 cells are wet. However, this may produce “bumpy” water level outputs along narrow streams in a direct rainfall SGS model between sheet flow cells and fully wet cells when using the new high resolution output. This is caused by applying a linear water level interpolation between SGS cells that have a non-linear minimum SGS elevation as illustrated below.

interpolated from the surrounding SGS sampled elevations.
*''under construction''

=Interpolation near Thin Breakline  =
*''under construction''
= Interpolation for Sheet Flow=
*''under construction''



Any further questions please email TUFLOW support: [mailto:support@tuflow.com?Subject=TUFLOW%201D%20gates%20help support@tuflow.com] 


{{Tips Navigation
|uplink=[[ TUFLOW 1D Channels and Hydraulic Structures | Back to 1D Channels and Hydraulic Structures]]
}}