obs <- readRDS("~/Dropbox/Dissertacao/data2sens&calib.rds") 

if(!"albedoPAR" %in% names(obs)){
albedoPAR <- obs %>%   
  select(date,PARin,PARref) %>%
  cleanQuantis(vars = c("PARin","PARref"),qqts = c(0.01,0.99),times = "%H") %>% 
  mutate(albedoPAR = PARref/PARin) %>%
  mutate(daytime = ifelse( fCalcPotRadiation( DoY.V.n = format(date,"%j") %>% as.numeric,
                                              Hour.V.n = format(date,"%H") %>% as.numeric,
                                              Lat_deg.n = -21.61,
                                              Long_deg.n = -47.63,
                                              TimeZone_h.n = -3,
                                              useSolartime.b = TRUE) > 0 , 
                           "day","night") ) %>%
  mutate(albedoPAR = ifelse(daytime == "day",albedoPAR,NA) ) %>% 
  mutate(albedoPAR = roll_mean(albedoPAR, n = 15*24,fill = c(0),align = "center",na.rm = TRUE)) %>%
  mutate(albedoPAR = ifelse( (albedoPAR > 0.1 | albedoPAR < 0.02 ), NA, albedoPAR)) %>%
  mutate(albedoPAR = ifelse(is.na(PARref),NA,albedoPAR)) %T>%
  filter(daytime == "day") %T>% 
  timePlot(c("PARin","PARref","albedoPAR")) %>%
    select(albedoPAR)
    
obs <- cbind(obs,albedoPAR)

saveRDS(obs, "~/Dropbox/Dissertacao/data2sens&calib.rds")
} else {
    
  obs %>% 
  select(date,PARin,PARref,albedoPAR) %>%
        timePlot(c("PARin","PARref","albedoPAR"))
}

dir_out <- paste0("/home/",system("echo $USER",intern = TRUE),"/ADBHM/trash/")
infile2 <- paste0("/home/",
                  system("echo $USER",intern = TRUE),
                  "/Dropbox/Dissertacao/TORRE.txt")