Skip to content

rothc_scenario seems to return incorrect data when scen = 'CLT' #23

@BrentWHR

Description

@BrentWHR

Code comments suggest that in the CLT scenario all cultivations are changed to either permanent grassland or winter wheat.
https://github.com/AgroCares/BLN/blob/21f2bb6272fe93c97a3f293084b3b28a96bf368b/R/bln_clim_rothc.R#L374

However, the object that will be returned still contains the EOM and HC data from the original crops.

Repex:

clt <- BLN::rothc_scenario(c(265, 256, 2014), 'CLT')
clt$rotation

This table clearly does not have the same data that I would expect, name that in:

bln_crops[crop_code == 233,.(crop_code, crop_name, B_LU_EOM, B_LU_EOM_RESIDUE, B_LU_HC)]

Furthermore, M_GREEN_TIMING is modified as in BAU. It does not consider that the crop rotation changes to permanent winter wheat.

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions