ORCHIDEE CARBON FLUXES

NPP = GPP - MAINT_RESP - GROWTH_RESP

NEP = NPP - HET_RESP - CO2_FIRE = GPP - MAINT_RESP - GROWTH_RESP - HET_RESP - CO2_FIRE

NEE = -NEP = -GPP + MAINT_RESP + GROWTH_RESP + HET_RESP + CO2_FIRE (positive is source to atm)

NBP = NEP - fLuc - fWoodharvest - fHarvest = GPP - MAINT_RESP - GROWTH_RESP - HET_RESP - CO2_FIRE - fLuc - fWoodharvest - fHarvest (positive is sink to land)

fLuc = CONVFLUX + CFLUX_PROD10 + CFLUX_PROD100

fWoodharvest = CONVFLUX_HARVEST + CFLUX_PROD10_HARVEST + CFLUX_PROD100_HARVEST

fHarvest = HARVEST_ABOVE

NBP_forest = ((GPP - MAINT_RESP - GROWTH_RESP - HET_RESP) * VEGET_MAX).sum(forest) - CONVFLUX * VEGET_MAX.sum(forest) / VEGET_MAX.sum(total) - CFLUX_PROD10 - CFLUX_PROD100 - CONVFLUX_HARVEST - CFLUX_PROD10_HARVEST - CFLUX_PROD100_HARVEST

NBP_grass = ((GPP - MAINT_RESP - GROWTH_RESP - HET_RESP) * VEGET_MAX).sum(grass) - CONVFLUX * VEGET_MAX.sum(grass) / VEGET_MAX.sum(total)

NBP_crops = ((GPP - MAINT_RESP - GROWTH_RESP - HET_RESP) * VEGET_MAX).sum(crops) - CONVFLUX * VEGET_MAX.sum(crops) / VEGET_MAX.sum(total) - HARVEST_ABOVE