Write output tables for Annual Sum (total water depth per water year)

[ChristinaB]

Water year = Oct 1- Sept 30

Sum 1: For each drainage, sum daily variable by water year. Results: Rows years; Columns Drainage

Sum 2: For each Sum1, average annual sum for each variable Results: Rows Drainage; Columns 1 value per drainage

Here is example code in R:

 `yearlySummary_Precipitation_in.txt` <- Precipitation_in.txt %>%
    # group by TimeStep and Month
    select(-c(TimeStep, yyyymmdd, year, day)) %>%
    melt(., id.vars = c('month', 'water_year'), variable.name = 'Drainages', value.name = 'Precipitation') 

    # calculate total by month, wateryear, and drainage
    group_by(month, water_year, Drainages) %>%
    summarise(sumByMonthByWateryearByDrainage = sum(Precipitation)) %>%

 # calculate mean by wateryear, and drainage
    group_by(water_year, Drainages) %>%
    summarise(meanByWateryearByDrainage = sum(sumByMonthByWateryearByDrainage)) %>%

 # calculate mean 
    group_by(water_year) %>%
    mutate(meanByWateryear = mean(meanByWateryearByDrainage),

# reformat the table to wide for each DrainageID
    dcast(., water_year + meanByWateryear ~ Drainages, value.var = 'meanByWateryearByDrainage') %>%
    data.table()`

Variables= Precip, ET, Recharge (Rd), Upwelling (r3), Return flow

and User Withdrawal

Slight change for users, group by user type

loop through user = 1:6

User Sum 1: For each drainage, sum daily variable by water year. Results: Rows years; Columns Drainage, append Source Mixing ID, ReturnFlowID

User Sum 2: For each Sum1, average annual sum for each variable Results: Rows Drainage; Columns 1 value per drainage, append Source Mixing ID, ReturnFlowID

print user withdrawal for each user type

`test_summary_wtryr <- UserWithdrawal_gpd.txt %>% 
    select(subset_usercols) %>%
    select(-c(TimeStep, yyyymmdd, month, year, day)) %>%
    melt(., id.vars = c('water_year'), variable.name = 'UserIDs', value.name = 'UserWithdrawal') %>%
    group_by(water_year, UserIDs) %>%
    summarise(annualUserWithdrawal = sum(UserWithdrawal))

test_usersWithdrawal_wtryr <- test_users %>%
  merge(., test_summary_wtryr, by='UserIDs')

fwrite(test_usersWithdrawal_wtryr, paste0(getwd(),'/AverageWithdrawalwateryearly_gpy_sourcesummary.txt'), sep = '\t', row.names = F, col.names = T, append = F)`

Ideal future outcome: Annual average sum of water depth Total P, ET, Rd, R2 per drainage
Total of user Groundwater withdrawal per drainage Total of user Surface water withddrawal per drainage Total return flow per drainage