alex/CREP_Policy_Interaction_Analysis
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Minor varible rename

Browse Source
This commit is contained in:
alex 2026-05-15 16:27:23 -06:00
parent 9787e54ea6
commit 8f386f0057
1 changed files with 1 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
2_Proc.r
View File

@ -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_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) 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. 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
STATIC_DATA <- WELLS%>% left_join(SBD_LINK) %>% mutate(SBD1_year=ifelse(is.na(SBD1_year),Inf,SBD1_year)) 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 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 ###########Collect Pumping Data
LIST <- WELL_DATA %>% pull(wdid) %>% unique LIST <- WELLS %>% pull(wdid) %>% unique
ALL <- length(LIST ) 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))]} 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))]}