fix mistake where conductance_radiative_heat was called with the wrong parameters

Author: s-kganz

R/figures.R

@@ -202,7 +202,7 @@ fig3_regression_slopes <- function(eb_regressions) {
                hjust="left", vjust="top") +
     scale_fill_viridis_c() +
     scale_y_discrete(limits=rev) +
-    scale_x_continuous(limits=c(0.92, 1.20)) +
+    scale_x_continuous(limits=c(0.98, 1.20)) +
     facet_grid(SITE_NEON ~ .,
                scales="free_y", space="free") +
     theme(strip.placement="outside",

R/gs_gbh_sensitivity.R

@@ -18,8 +18,8 @@ energy_balance_conductance_error <- function(Tl, Ta, gs, gbH, Pa, RH,
   desat   <- esat_slope(Ta)
 
   # Conductances - don't need to calculate gbH since that is provided
-  gR    <- conductance_radiative_heat(Ta, Pa, rho, rho_mol, a_lw, 
-                                      eb_constants_$cp)
+  gR    <- conductance_radiative_heat(Ta, rho, rho_mol, a_lw, 
+                                      eb_constants_$cp, eb_constants_$sb)
   gtot  <- total_conductance(gs, gbH)
 
   # Net radiation - here we don't consider canopy effects, so assume

R/leaf_eb_numeric.R

@@ -35,8 +35,8 @@ energy_balance_error <- function(Tl, Ta, u, gs, Pa, RH,
 
   # Conductances
   gbH   <- conductance_boundary_heat_needleleaf(u, d, rho_mol)
-  gR    <- conductance_radiative_heat(Ta, Pa, rho, rho_mol, a_lw, 
-                                      eb_constants_$cp)
+  gR    <- conductance_radiative_heat(Ta, rho, rho_mol, a_lw, 
+                                      eb_constants_$cp, eb_constants_$sb)
   gtot  <- total_conductance(gs, gbH)
 
   # Net radiation