AlltdoteriidBLHofl

library("ggplot2")

Data1 <- read_csv("Allt_202021.csv") %>% 
  
  
  
  


ggplot(data = Allt_202021, aes(x=CH4, y=OH_proxy)) + geom_point(aes(colour=BLH)) + 
  ggtitle("CH4 vs. OH colored by BLH, 2019-2020") +
  xlab("CH4") + ylab("OH [mol/cm3]") +
  scale_color_gradient2(low = "darkblue",
                        mid = "red",
                        high = "gold",
                        midpoint = 750,
                        breaks=seq(0,1500,by=250),
                        limits=c(0,1500)) +
  labs(col="BLH [m]") + #color bar title
  theme(axis.text.y = element_text(size = 12),
        axis.text.x = element_text(size = 12),
        axis.title.y = element_text(size = 14),
        axis.title.x = element_text(size = 14))

xlim(-10,10) + ylim(5,43) +