manuscript accepted

robertladwig · Sep 13, 2021 · 0ca15e4 · 0ca15e4
1 parent e3970f7
commit 0ca15e4
Show file tree

Hide file tree

Showing 31 changed files with 325,408 additions and 0 deletions.
diff --git a/2_createFigures_HD.R b/2_createFigures_HD.R
@@ -347,6 +347,7 @@ g <- g_density / g_icestrat_scaled / g_mixing_hd / g_lakenumber / g_summerstrat
  theme(legend.position = 'bottom', plot.tag = element_text(size = 8)); g
 
 ggsave('figs_HD/scenarios_Ensemble.png', g, dpi = 500, width = 165,height = 250, units = 'mm')
+ggsave('figs_HD/figure05.pdf', g, dpi = 500, width = 165,height = 250, units = 'mm')
 
 
 # Approach
@@ -521,6 +522,7 @@ w3 = plotWeather(df = df_mo20, icedate = '2021-03-22', usetitle = 'C) Monona 202
 w1 / g1 / w2 / g2 / w3/ g3 + plot_layout(heights = c(1,3,1,3,1,3))
 
 ggsave('figs_HD/fieldmonitoring_model2.png', dpi = 700, width = 3.25,height = 5, units = 'in')
+ggsave('figs_HD/figure03.pdf', dpi = 700, width = 3.25,height = 5, units = 'in')
 
 
 ### Other plots 

diff --git a/3_plotModelOutputs_plots.R b/3_plotModelOutputs_plots.R
@@ -90,6 +90,7 @@ me.t.glm.1.5 + mo.t.glm.1.5 +
 me.s.glm.1.5 + mo.s.glm.1.5 +
 plot_layout(design = layout, guides = 'collect')
 ggsave(paste0('figs_HD/plotModelOutputs_',startyear,'_GLM.png'), width = 6.5, height = 5, units = 'in', dpi = 500)
+ggsave(paste0('figs_HD/figureS03.pdf'), width = 6.5, height = 5, units = 'in', dpi = 500)
 
 #### GOTM #####
 me.t.GOTM.NULL = plot.output(read.output('output_modelruns/NULL/Mendota_GOTM_temp.csv', startyear), 
@@ -142,6 +143,7 @@ me.obs + mo.obs +
  me.s.GOTM.1.5 + mo.s.GOTM.1.5 +
  plot_layout(design = layout, guides = 'collect')
 ggsave(paste0('figs_HD/plotModelOutputs_',startyear,'_GOTM.png'), width = 6.5, height = 5, units = 'in', dpi = 500)
+ggsave(paste0('figs_HD/figureS04.pdf'), width = 6.5, height = 5, units = 'in', dpi = 500)
 
 #### Simstrat #####
 me.t.Simstrat.NULL = plot.output(read.output('output_modelruns/NULL/Mendota_Simstrat_temp.csv', startyear), 
@@ -194,6 +196,7 @@ me.t.Simstrat.NULL + mo.t.Simstrat.NULL +
  me.s.Simstrat.1.5 + mo.s.Simstrat.1.5 +
  plot_layout(design = layout, guides = 'collect')
 ggsave(paste0('figs_HD/plotModelOutputs_',startyear,'_Simstrat.png'), width = 6.5, height = 5, units = 'in', dpi = 500)
+ggsave(paste0('figs_HD/figureS05.pdf'), width = 6.5, height = 5, units = 'in', dpi = 500)
 
 
 

diff --git a/analytical/empirSlt.R b/analytical/empirSlt.R
@@ -282,4 +282,5 @@ g2 <- ggplot(m.out.contour, aes(x = Salt, y = Wind[value], col = variable)) +
 
 g <- g1 + g2 + plot_annotation(tag_levels = 'A') ;g 
 ggsave(file = paste0(name_lake,'_theoretical.png'), g, dpi = 500, width = 6.5, height = 2.5, units = 'in')
+ggsave(file = paste0(name_lake,'_theoretical.eps'), g, dpi = 500, width = 6.5, height = 2.5, units = 'in')