Skip to content

Inconsistent calculation of total_tree_area #1481

@samsrabin

Description

@samsrabin

Originally posted by @ckoven in #1432:

The function currentPatch%UpdateTreeGrassArea() is outside the sope of this PR, but I just realized that the variable currentPatch%total_tree_area is being defined in two subtly different ways in different parts of the code. I don't understand why it should be that way, so unless someone else does, we should align the two definitions.

In currentPatch%UpdateTreeGrassArea(), the total_tree_area and total_grass_area are calculated for both canopy and understory plants, but then clipped so that they cannot exceed the patch area.

Whereas in canopy_summarization, total_tree_areais calculated only for canopy trees (and thus cannot exceed patch area, unless something is wrong).

If PPA logic is strictly rank-ordered and trees are always taller than grasses, then this should not make much of a difference. It seems to me like we should base all of these definitions only on canopy plants, so I'd advocate changing the logic in currentPatch%UpdateTreeGrassArea() to align with the current logic in canopy_summarization. But curious to understand the difference.

UpdateEffectiveWindSpeed should also be changed to ensure weights always sum to 1.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    Projects

    Status

    ❕Todo

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions