exec/VWC_Diurnal_forEach_Rep.R

In bmcnellis/SDEF.analysis: Data Analysis for 2016 San Dimas Experimental Forest Project

#####VWC analysis for SDEF Project (chapter 1)
####Goals: Calculate daily max, daliy min and diurnal values for soils moisture
####Author: Michala Phillips

setwd("D:/R_Workspace/SDEF/2017_6_Analysis_scripts/Processed_VWC_Data")
adefas<-read.csv(file = " Adenostoma_VWC_2017_02_07.csv")
bare<-read.csv(file = " Bare_VWC_2017_02_07.csv")
grass<-read.csv(file = " Grass_VWC_2017_02_07.csv")


#calculates diurnal flux for ADEFAS
time<-as.POSIXlt(adefas[,2])
diurnal.time<-time[which(time$hour==12&time$min==0)]#sets all hours to noon
diurnal.vwc<-vector(mode="numeric",length=length(diurnal.time))
for (i in 1:length(diurnal.time)){
  xyear<-diurnal.time[i]$year
  xmonth<-diurnal.time[i]$mon
  day<-diurnal.time[i]$mday
  sub<-mean(adefas[which(time$year==xyear&time$mon==xmonth&time$mday==day), 5],na.rm=TRUE)
  diurnal.vwc[i]<-sub
}  #loop makes a vector of time that is y,d,m and vwc day means
diurnal.adefas<-data.frame(diurnal.time,diurnal.vwc)
colnames(diurnal.adefas)<-c("time","vwc3")
rm(day,xmonth,xyear,sub,i,time,diurnal.time,diurnal.vwc)
#adeno <- diurnal.adefas
#adeno <- merge(diurnal.adefas, adeno, by = "time")
#write.csv(adeno, file="Adenostoma_VWC_withreps_2019_04_01.csv")

#calculates diurnal flux for Grass
time<-as.POSIXlt(grass[,2])
diurnal.time<-time[which(time$hour==12&time$min==0)]#sets all hours to noon
diurnal.vwc<-vector(mode="numeric",length=length(diurnal.time))
for (i in 1:length(diurnal.time)){
  xyear<-diurnal.time[i]$year
  xmonth<-diurnal.time[i]$mon
  day<-diurnal.time[i]$mday
  sub<-mean(grass[which(time$year==xyear&time$mon==xmonth&time$mday==day), 5],na.rm=TRUE)
  diurnal.vwc[i]<-sub
}  #loop makes a vector of time that is y,d,m and vwc day means
diurnal.grass<-data.frame(diurnal.time,diurnal.vwc)
colnames(diurnal.grass)<-c("time","vwc3")
rm(day,xmonth,xyear,sub,i,time,diurnal.time,diurnal.vwc)
#g <- diurnal.grass
g <- merge(diurnal.grass, g, by = "time")

#write.csv(g, file="Grass_VWC_withreps_2019_04_01.csv")

#calculates diurnal flux for Bare
time<-as.POSIXlt(bare[,2])
diurnal.time<-time[which(time$hour==12&time$min==0)]#sets all hours to noon
diurnal.vwc<-vector(mode="numeric",length=length(diurnal.time))
for (i in 1:length(diurnal.time)){
  xyear<-diurnal.time[i]$year
  xmonth<-diurnal.time[i]$mon
  day<-diurnal.time[i]$mday
  sub<-mean(bare[which(time$year==xyear&time$mon==xmonth&time$mday==day),5],na.rm=TRUE)
  diurnal.vwc[i]<-sub
}  #loop makes a vector of time that is y,d,m and vwc day means
diurnal.bare<-data.frame(diurnal.time,diurnal.vwc)
colnames(diurnal.bare)<-c("time", "vwcbare3")
rm(day,xmonth,xyear,sub,i,time,diurnal.time,diurnal.vwc)
#b <- diurnal.bare
b <- merge(diurnal.bare, b, by = "time")
write.csv(b, file="Bare_VWC_withreps_2019_04_01.csv")