A further canopy productivity model (Two leaf) is desired to complement existing canopy models under development and further extend pyrealm functionality.

The proposed two leaf model is to be migrated from an existing R script and integrated with Pyrealm code and further canopy classes. The R script leans on libraries that are not available in Python, therefore key functionality must be implemented using existing Python packages or with custom implementations.

It is proposed that the existing Pyrealm Splash solar radiation implementation can be modified to provide some of the necessary functionality for the new canopy class.

Key tasks:

• refactor pyrealm.splash.solar to break out key calculations into pyrealm.core.solar
• Define new radiation calculation methodology for subdaily radiation estimates for two leaf canopy model
• Define two leaf canopy model
• Integrate two leaf canopy model with pyrealm single leaf canopy productivity model