powered by
Extract forest floor (internal)
extract_forest_floor(a, soil_dim, th = 20, N = 500, output_path, plot)