make.entrainment

make.entrainment generates files with a uniform thickness of entrainable material. The typical usage for make.entrainment differs from make.scoops and make.slabs in that the user specifies the extent of the area not covered by material that can be entrained rather than the area within the landslide extent. This method was developed so the same file indicating the initial landslide extent could be used for both make.scoops or make.slabs and make.entrainment.

  1. Prepare a topography file as a geotif. If the area of interest includes submarine regions, the topography and bathymetry should be represented as a seamless surface in a single file.

  2. Prepare or obtain a shapefile containing one or more polygons indicating landslide extents. Ensure that the geotif and shapefile have the same coordinate system.

  3. Run make.entrainment to generate input files for D-Claw. Inspect the output and adjust any input parameters as needed.

  4. Using these files, set up a D-Claw simulation to run and analyze. digger will not do this for you.

A code snippet that uses make.entrainment may be found in the file digger/examples/pre-run/synthetic/make_entrainment_synthetic_example.py,

"""An example of using digger.make.entrainment.

This example uses synthetic topography and provides an example of using
digger.make.entrainment.

It will generate tifs and topotype3 files with 5 m of material
everywhere except within the polygon "ScoopsPoly.shp".
"""

from digger import make

make.entrainment(
    topo_path="../../../data/topo_tin5_raster.tif",
    shapefile_path="../../../data/Scoop_Poly.shp",
    depth=5.0,
    write_tt3=True,
    write_tif=True,
    fig_path="entrainment_source.png",
)

Attention

This code snippet is not fully self-sufficient. To reproduce the example, execute the code such that the relative path is correct. This may be done by executing the example script from within the directory in which it is located.

After this code runs, it will produce geotif and topotype3 files that specify the spatial extent and initial depth of erodible material (aux5).

It will also create diagnostic figures to assist the user.

Example digger.make.entrainment summary figure.

Fig. 4 An example of the diagnostic output provided by digger.make.entrainment.