R/init_mylake.R

In aemon-j/MyLakeR: Work with MyLake in R

init_mylake <- function(model = 'MyLake', wtemp_file, date, folder = '.',
                        config_dat = "mylake_config_final.Rdata",
                        init_dat = "mylake_init.Rdata"){
  
  setwd(folder)
  obs <- read.csv(wtemp_file)
  
  dat = which(obs[,1] == date)
  ndeps = length(dat)
  deps = obs[dat,2]
  tmp = obs[dat,3]
  deps <- signif(deps,4)
  tmp <- signif(tmp,4)
  
  if("MyLake" %in% model){
    
    if(!dir.exists('MyLake')){
      dir.create('MyLake')
    }
    
    load(file.path(".","MyLake", config_dat))
    
    mylake_init <- list()
    
    deps_Az <- data.frame("Depth_meter"=mylake_config[["In.Z"]],
                          "Az"=mylake_config[["In.Az"]])
    
    temp_interp1 <- dplyr::full_join(deps_Az,
                                     data.frame("Depth_meter"=deps,
                                                "Water_Temperature_celsius"=tmp))
    temp_interp2 <- dplyr::arrange(temp_interp1,Depth_meter)
    temp_interp3 <- dplyr::mutate(temp_interp2,
                                  TempInterp=approx(x=Depth_meter,
                                                    y=Water_Temperature_celsius,
                                                    xout=Depth_meter,
                                                    yleft=dplyr::first(na.omit(Water_Temperature_celsius)),
                                                    yright=dplyr::last(na.omit(Water_Temperature_celsius)))$y)
    temp_interp <- dplyr::filter(temp_interp3, !is.na(Az))
    
    mylake_init[["In.Tz"]]=as.matrix(temp_interp$TempInterp)
    
    mylake_init[["In.Z"]]=as.matrix(temp_interp$Depth_meter)
    
    mylake_config[["Phys.par"]][1]=median(diff(mylake_init$In.Z))
    
    save(mylake_config,file=file.path(folder,"MyLake", config_dat))
    
    save(mylake_init,file=file.path(folder,"MyLake", init_dat))
    
  }
}