##################################
rm(list = ls()) # Clean the memory
######### DEFINE THE WORKING DIRECTORY
setwd("C:\\Users\\Lenovo\\Desktop\\APPROACH1\\Inputs") # or setwd("C:/Users/Lenovo/Desktop/APPROACH1/Inputs")
##################################################################
################################################# READ THE DATABSE
######## IMPORT TTHE DATA
Data <- read.csv("Approach1.csv") # Read the database
AD <- Data[, 1] # Column 1: Activity data (area)
U_AD <- Data[, 2] # Column 2: U_Activity data (area)
Gw <- Data[, 3] # Column 3: Average annual above-ground biomass growth
U_Gw <- Data[, 4] # Column 4: U_Average annual above-ground biomass growth
R <- Data[, 5] # Column 5: Ratio of below-ground biomass to above-ground biomass
U_R <- Data[, 6] # Column 6: U_Ratio of below-ground biomass to above-ground biomass
CF <- Data[, 7] # Column 7: Carbon fraction of dry matter
U_CF <- Data[, 8] # Column 8: U_Carbon fraction of dry matter
n <- length(AD) # or length(Any variable) # Sample size
####################################################################
####################################################################
############################## EMISSION FACTORS /U_EMISSION FACTORS
EF <- Gw * (1 + R) * CF
U_EF <- sqrt(((Gw * CF) * sqrt(U_Gw^2 + U_CF^2))^2 + ((Gw * R * CF)*sqrt(U_Gw^2 + U_R^2 + U_CF^2))^2)/abs((Gw * CF) + (Gw * R * CF))
####################################################
############################## REMOVALS /U_REMOVALS
REMOVALS <- AD * EF * (-1) * (44/12)
U_REMOVALS <- sqrt(U_AD^2 + U_EF^2)
################################################################
############################## TOTAL REMOVALS /U_TOTAL REMOVALS
TOTAL_REMOVALS <- sum(REMOVALS)
U_TOTAL_REMOVALS <- sqrt(sum((REMOVALS * U_REMOVALS)^2))/abs(TOTAL_REMOVALS)
##########################################################
### CONCATENATE THE DATABASE AND THE ESTIMATED QUANTITIES
Data_Extended <- data.frame(AD, U_AD, Gw, U_Gw, R, U_R, CF, U_CF, EF, U_EF, REMOVALS, U_REMOVALS)
#####################################################
### ADD A ROW OF TOTAL_REMOVALS AND U_TOTAL_REMOVALS
TOTAL_REMOVALS_U_TOTAL_REMOVALS <- c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, TOTAL_REMOVALS, U_TOTAL_REMOVALS)
###########################################################################
### FINAL DATA FRAME #####################################################
Final.data.frame <- rbind(Data_Extended,TOTAL_REMOVALS_U_TOTAL_REMOVALS)
print(Final.data.frame)
## AD U_AD Gw U_Gw R U_R CF U_CF EF U_EF
## 1 940357.62 4.96 3.18 109.94 0.37 5.40 0.47 5.32 2.047602 85.677562
## 2 5830.74 4.96 3.18 109.94 0.37 5.40 0.47 5.32 2.047602 85.677562
## 3 11460.42 4.96 0.18 472.22 0.37 5.40 0.47 5.32 0.115902 367.549571
## 4 66148.74 4.96 9.90 5.05 0.49 4.08 0.47 5.32 6.932970 5.643982
## 5 904.77 4.96 15.00 136.67 0.37 5.40 0.47 5.32 9.658500 106.459223
## 6 NA NA NA NA NA NA NA NA NA NA
## REMOVALS U_REMOVALS
## 1 -7060086.526 85.82101
## 2 -43776.461 85.82101
## 3 -4870.381 367.58304
## 4 -1681559.843 7.51373
## 5 -32041.977 106.57471
## 6 -8822335.188 68.69603
##################################################################
#################################### FINAL RESULTS ##############
##################################### DEFINE THE OUTPUT DIRECTORY
setwd("C:\\Users\\Lenovo\\Desktop\\APPROACH1\\Outputs") # or setwd("C:/Users/Lenovo/Desktop/APPROACH1/Outputs")
# SAVE THE FINAL RESULTS
write.csv(Final.data.frame,file="Belize_Forest land remaining forest land (undisturbed).csv")
################################################### END ########################################