R/SWATProfileDailyToVegper.R

In rhabel/modLWFB90: R-package for applying LWF-Brook90 directly from FVA-BW servers

#' Internal Helper SWATProfile.DailyToVegper
#'
#' ...
#'
#' @param ...
#'
#' @return Returns the desired output
#'
#'
#' @import data.table
#' @export


SWATProfile.DailyToVegper <- function(dat, vp.year, vp.start, vp.end) {

  setDT(dat)
  setkey(dat, YR)
  vp <- data.table::data.table(YR = vp.year, start = vp.start, end = vp.end)
  data.table::setkey(vp, YR)

  dat <- merge(dat,vp, by = "YR")
  dat_vp <- dat[which(DOY >= start & DOY <= end),list(
    #SWAT
    SWAT_prf_avg = round(mean(SWAT_prf),1),
    SWAT_prf_min = round(min(SWAT_prf),1),
    SWAT_we_avg = round(mean(SWAT_we),1),
    SWAT_we_min = round(min(SWAT_we),1),
    SWAT_0100_avg = round(mean(SWAT_0100),1),
    SWAT_0100_min = round(min(SWAT_0100),1),
    SWAT_090_avg = round(mean(SWAT_090),1),
    SWAT_090_min = round(min(SWAT_090),1),
    SWAT_060_avg = round(mean(SWAT_060),1),
    SWAT_060_min = round(min(SWAT_060),1),
    SWAT_030_avg = round(mean(SWAT_030),1),
    SWAT_030_min = round(min(SWAT_030),1),

    #AWAT
    AWAT_we_avg = round(mean(AWAT_we),1),
    AWAT_we_min = round(min(AWAT_we),1),
    AWAT_0100_avg = round(mean(AWAT_0100),1),
    AWAT_0100_min = round(min(AWAT_0100),1),
    AWAT_090_avg = round(mean(AWAT_090),1),
    AWAT_090_min = round(min(AWAT_090),1),
    AWAT_060_avg = round(mean(AWAT_060),1),
    AWAT_060_min = round(min(AWAT_060),1),
    AWAT_030_avg = round(mean(AWAT_030),1),
    AWAT_030_min = round(min(AWAT_030),1),

    #RELAWAT
    RELAWAT_we_avg = round(mean(RELAWAT_we),4),
    RELAWAT_we_min = round(min(RELAWAT_we),4),
    RELAWAT_0100_avg = round(mean(RELAWAT_0100),4),
    RELAWAT_0100_min = round(min(RELAWAT_0100),4),
    RELAWAT_090_avg = round(mean(RELAWAT_090),4),
    RELAWAT_090_min = round(min(RELAWAT_090),4),
    RELAWAT_060_avg = round(mean(RELAWAT_060),4),
    RELAWAT_060_min = round(min(RELAWAT_060),4),
    RELAWAT_030_avg = round(mean(RELAWAT_030),4),
    RELAWAT_030_min = round(min(RELAWAT_030),4),

    #PSI
    PSIlogmean_we_avg = round(mean(PSIlogmean_we),1),
    PSIlogmean_we_min = round(min(PSIlogmean_we),1),
    PSIlogmean_0100_avg = round(mean(PSIlogmean_0100),1),
    PSIlogmean_0100_min = round(min(PSIlogmean_0100),1),
    PSIlogmean_090_avg = round(mean(PSIlogmean_090),1),
    PSIlogmean_090_min = round(min(PSIlogmean_090),1),
    PSIlogmean_060_avg = round(mean(PSIlogmean_060),1),
    PSIlogmean_060_min = round(min(PSIlogmean_060),1),
    PSIlogmean_030_avg = round(mean(PSIlogmean_030),1),
    PSIlogmean_030_min = round(min(PSIlogmean_030),1),

    # #StagWat
    #StagWat-Stress
    WaterTableDepth_min = round(min(WaterTableDepth, na.rm = T),3),
    WaterTableDepth_max = round(max(WaterTableDepth, na.rm=T),3),
    WaterTableDepth_avg = round(mean(WaterTableDepth, na.rm=T),3),

    Days_SaturatedCond30cm = sum(WaterTableDepth>=-0.3),
    Days_SaturatedCond60cm = sum(WaterTableDepth>=-0.6),
    Days_SaturatedCond90cm = sum(WaterTableDepth>=-0.9),
    Durations_Saturated30cm = paste((rle(WaterTableDepth >= -0.3)$lengths[rle(WaterTableDepth >= -0.3)$values]), collapse = " "),
    Durations_Saturated60cm = paste((rle(WaterTableDepth >= -0.6)$lengths[rle(WaterTableDepth >= -0.6)$values]), collapse = " "),
    Durations_Saturated90cm = paste((rle(WaterTableDepth >= -0.9)$lengths[rle(WaterTableDepth >= -0.9)$values]), collapse = " "),
    # Days_StagWat30cm = sum(DepthStagWater>=-0.3),
    # Durations_StagWat30cm = paste((rle(DepthStagWater >= -0.3)$lengths[rle(DepthStagWater >= -0.3)$values]), collapse = " "),
    #PSI-Stress
    Days_PsiWe_lower1200 = sum(PSIlogmean_we < -1200),
    Durations_PsiWe_lower1200 = paste((rle(PSIlogmean_we < -1200)$lengths[rle(PSIlogmean_we < -1200)$values]), collapse = " "),
    Defsum_PsiWe_lower1200 = round(sum( (PSIlogmean_we+1200) * (PSIlogmean_we < -1200)),2),

    Days_Psi0100_lower1200 = sum(PSIlogmean_0100 < -1200),
    Durations_Psi0100_lower1200 = paste((rle(PSIlogmean_0100 < -1200)$lengths[rle(PSIlogmean_0100 < -1200)$values]), collapse = " "),
    Defsum_Psi0100_lower1200 = round(sum( (PSIlogmean_0100+1200) * (PSIlogmean_0100 < -1200)),2),

    Days_Psi090_lower1200 = sum(PSIlogmean_090 < -1200),
    Durations_Psi90_lower1200 = paste((rle(PSIlogmean_090 < -1200)$lengths[rle(PSIlogmean_090 < -1200)$values]), collapse = " "),
    Defsum_Psi090_lower1200 = round(sum( (PSIlogmean_090+1200) * (PSIlogmean_090 < -1200)),2),

    Days_Psi060_lower1200 = sum(PSIlogmean_060 < -1200),
    Durations_Psi60_lower1200 = paste((rle(PSIlogmean_060 < -1200)$lengths[rle(PSIlogmean_060 < -1200)$values]), collapse = " "),
    Defsum_Psi060_lower1200 = round(sum( (PSIlogmean_060+1200) * (PSIlogmean_060 < -1200)),2),

    Days_Psi030_lower1200 = sum(PSIlogmean_030 < -1200),
    Durations_Psi30_lower1200 = paste((rle(PSIlogmean_030 < -1200)$lengths[rle(PSIlogmean_030 < -1200)$values]), collapse = " "),
    Defsum_Psi030_lower1200 = round(sum( (PSIlogmean_030+1200) * (PSIlogmean_030 < -1200)),2),

    #REW-STRESS
    Days_RELAWATwe_lower40 = sum(RELAWAT_we < 0.4),
    Durations_RELAWATwe_lower40 = paste((rle(RELAWAT_we < 0.4)
                                         $lengths[rle(RELAWAT_we < 0.4)$values]), collapse = " "),
    Defsum_RELAWATwe_lower40 = round(sum((1 - (RELAWAT_we/0.4))*(RELAWAT_we < 0.4)),2),

    Days_RELAWAT0100_lower40 = sum(RELAWAT_0100 < 0.4),
    Durations_RELAWAT0100_lower40 = paste((rle(RELAWAT_0100 < 0.4)
                                           $lengths[rle(RELAWAT_0100 < 0.4)$values]), collapse = " "),
    Defsum_RELAWAT0100_lower40 = round(sum((1 - (RELAWAT_0100/0.4))*(RELAWAT_0100 < 0.4)),2),

    Days_RELAWAT090_lower40 = sum(RELAWAT_090 < 0.4),
    Durations_RELAWAT090_lower40 = paste((rle(RELAWAT_090 < 0.4)
                                          $lengths[rle(RELAWAT_090 < 0.4)$values]), collapse = " "),
    Defsum_RELAWAT090_lower40 = round(sum((1 - (RELAWAT_090/0.4))*(RELAWAT_090 < 0.4)),2),

    Days_RELAWAT060_lower40 = sum(RELAWAT_060 < 0.4),
    Durations_RELAWAT060_lower40 = paste((rle(RELAWAT_060 < 0.4)
                                          $lengths[rle(RELAWAT_060 < 0.4)$values]), collapse = " "),
    Defsum_RELAWAT060_lower40 = round(sum((1 - (RELAWAT_060/0.4))*(RELAWAT_060 < 0.4)),2),

    Days_RELAWAT030_lower40 = sum(RELAWAT_030 < 0.4),
    Durations_RELAWAT030_lower40 = paste((rle(RELAWAT_030 < 0.4)
                                          $lengths[rle(RELAWAT_030 < 0.4)$values]), collapse = " "),
    Defsum_RELAWAT030_lower40 = round(sum((1 - (RELAWAT_030/0.4))*(RELAWAT_030 < 0.4)),2)
  ),
  by = list(YR)]
  return(dat_vp)
}