computeCritCutLength: Compute critical crack length

Description

This function implements Bettina Richter's (2019) parametrization for the critical crack length for flat simulations based on density, grain size, and shear strength. The parametrization also needs the mean density of the slab, which can be computed automatically if a snowprofile object is provided. In case the functions gets a snowprofileLayers object it expects slab_rho being precomputed. This acts as a safety mechanism to ensure that slab_rho is computed over one profile and not over a stacked layers data.frame containing multiple profiles. Note that the critical crack length can be computed alongside the layer probabilities for instability p_unstable in computePunstable.

Usage

computeCritCutLength(x)
# S3 method for snowprofileSet
computeCritCutLength(x)
# S3 method for snowprofile
computeCritCutLength(x)
# S3 method for snowprofileLayers
computeCritCutLength(x)

Value

Input object is returned with $crit_cut_length (and potentially $slab_rho) appended to the layers object.

Arguments

x: sarp.snowprofile::snowprofileSet, sarp.snowprofile::snowprofile, or sarp.snowprofile::snowprofileLayers object

Methods (by class)

computeCritCutLength(snowprofileSet): for sarp.snowprofile::snowprofileSets
computeCritCutLength(snowprofile): for sarp.snowprofile::snowprofiles
computeCritCutLength(snowprofileLayers): for sarp.snowprofile::snowprofileLayers

Author

fherla based on the python function by smayer and brichter

References

Richter, B., Schweizer, J., Rotach, M. W., & Van Herwijnen, A. (2019). Validating modeled critical crack length for crack propagation in the snow cover model SNOWPACK. The Cryosphere, 13(12), 3353–3366. https://doi.org/10.5194/tc-13-3353-2019