气孔导度测试基本完成

Author: kongdd

Project_PhotoPeriod/Figure2_Stomotal_PC.R

@@ -0,0 +1,78 @@
+pacman::p_load(
+  Ipaper, data.table, dplyr, lubridate, 
+  ggplot2, gg.layers, ggrepel
+)
+
+# f = "./OUTPUT/PMLV2China_flux37_LAI_glass,ConstPC_gof.csv"
+# d = fread(f)
+indir = "./Project_PhotoPeriod/OUTPUT/"
+
+sites = c("儋州", "哀牢山", "西双版纳-Rubber", "西双版纳-EBF", "鼎湖山")
+# d[site %in% sites, ]
+
+# dat[GLASS <= -2]
+set_font()
+
+plot_Figure1 <- function(VAR="GPP", type_lai = "whit", 
+  x = "g1=false", y = "g1=true", 
+  lims = c(0.2, 1), thr = 0.02) {
+  
+  # file = "PMLV2China_flux37_LAI_glass,ConstPC_gof.csv"
+  # dat_x = fread(glue("./Project_PhotoPeriod/OUTPUT/{file}"))
+  # dat_y = fread(glue("./Project_PhotoPeriod/OUTPUT/Medlyn2011_V2/{file}"))
+  dat_x = fread(glue("Project_PhotoPeriod/OUTPUT/Yu2004_(PC_photo=true,g1=false,g0=false)/PMLV2China_flux37_LAI_{type_lai},ConstPC_gof.csv")) # 对照组
+
+  dat_y = fread(glue("Project_PhotoPeriod/OUTPUT/Yu2004_(PC_photo=true,g1=true,g0=false)/PMLV2China_flux37_LAI_{type_lai},ConstPC_gof.csv")) # 实验组
+
+  ## Medlyn2011
+  # dat_x = fread("Project_PhotoPeriod/OUTPUT/Medlyn2011_(PC_photo=true,g1=false,g0=false)/PMLV2China_flux37_LAI_whit,ConstPC_gof.csv")
+
+  # dat_y = fread("Project_PhotoPeriod/OUTPUT/Medlyn2011_(PC_photo=true,g1=false,g0=true)/PMLV2China_flux37_LAI_whit,ConstPC_gof.csv")
+
+  # dat_x = fread(glue("{indir}/PMLV2China_flux37_LAI_{type_lai},{x}_gof.csv"))
+  # dat_y = fread(glue("{indir}/PMLV2China_flux37_LAI_{type_lai},{y}_gof.csv"))
+
+  df = listk(x = dat_x, y = dat_y) %>% melt_list("type")
+  dat = df %>%
+    select(-n_valid) %>%
+    melt(c("type", "site", "var")) %>%
+    dcast(site + var + variable ~ type, value.var = "value") %>%
+    .[variable %in% c("NSE", "KGE", "R2"), ]
+
+  brks <- c(thr, Inf) %>% c(-rev(.), .)
+  pdat = dat[var == VAR & x >= lims[1] & y >= lims[1]] %>%
+    mutate(diff = y - x, diff_lev = cut(diff, brks))
+
+  p <- ggplot(pdat, aes(x, y)) +
+    geom_point() +
+    geom_point(data = pdat[abs(diff) >= thr], aes(color = diff_lev), alpha = 0.8) +
+    # geom_point(data = pdat[(WithPC - NonPC) <= -0.02,], color = "red", alpha = 0.8) +
+    # geom_text_repel(data = pdat[abs(NonPC - WithPC) >= 0.02,],
+    #   aes(label = site), hjust = 0, vjust = 0, color = "blue") +
+    geom_text_repel(
+      data = pdat[abs(diff) >= thr], aes(color = diff_lev, label = site), min.segment.length = 0.2, 
+      hjust = 0, vjust = 0, alpha = 0.8) +
+    # geom_text_repel(data = pdat[(WithPC - NonPC) <= -0.02,],
+    #   aes(label = site), hjust = 0, vjust = 0, color = "red", alpha = 0.8) +
+    geom_abline(intercept = 0, slope = 1, linetype = 2, color = "red", alpha = 0.8) +
+    coord_cartesian(xlim = lims, ylim = lims, expand = FALSE) +
+    scale_color_manual(values = c("red", "darkgreen")) +
+    theme(legend.position = "none") +
+    labs(x = x, y = y) +
+    # lims(x = lims, y = lims) +
+    facet_wrap(~variable, scales = "free") 
+  prefix = toupper(type_lai)
+  write_fig(p, glue("Figure3_{x}&{y}_{prefix}_{VAR}_v2.pdf"), 10, 4.5)
+}
+
+## LAI
+# type_lai = "whit"
+type_lai = "glass"
+# x = "ConstPC"
+# x = "NonPC"
+# y = "ConstPC"
+# y = "WithPC"
+
+# plot_Figure1("GPP", type_lai, thr = 0.02)
+plot_Figure1("ET", type_lai = "whit", lims = c(0.2, 0.9), thr = 0.02)
+plot_Figure1("ET", type_lai = "glass", lims = c(0.2, 0.9), thr = 0.02)

Project_PhotoPeriod/s1_calib_ChinaPMLV2_stomatal.jl

@@ -1,13 +1,13 @@
 include("main_pkgs.jl")
 
-function runtests(config; prefix="Medlyn2011")
+function runtests(config; prefix="Medlyn2011", type_lai="whit")
   for i in eachindex(config)
     kw = config[i]
     (; PC_photo, PC_g1, PC_g0) = kw
     _prefix_PC = "(PC_photo=$PC_photo,g1=$PC_g1,g0=$PC_g0)"
     outdir = "./Project_PhotoPeriod/OUTPUT/$(prefix)_$(_prefix_PC)"
     println(outdir)
-    r = process(; type_lai="whit", optim_PC=false, outdir, kw...)
+    r = process(; type_lai, optim_PC=false, outdir, kw...)
   end
 end
 
@@ -29,11 +29,12 @@ types_lai = ["whit", "glass"]
 stomatal = Stomatal_Medlyn2011
 configs = [
   (; stomatal, PC_photo=false, PC_g1=false, PC_g0=false), # 对照组
+  (; stomatal, PC_photo=true, PC_g1=false, PC_g0=false),    # A1
   (; stomatal, PC_photo=true, PC_g1=true, PC_g0=true),    # A1
   (; stomatal, PC_photo=true, PC_g1=true, PC_g0=false),
   (; stomatal, PC_photo=true, PC_g1=false, PC_g0=true),
 ]
-runtests(configs)
+runtests(configs; type_lai="glass")
 
 
 ## 方案2: 
@@ -43,4 +44,4 @@ configs = [
   (; stomatal, PC_photo=true, PC_g1=true, PC_g0=false),   # A1
   (; stomatal, PC_photo=true, PC_g1=false, PC_g0=false),
 ]
-runtests(configs; prefix="Yu2004")
+runtests(configs; prefix="Yu2004", type_lai="glass")