update of the figure following Denis advises. (commit before I break everything for some strips)

scripts/figure_cam_target.R

@@ -22,12 +22,11 @@ TIR_adequacy <- vroom::vroom(file.path("data","cam_adequacy_target.csv"),
 summary(lm(T_target_IR~ T_target_m, data = TIR_adequacy)) -> lm_Tcam_ttarget
 ################################################################################
 #### figure camera adequacy on dates kept ####
-lm_Tcam_ttarget$coefficients
 (TIR_adequacy %>% 
    ggplot()+
    aes(x = T_target_m, y = T_target_IR)+
    geom_point(na.rm = T, alpha = .5)+
-   geom_abline(slope = 1, intercept = 0, colour = "blue", linetype = "dashed")+
+   geom_abline(slope = 1, intercept = 0, colour = "blue")+
    geom_abline(slope = lm_Tcam_ttarget$coefficients[2], intercept = lm_Tcam_ttarget$coefficients[1], colour = "red", linetype = "dashed")+
    # annotate("text",x = 5, y = 45, label = bquote({T[s] == .(round(lm_Tcam_ttarget$coefficients[2],2))*T[TIR] + .(round(lm_Tcam_ttarget$coefficients[1],2))}),
    #          size = textsize/3, alpha = 0.8, colour = " dark red", hjust = 0)+
@@ -35,8 +34,10 @@ lm_Tcam_ttarget$coefficients
             size = textsize/3, alpha = 0.8, colour = "red",hjust = 0)+
    # annotate("text",x = 5, y = 37, label = bquote({RMSE == .(format(sqrt(mean(lm_Tcam_ttarget$residuals**2)),digits = 2)) .(paste("°C"))}),
    #          size = textsize/3, alpha = 0.8, colour = " dark red",hjust = 0)+
+   theme_cowplot()+
+   theme(panel.border = element_rect(colour="black", fill=NA, size=1), axis.line = element_blank())+
    scale_y_continuous(expression(T[TIR]~paste("(°C)")), limits = c(5,55))+
-   scale_x_continuous(expression(T[thermocouples]~paste("(°C)")), limits = c(5,55))+theme_cowplot()+
+   scale_x_continuous(expression(T[thermocouples]~paste("(°C)")), limits = c(5,55))+
    annotate("text", x = 53,y=50, angle = 45, label ="1:1", size = textsize/3.5, colour ="blue")) %>% 
  
   ggsave2(filename =  file.path("graphes","figure_TIR_accuracy.jpeg"),plot = ., dpi= 300, bg="white",width= 70,height = 70,units = "mm" )

scripts/figure_daily_gs.R

@@ -52,8 +52,8 @@ if (sys.nframe() == 0) {# must be run manually if this is the file on which we a
             axis.ticks = element_line(size = .3),
             axis.ticks.length = unit(-.15, units = "cm"),
             
-            axis.title = element_text(size = textsize),
-            axis.title.x.bottom = element_blank(),
+            axis.title = element_text(size = textsize*1.1),
+            axis.title.x.bottom = element_text(size = textsize*1.1, vjust = 0),
             axis.text.x = element_text(size =textsize , vjust = -1),
             axis.text.x.top = element_blank(),
             axis.ticks.x.top = element_blank(), 
@@ -73,8 +73,8 @@ if (sys.nframe() == 0) {# must be run manually if this is the file on which we a
       scale_linetype_manual(name = expression(g[s]~ method),
                             limits = c("GT","GT'","PMT"),values = c("solid","dashed","dotdash"),
                             labels = c(expression(g[paste("s",",GT")]),expression(g[paste("s",",GT'")]),expression(g[paste("s",",PMT")])))+
-      scale_x_datetime( date_breaks = "3 hours", date_labels = "%HH", date_minor_breaks = "1 hour")+
-      scale_y_continuous(parse(text="stomatal~conductance~(m.s^-1)"),labels = fmt_dcimals(2),
+      scale_x_datetime("Time (UT+1)", date_breaks = "3 hours", date_labels = "%HH", date_minor_breaks = "1 hour")+
+      scale_y_continuous(parse(text="canopy~conductance~(m.s^-1)"),labels = fmt_dcimals(2),
                          sec.axis = dup_axis()) } %>% 
     
     ggsave2("graphes/figure_daily_gs.jpeg",., 
@@ -144,7 +144,7 @@ if (sys.nframe() == 0) {# must be run manually if this is the file on which we a
                             limits = c("GT","GT'","PMT","PMEC"),values = c("solid","dashed","dotdash","twodash"),
                             labels = c(expression(g[paste("s",",GT")]),expression(g[paste("s",",GT'")]),expression(g[paste("s",",PMT")]),expression(g[paste("s",",PMEC")])))+
       scale_x_datetime( date_breaks = "3 hours", date_labels = "%HH", date_minor_breaks = "1 hour")+
-      scale_y_continuous(parse(text="stomatal~conductance~(m.s^-1)"),labels = fmt_dcimals(2),
+      scale_y_continuous(parse(text="canopy~conductance~(m.s^-1)"),labels = fmt_dcimals(2),
                          sec.axis = dup_axis())} %>% 
     
     ggsave2("graphes/SA_daily_gs_FULL.jpeg",., 

scripts/figure_energy_canopy.R

@@ -109,8 +109,8 @@ dates_Denis %>%
         scale_y_continuous(expression(Energy~flux~(W~m^-2)),limits = c(-100,450),
                            breaks = -1:4*100,expand = expansion(add=0),
                            sec.axis = sec_axis(expression(paste(T[canopy] - T[8]," (°C)")),trans=~(.-150)/50))+
-        scale_x_datetime("Hours (UT+1)",date_labels = "%H", date_breaks = "4 hours", expand = expansion(add = .01), position ="top")+
-        scale_fill_manual(values=c(Rn_canopy ="#cccccc",LE_canopy = "#818181",S_canopy = "#333333",Tc_Ta = "red")),x=-.005,y=.65,width = 1.02,height = .35,scale=1)) %>% 
+        scale_x_datetime("Time (UT+1)",date_labels = "%H", date_breaks = "4 hours", expand = expansion(add = .01), position ="top")+
+        scale_fill_manual(values=c(Rn_canopy ="#cccccc",LE_canopy = "#818181",S_canopy = "#333333",Tc_Ta = "red")), x=-.005,y=.65,width = 1.02,height = .35,scale=1)) %>% 
   
   ggsave2(file.path("graphes","figure_energy_canopy.jpeg"),.,dpi = 500, bg= "white", width = 140,height = 100, units = "mm" )
 }