Add biomass and HI to final_stats

Added biomass, biomass_ns, harvest_index and harvest_index_adj to the final_stats output for OCI runs.

Author: chris-s-bowden

aquacrop/entities/output.py

@@ -32,6 +32,10 @@ def __init__(self, time_span, initial_th):
                 "crop Type",
                 "Harvest Date (YYYY/MM/DD)",
                 "Harvest Date (Step)",
+                'biomass',
+                'biomass_ns',
+                'harvest_index',
+                'harvest_index_adj',
                 "Dry yield (tonne/ha)",
                 "Fresh yield (tonne/ha)",
                 "Yield potential (tonne/ha)",

aquacrop/timestep/run_single_timestep.py

@@ -547,6 +547,10 @@ def solution_single_time_step(
                 Crop_Name,
                 clock_struct.step_end_time,
                 clock_struct.time_step_counter,
+                NewCond.biomass,
+                NewCond.biomass_ns,
+                NewCond.harvest_index,
+                NewCond.harvest_index_adj,
                 NewCond.DryYield,
                 NewCond.FreshYield,
                 NewCond.YieldPot,