################################################################################
#### Project: Lowland plant migrations alpine soil C loss
#### Title:   Function | Load soil data
#### Author:  Tom Walker (thomas.walker@usys.ethz.ch)
#### Date:    26 May 2021
#### ---------------------------------------------------------------------------

load_soil <- function(){
  # load pools
  pools <- fread("./data/from_initial_results/ready_all_data.csv") %>%
    # edit treatment and column names, remove species data
    mutate(Treatment = substr(Treatment, nchar(Treatment), nchar(Treatment))) %>%
    mutate(Site = tolower(Site)) %>%
    # add leading zero to block
    mutate(Block = ifelse(
      nchar(Block) == 3,
      Block,
      paste0(substr(Block, 1, 1), "0", substr(Block, 2, 2))
    )) %>%
    # sort column names and select variables
    dplyr::rename_with(.cols = Sample:Treatment, tolower) %>%
    select(-sample, -spp_PC1, -spp_PC2) %>%
    as.data.frame
  # load lavey extra data and combine
  lav <- fread("./data/from_initial_results/ready_lav_data.csv") %>%
    as.data.frame %>%
    select(Block:Elevation, MR.ugC.g.h:CUE) %>%
    rename_with(.cols = Block:Elevation, tolower) %>%
    mutate(site = "lavey") %>%
    # add leading zero to block
    mutate(block = ifelse(
      nchar(block) == 3,
      block,
      paste0(substr(block, 1, 1), "0", substr(block, 2, 2))
    )) %>%
    mutate(treatment = substr(treatment, nchar(treatment), nchar(treatment)))
  # combine
  out <- pools %>%
    left_join(., lav, by = c("site", "block", "treatment", "elevation")) %>%
    # setup proper units
    rename(Rm = MR.ugC.g.h, 
           Gm = MG.ugC.g.h, 
           RmM = MRm.ugCmic.g.h,
           GmM = MGm.ugCmic.g.h) %>%
    mutate(Cmic = Cmic/1000,
           CUE = CUE*100,
           RmM = RmM * 1000,
           GmM = GmM * 1000)
  # return
  return(out)
}