NAME
m.swim.glaciers.py - Preprocessor for the Soil and Water Integrated Model - Glacier dynamics (SWIM-G), requires pandas.KEYWORDS
SYNOPSIS
m.swim.glaciers.py
m.swim.glaciers.py --helpm.swim.glaciers.py [-dsak] elevation=name accumulation=name glacierarea=name subbasins=name landuse=name soil=name [minarea=sq. m] [valleythreshold=dg] [valleysmoothing=m] [contoursvalleys=m] [contoursslopes=m] [naspectclasses=n classes] [avalanchethreshold=degrees] [hemisphere=N|S] gunits=name downstreamgunits=name slope_dg=name aspect=name avalanchearea=name gunitsglacierarea=name sunhoursprefix=name valleys=name glaciercontours=name slopeaspect=name gunitssoil=name gunitslanduse=name hydrotopes=name contours=int m [contourrast=string] [management=name] [wetland=name] [strfilepath=path] [initialdebris=name] [routingratio=name] [mswimhydrotopespath=name] [--overwrite] [--help] [--verbose] [--quiet] [--ui]
Flags:
- -d
- Skip recreating slope_dg, aspect and sunhours (must all exist)
- -s
- Only rewrite structure file, all maps must exist
- -a
- Add new soil ID for slope soils (recommended)
- -k
- Keep intermediat files (include __ in names)
- --overwrite
- Allow output files to overwrite existing files
- --help
- Print usage summary
- --verbose
- Verbose module output
- --quiet
- Quiet module output
- --ui
- Force launching GUI dialog
Parameters:
- elevation=name [required]
- Elevation raster
- accumulation=name [required]
- Accumulation raster, e.g. from m.swim.subbasins
- glacierarea=name [required]
- Glacier modelling domain raster, e.g. created using a min. elevation
- subbasins=name [required]
- Subbasin raster map
- landuse=name [required]
- Land use raster map in SWIM land use classification
- soil=name [required]
- Soil raster map classified to SWIM soil files
- minarea=sq. m
- Minimum unit area in sq. m
- Default: 100000
- valleythreshold=dg
- Valley slope threshold in dg
- Default: 19
- valleysmoothing=m
- Valley smoothing radius in m
- Default: 200
- contoursvalleys=m
- Elevation contour intervals in valleys in m
- Default: 30
- contoursslopes=m
- Elevation contour intervals on slopes in m
- Default: 200
- naspectclasses=n classes
- Number of aspect classes
- Default: 4
- avalanchethreshold=degrees
- Slope threshold to delineate avalanche area in dg. If not given, avalanchearea raster must exist.
- Default: 45
- hemisphere=N|S
- Summer/winter differentiation
- Default: N
- gunits=name [required]
- Name of glacier units raster
- Default: gunits
- downstreamgunits=name [required]
- Name of next downstream glacier unit ID raster
- Default: downstreamgunits
- slope_dg=name [required]
- Name of slope raster
- Default: slope_dg
- aspect=name [required]
- Name of aspect raster
- Default: aspect
- avalanchearea=name [required]
- Name of raster with avalanche area. Must exist if avalanchethreshold not given.
- Default: avalanchearea
- gunitsglacierarea=name [required]
- Cleaned glacier area raster fitting to gunits
- Default: gunitsglacierarea
- sunhoursprefix=name [required]
- Prefix for sun hours rasters *_summer and *_winter
- Default: sun_hr
- valleys=name [required]
- Name of valleys raster
- Default: valleys
- glaciercontours=name [required]
- Name of glacier valley and slope classed contours raster
- Default: glaciercontours
- slopeaspect=name [required]
- Name of slope aspect raster
- Default: slopeaspect
- gunitssoil=name [required]
- Name of soils raster mapped to glacier units
- Default: gunitssoil
- gunitslanduse=name [required]
- Name of landuse raster mapped to glacier units
- Default: gunitslanduse
- hydrotopes=name [required]
- Name of hydrotope raster to be created
- Default: hydrotopes
- contours=int m [required]
- Elevation contour intervals outside of glacier area
- Default: 100
- contourrast=string
- Elevation contours to be created with glacier units and contours outside of glacier area
- Default: contours
- management=name
- Management raster name (if not given, column is filled with default)
- wetland=name
- Wetlands raster name (if not given, column is filled with default)
- strfilepath=path
- Path where structure file will be written if given, glacier input files will be created in the same directory as glaciers.str.
- initialdebris=name
- Initial debris concentration (fraction of total volume), if not given set to 0
- routingratio=name
- Create raster with number of cells upstream/downstream raster
- mswimhydrotopespath=name
- Path to custom m.swim.hydrotopes script or if not installed (e.g. in testing)
DESCRIPTION
The m.swim.glaciers module creates all necessary input for the SWIM glacier dynamics module, aka SWIM-G. It creates the glacier unit input file (glaciers.str) and the hydrotope "structure" file.For a detailed description, refer to:
Wortmann, M., Bolch, T., Menz, C., Tong, J., Krysanova, V., 2018. Comparison and Correction of High-Mountain Precipitation Data Based on Glacio-Hydrological Modeling in the Tarim River Headwaters (High Asia). J. Hydrometeor. 19, 777-801. https://doi.org/10.1175/JHM-D-17-0106.1
Wortmann, M., Bolch, T., Su, B., Krysanova, V., 2019. An efficient representation of glacier dynamics in a semi-distributed hydrological model to bridge glacier and river catchment scales. Journal of Hydrology 573, 136-152. https://doi.org/10.1016/j.jhydrol.2019.03.006
Wortmann, M.S.J., 2017. Integrating glacier dynamics into hydrological modelling: A climate change impact assessment of the data-scarce headwaters of the Tarim River, Central Asia (Ph.D.). University College London, London, UK. https://discovery.ucl.ac.uk/id/eprint/10052031/
EXAMPLE
Refer to the test script of the module in the m.swim repository under:test/test_glaciers.grass
SEE ALSO
m.swim.subbasins,Documentation [external]: m.swim
AUTHORS
Michel Wortmann, Potsdam Institute for Climate Impact Research (wortmann@pik-potsdam.de)Last updated (docs/code): Tue Feb 16 11:51:42 2021 +0100 / Wed Jan 19 23:04:38 2022 +0000