data <- read.csv("./data/n_transformation.csv", na.strings="N.D.")
data <- na.omit(data)
meta_data <- read.csv("./data/site_metadata.csv")

d <- merge(x=data, y=meta_data, by.x=c("Site"),by.y=c("Abbr"))
plot(d$Lon[d$Depth == "0-10_cm"],d$Nmin20[d$Depth == "0-10_cm"])