diff --git a/2_Proc.r b/2_Proc.r index 0c57691..7d2e774 100644 --- a/2_Proc.r +++ b/2_Proc.r @@ -28,7 +28,6 @@ SBD1_FROM_MAP <- read_csv("Data/SBD_Data/SBD1_Map_Wells.csv") %>% mutate(wdid=as SBD1_WELLS <- SBD1_FROM_MAP %>% full_join(SBD1_WELLS) %>% mutate(IN_EXTENT=ifelse(is.na(IN_EXTENT),0,1)) SBD1_WELLS <- SBD1_WELLS %>% filter(!is.na(SBD1_year))%>% mutate(SBD1_year=ifelse(IN_EXTENT==1,2011,SBD1_year)) %>% select(-IN_EXTENT) SBD_LINK <- SBD1_WELLS %>% full_join(SBD_LINK) %>% mutate(SBD1_year=ifelse(is.na(SBD1_year),Inf,SBD1_year),SBD1=ifelse(SBD1_year==Inf,0,1)) #Make sure both data sets conform. -WELLS ######################Static Data STATIC_DATA <- WELLS%>% left_join(SBD_LINK) %>% mutate(SBD1_year=ifelse(is.na(SBD1_year),Inf,SBD1_year)) STATIC_DATA[,-1:-5] <- STATIC_DATA[,-1:-5] %>% replace(is.na(.), 0) #Include subdistrict indicators @@ -58,7 +57,7 @@ saveRDS(STATIC_DATA,"Data/Output_Data/Well_Level_Static_Data.rds") ###########Collect Pumping Data -LIST <- WELL_DATA %>% pull(wdid) %>% unique +LIST <- WELLS %>% pull(wdid) %>% unique ALL <- length(LIST ) if(file.exists("./Data/Output_Data/Div3_Pumping_Data.csv")){LIST <-LIST[-which(LIST %in% t(read.csv("./Data/Output_Data/Div3_Pumping_Data.csv")[,1] %>% unique))]}