Skip to contents

Compute all hydraulic data using new Morel equations

Source: R/TnetHydraulic.R
TnetHydraulic_MorelNew.Rd

Calculate rivers width (\(Bm\)), depth (\(H\)) and water travel time (\(TPS\)) with Discharge as input data, using Morel formulas with simplified parameters. TnetHydraulic_MorelNewRec() can be use to recalibrate some river segments.

Usage

TnetHydraulic_MorelNew(
  path_data,
  shapefile,
  export_files = c("Bm", "H", "TPS")
)

Arguments

path_data

Path to the folder containing Q created by TNET_computeQ()

shapefile

Path to the shapefile with all info on segments. Columns needed are detailed in the Shapefile columns section.

export_files

vector of files that will be exported. It can contain "Bm", "H", "CV" and "TPS", but only "Bm", "H" and "TPS" are needed for T-NET

Value

All hydraulic data (H, Bm, TPS) as NetCDF files in the folder path_data

Equations used

River segment width (\(Bm\))

$$Bm = ad \cdot (Qmean)^{bd} \cdot \left(\frac{Qaval}{Qmean}\right)^{b}$$

with \(ad = 7.482\), \(bd = 0.477\), \(b = 0.148\)

River segment depth (\(H\))

$$H = cd \cdot (Qmean)^{fd} \cdot \left(\frac{Qaval}{Qmean}\right)^{f}$$

with \(cd = 0.340\), \(fd = 0.259\), \(f = 0.292\)

Water speed in the river segment (\(CV\))

$$CV = \frac{Qaval}{H \cdot Bm}$$

Water travel time in river segment (\(TPS\))

$$TPS = \frac{Longueur\_m}{CV} \div 3600$$

where:

  • \(Longueur\_m\) is the length of the river segment in meters.

  • \(Qaval\) is the flow rate at the segment's outlet.

  • \(Qmean\) is the mean flow rate on the entire period.



Shapefiles columns

\(gid\_new\)ID of the Topage segment
\(Longueur\_m\)Length of the river segment (in meter)

References

Maxime Morel, Doug J. Booker, Frédéric Gob, Nicolas Lamouroux, Intercontinental predictions of river hydraulic geometry from catchment physical characteristics, Journal of Hydrology, https://doi.org/10.1016/j.jhydrol.2019.124292.

See also

TnetHydraulic_MorelNewRec()

Examples


TnetHydraulic_MorelNew(path_data = path/to/data,
                      shapefile = path/to/shapefile.shp)
#> Error: object 'to' not found

############################################################
## Using this function with TNET_initializeSim() function ##
############################################################
infoSimu <- TNET_initializeSim(...)
#> Error: '...' used in an incorrect context

#Do discharge computation needed using a Tnet*_computeQ() function

TnetHydraulic_MorelNew(path_data = infoSimu$hydroPath, 
                      shapefile = infoSimu$TOPAGE_shape)
#> Error: object 'infoSimu' not found