Recover the biomass_max a candidate implies

A ForCS cohort never quite reaches its biomass_max: it approaches it asymptotically while mortality is already removing biomass, so the plateau it actually holds is some fraction of the parameter. The current parameter set was built as biomass_max = reference curve maximum x 200, which assumes that fraction is 1. It is not: it runs from 0.90 (trembling aspen) to 1.00 (lodgepole pine) across these species, so that rule undershoots the intended plateau by up to 11%.

Usage

growth_inflation_factor(achieved, biomass_max, level, biomass_max_scale = 200)

Arguments

achieved

Numeric. Plateau the simulated curve actually reaches, in whatever units the curves are expressed in.

biomass_max

Numeric. The maximum-biomass parameter that simulation ran at.

level

Numeric. Plateau the curve should reach; same units as achieved.

biomass_max_scale

Numeric. Divide biomass_max by this to express it in the curve's own units. The default 200 is the ForC Succession convention: biomass_max is g m^-2 of biomass while the summary log reports g C m^-2, and 1 Mg C ha^-1 corresponds to 200 g m^-2 of biomass. For an extension whose parameter and output share units – Biomass Succession reports g m^-2 of biomass against a maxBiomass in g m^-2 – pass 1.

Value

A list with achieved_frac, inflation, and biomass_max_est.

Details

The fraction depends only on the growth and mortality shapes and on the ratio of anpp_max to biomass_max, not on the absolute level – verified across a ForC Succession sweep, where combinations sharing a ratio but differing in absolute biomass_max agree on the achieved fraction to within 0.04%. That invariance is what lets the level be recovered arithmetically instead of searched, and it is why biomass_max is held fixed across the factorial.

See also

Other growth calibration helpers: extract_landscape_cohort_structures(), growth_add_objective(), growth_auto_window(), growth_best_candidates(), growth_bin_observations(), growth_calibration_design(), growth_calibration_partition(), growth_calibration_work_root(), growth_climatic_distance(), growth_climatic_weight(), growth_expand_over_pseudo_species(), growth_factorial_ratio_grid(), growth_fitting_windows(), growth_identifiability(), growth_mortality_onset_frac(), growth_plot_palette(), growth_pseudo_species_name(), growth_reference_curves(), growth_reference_inflection(), growth_score_fit(), growth_scoring_for(), growth_smooth_observations(), growth_structure_cell_curves(), growth_structure_cohort_table(), growth_structure_design(), growth_structure_summary(), growth_window_for(), plot_growth_calibration(), plot_growth_candidate(), plot_growth_factorial_sensitivity(), plot_growth_structures(), read_growth_scoring(), read_landscape_cohort_structures(), scale_linetype_growth_reference(), write_growth_review_bundle()