From b82e2de6c766329401d51e1262ec1db2e83a2887 Mon Sep 17 00:00:00 2001
From: Keith Oleson <oleson@ucar.edu>
Date: Fri, 21 Jun 2024 09:58:31 -0600
Subject: [PATCH 1/4] Initial commit to branch, executable permission for
 scripts

---
 bin/ilamb-fetch | 0
 bin/ilamb-mean  | 0
 bin/ilamb-run   | 0
 bin/ilamb-setup | 0
 4 files changed, 0 insertions(+), 0 deletions(-)
 mode change 100644 => 100755 bin/ilamb-fetch
 mode change 100644 => 100755 bin/ilamb-mean
 mode change 100644 => 100755 bin/ilamb-run
 mode change 100644 => 100755 bin/ilamb-setup

diff --git a/bin/ilamb-fetch b/bin/ilamb-fetch
old mode 100644
new mode 100755
diff --git a/bin/ilamb-mean b/bin/ilamb-mean
old mode 100644
new mode 100755
diff --git a/bin/ilamb-run b/bin/ilamb-run
old mode 100644
new mode 100755
diff --git a/bin/ilamb-setup b/bin/ilamb-setup
old mode 100644
new mode 100755

From 5d0d616c2594f367c49c0157541ac25af2d31bef Mon Sep 17 00:00:00 2001
From: Keith Oleson <oleson@ucar.edu>
Date: Fri, 21 Jun 2024 10:02:42 -0600
Subject: [PATCH 2/4] Add alt_var=TSOI to permafrost confrontation

---
 src/ILAMB/ConfPermafrost.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/ILAMB/ConfPermafrost.py b/src/ILAMB/ConfPermafrost.py
index 95ec261d..8baa52b9 100644
--- a/src/ILAMB/ConfPermafrost.py
+++ b/src/ILAMB/ConfPermafrost.py
@@ -191,7 +191,7 @@ def stageData(self, m):
         dmax = float(self.keywords.get("dmax", 3.5))
         Teps = float(self.keywords.get("Teps", 273.15))
         tsl = m.extractTimeSeries(
-            "tsl", initial_time=(y0 - 1850) * 365, final_time=(yf - 1850) * 365
+            "tsl", alt_vars = ["TSOI"], initial_time=(y0 - 1850) * 365, final_time=(yf - 1850) * 365
         )
         tsl = tsl.trim(lat=[max(obs.lat.min(), tsl.lat.min()), 90])
 

From 17a92960a7a11fd6f0e7aa497cfbafc1281b169b Mon Sep 17 00:00:00 2001
From: Keith Oleson <oleson@ucar.edu>
Date: Fri, 21 Jun 2024 10:06:04 -0600
Subject: [PATCH 3/4] Add CLM variables to soil carbon confrontation

---
 src/ILAMB/ConfSoilCarbon.py | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/ILAMB/ConfSoilCarbon.py b/src/ILAMB/ConfSoilCarbon.py
index 1476291b..2933690c 100644
--- a/src/ILAMB/ConfSoilCarbon.py
+++ b/src/ILAMB/ConfSoilCarbon.py
@@ -114,7 +114,7 @@ def confront(self, m):
         mod_soilc = (
             m.extractTimeSeries(
                 "cSoilAbove1m",
-                alt_vars=["soilc", "cSoil"],
+                alt_vars=["soilc", "cSoil", "TOTSOMC_1m"],
                 initial_time=t0,
                 final_time=tf,
             )
@@ -123,18 +123,18 @@ def confront(self, m):
         )
         mod_npp = (
             m.extractTimeSeries(
-                "npp", initial_time=t0, final_time=tf, expression="gpp-ra"
+                "npp", alt_vars=["NPP"], initial_time=t0, final_time=tf, expression="gpp-ra"
             )
             .integrateInTime(mean=True)
             .convert("kg m-2 yr-1")
         )
         mod_tas = (
-            m.extractTimeSeries("tas", initial_time=t0, final_time=tf)
+            m.extractTimeSeries("tas", alt_vars=["TSA"], initial_time=t0, final_time=tf)
             .integrateInTime(mean=True)
             .convert("degC")
         )
         mod_pr = (
-            m.extractTimeSeries("pr", initial_time=t0, final_time=tf)
+            m.extractTimeSeries("pr", initial_time=t0, final_time=tf, expression   = "RAIN+SNOW")
             .integrateInTime(mean=True)
             .convert("mm yr-1")
         )

From d3f502a47d9d0cbff462e067269bc2587b6ef3b3 Mon Sep 17 00:00:00 2001
From: Keith Oleson <oleson@ucar.edu>
Date: Fri, 21 Jun 2024 10:08:01 -0600
Subject: [PATCH 4/4] Add CLM SP and BGC configuration files

---
 src/ILAMB/data/ilamb_nohoff_final_CLM.cfg    | 836 +++++++++++++++++++
 src/ILAMB/data/ilamb_nohoff_final_CLM_SP.cfg | 610 ++++++++++++++
 2 files changed, 1446 insertions(+)
 create mode 100755 src/ILAMB/data/ilamb_nohoff_final_CLM.cfg
 create mode 100755 src/ILAMB/data/ilamb_nohoff_final_CLM_SP.cfg

diff --git a/src/ILAMB/data/ilamb_nohoff_final_CLM.cfg b/src/ILAMB/data/ilamb_nohoff_final_CLM.cfg
new file mode 100755
index 00000000..63c26346
--- /dev/null
+++ b/src/ILAMB/data/ilamb_nohoff_final_CLM.cfg
@@ -0,0 +1,836 @@
+# This configure file uses observational data which can be obtained by
+# running the following command after exporting ILAMB_ROOT to the
+# appropriate location.
+#
+#   ilamb-fetch --remote_root https://www.ilamb.org/ILAMB-Data
+#
+# Extra ilamb-run options are defined here with the '#!' symbols:
+#
+#! define_regions = DATA/regions/GlobalLand.nc
+#
+[h1: Ecosystem and Carbon Cycle]
+bgcolor = "#ECFFE6"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Biomass]
+variable       = "TOTVEGC"
+alternate_vars = "biomass,cVeg"
+weight         = 5
+skip_rmse      = True
+mass_weighting = True
+
+[Tropical]
+source     = "DATA/biomass/Tropical/biomass_0.5x0.5.nc"
+weight     = 20
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+
+[GlobalCarbon]
+source     = "DATA/biomass/GLOBAL.CARBON/biomass_0.5x0.5.nc"
+weight     = 16
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+
+[NBCD2000]
+source     = "DATA/biomass/NBCD2000/biomass_0.5x0.5.nc"
+weight     = 8
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+
+[USForest]
+source     = "DATA/biomass/USForest/biomass_0.5x0.5.nc"
+weight     = 8
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+
+[Thurner]
+source = "DATA/biomass/Thurner/biomass_0.5x0.5.nc"
+weight     = 20
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+
+[ESACCI]
+source = "DATA/biomass/ESACCI/biomass.nc"
+weight     = 20
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Burned Area]
+ctype          = "ConfBurntArea"
+variable       = "FAREA_BURNED"
+alternate_vars = "burntArea"
+weight         = 4
+cmap           = "OrRd"
+mass_weighting = True
+
+[GFED4.1S]          
+source        = "DATA/burntArea/GFED4.1S/burntArea.nc"
+weight        = 20
+relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Carbon Dioxide]
+variable      = "co2"
+ctype         = "ConfCO2"
+weight        = 5
+
+[NOAA.emulated]
+source        = "DATA/co2/NOAA.GMD/co2.nc"
+emulated_flux = "NBP"
+sites         = "alt,asc,azr,bhd,bmw,brw,cba,cgo,chr,crz,gmi,hba,ice,key,kum,mhd,mid,pocs35,pocs30,pocs25,pocs20,pocs15,pocs10,pocs05,poc000,pocn05,pocn10,pocn15,pocn20,pocn25,psa,rpb,sey,shm,smo,spo,syo,zep"
+lat_bands     = "-90,-60,-23,0,+23,+60,+90"
+force_emulation = True
+
+[HIPPOAToM]
+ctype      = "ConfGSNF"
+variable   = "gsnf"
+source     = "DATA/co2/HIPPOATom/HIPPO_ATom.nc"
+model_flux = "AR+HR-GPP"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Gross Primary Productivity]
+variable       = "GPP"
+alternate_vars = "gpp"
+cmap           = "Greens"
+weight         = 5
+mass_weighting = True
+
+[FLUXNET2015]
+source     = "DATA/gpp/FLUXNET2015/gpp.nc"
+weight     = 9
+table_unit = "g m-2 d-1"
+plot_unit  = "g m-2 d-1"
+
+[FLUXCOM]
+source        = "DATA/gpp/FLUXCOM/gpp.nc"
+weight        = 15
+table_unit    = "Pg yr-1"
+plot_unit     = "g m-2 d-1"
+space_mean    = False
+skip_iav      = True
+relationships = "Evapotranspiration/GLEAMv3.3a","Precipitation/GPCPv2.3","SurfaceDownwardSWRadiation/CERESed4.1","SurfaceNetSWRadiation/CERESed4.1","SurfaceAirTemperature/CRU4.02"
+
+[WECANN]
+source        = "DATA/gpp/WECANN/gpp.nc"
+weight        = 15
+table_unit    = "Pg yr-1"
+plot_unit     = "g m-2 d-1"
+space_mean    = False
+
+[GBAF]
+source        = "DATA/gpp/GBAF/gpp_0.5x0.5.nc"
+weight        = 15
+table_unit    = "Pg yr-1"
+plot_unit     = "g m-2 d-1"
+space_mean    = False
+skip_iav      = True
+relationships = "Evapotranspiration/GLEAMv3.3a","Precipitation/GPCPv2.3","SurfaceDownwardSWRadiation/CERESed4.1","SurfaceNetSWRadiation/CERESed4.1","SurfaceAirTemperature/CRU4.02"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Leaf Area Index]
+variable       = "TLAI"
+alternate_vars = "lai"
+cmap           = "Greens"
+weight         = 3
+mass_weighting = True
+
+[AVHRR]
+source        = "DATA/lai/AVHRR/lai_0.5x0.5.nc"
+weight        = 15
+relationships = "Precipitation/GPCPv2.3"
+
+[AVH15C1]
+source        = "DATA/lai/AVH15C1/lai.nc"
+weight        = 15
+relationships = "Precipitation/GPCPv2.3"
+
+[MODIS]
+source        = "DATA/lai/MODIS/lai_0.5x0.5.nc"
+weight        = 15
+relationships = "Precipitation/GPCPv2.3"
+skip_iav      = True
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Global Net Ecosystem Carbon Balance]
+variable = "NBP"
+alternate_vars = "nbp"
+weight   = 5
+ctype    = "ConfNBP"
+
+[GCP]      
+source   = "DATA/nbp/GCP/nbp_1959-2016.nc"
+weight   = 20
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Net Ecosystem Exchange]
+variable       = "NEE"
+alternate_vars = "nee"
+derived        = "gpp-ra-rh,GPP-AR-HR"
+weight         = 5
+mass_weighting = True
+
+[FLUXNET2015]
+source     = "DATA/nee/FLUXNET2015/nee.nc"
+weight     = 9
+table_unit = "g m-2 d-1"
+plot_unit  = "g m-2 d-1"
+
+[GBAF]
+source     = "DATA/nee/GBAF/nee_0.5x0.5.nc"
+weight     = 4
+table_unit = "Pg yr-1"
+plot_unit  = "g m-2 d-1"
+space_mean = False
+skip_iav   = True
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Ecosystem Respiration]
+variable       = "reco"
+derived        = "AR+HR"
+weight         = 4
+mass_weighting = True
+
+[FLUXNET2015]
+source     = "DATA/reco/FLUXNET2015/reco.nc"
+weight     = 6
+table_unit = "g m-2 d-1"
+plot_unit  = "g m-2 d-1"
+
+[FLUXCOM]
+source     = "DATA/reco/FLUXCOM/reco.nc"
+weight     = 4
+table_unit = "Pg yr-1"
+plot_unit  = "g m-2 d-1"
+space_mean = False
+skip_iav   = True
+
+[GBAF]
+source     = "DATA/reco/GBAF/reco_0.5x0.5.nc"
+weight     = 4
+table_unit = "Pg yr-1"
+plot_unit  = "g m-2 d-1"
+space_mean = False
+skip_iav   = True
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Soil Carbon]
+variable       = "TOTSOMC_1m"
+alternate_vars = "cSoil,cSoilAbove1m"
+weight         = 5
+mass_weighting = True
+
+[HWSD]
+source     = "DATA/soilc/HWSD_M/soilc_0.5x0.5.nc"
+weight     = 15
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+skip_rmse  = True
+
+[NCSCDV22]
+source     = "DATA/cSoil/NCSCDV22/soilc_0.5x0.5.nc"
+weight     = 12
+table_unit = "Pg"
+plot_unit  = "kg m-2"
+space_mean = False
+skip_rmse  = True
+
+[Koven]
+ctype        = "ConfSoilCarbon"
+source       = "DATA/cSoil/NCSCDV22/soilc_0.5x0.5.nc"
+weight       = 15
+soilc_source = "DATA/cSoil/NCSCDV22/soilc_0.5x0.5.nc, DATA/cSoil/HWSD/soilc_0.5x0.5.nc"
+tas_source   = "DATA/tas/CRU4.02/tas.nc"
+pr_source    = "DATA/pr/GPCCv2018/pr.nc"
+npp_source   = "DATA/cSoil/Koven/npp_0.5x0.5.nc"
+pet_source   = "DATA/cSoil/Koven/pet_0.5x0.5.nc"
+fracpeat_source = "DATA/cSoil/Koven/fracpeat_0.5x0.5.nc"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Nitrogen Fixation]
+variable = "FFIX_TO_SMINN"
+alternate_vars = "fBNF"
+cmap = "Greens"
+weight = 3
+
+[Davies-Barnard]
+source     = "DATA/fBNF/DaviesBarnard/fBNF_0.5x0.5.nc"
+table_unit = "Tg yr-1"
+plot_unit  = "kg ha-1 yr-1"
+space_mean = False
+weight     = 16
+
+###########################################################################
+
+[h1: Hydrology Cycle]
+bgcolor = "#E6F9FF"
+
+[h2: Evapotranspiration]
+variable       = "et"
+alternate_vars = "evspsbl"
+derived        = "QVEGE+QVEGT+QSOIL"
+cmap           = "Blues"
+weight         = 5
+mass_weighting = True
+
+[GLEAMv3.3a]
+source        = "DATA/evspsbl/GLEAMv3.3a/et.nc"
+weight        = 15
+table_unit    = "mm d-1"
+plot_unit     = "mm d-1"
+relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02"
+
+[MODIS]
+source        = "DATA/evspsbl/MODIS/et_0.5x0.5.nc"
+weight        = 15
+table_unit    = "mm d-1"
+plot_unit     = "mm d-1"
+relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02"
+
+[MOD16A2]
+source        = "DATA/evspsbl/MOD16A2/et.nc"
+weight        = 15
+table_unit    = "mm d-1"
+plot_unit     = "mm d-1"
+relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Evaporative Fraction]
+variable       = "EvapFrac"
+weight         = 5
+mass_weighting = True
+ctype          = "ConfEvapFraction"
+cmap           = "BuGn"
+
+[FLUXCOM]
+hfss_source = "DATA/hfss/FLUXCOM/sh.nc"
+hfls_source = "DATA/hfls/FLUXCOM/le.nc"
+weight      = 9
+skip_rmse   = True
+
+[CLASS]
+variable    = "evapfrac"
+hfss_source = "DATA/hfss/CLASS/hfss.nc"
+hfls_source = "DATA/hfls/CLASS/hfls.nc"
+skip_rmse   = True
+weight      = 25
+
+[WECANN]
+hfss_source = "DATA/hfss/WECANN/hfss.nc"
+hfls_source = "DATA/hfls/WECANN/hfls.nc"
+weight      = 9
+skip_rmse   = True
+
+[GBAF]
+hfss_source     = "DATA/sh/GBAF/sh_0.5x0.5.nc"
+hfls_source     = "DATA/le/GBAF/le_0.5x0.5.nc"
+weight     = 9
+skip_rmse  = True
+skip_iav   = True
+limit_type = "99per"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Latent Heat]
+variable       = "hfls"
+alternate_vars = "le,EFLX_LH_TOT"
+cmap           = "Oranges"
+weight         = 5
+mass_weighting = True
+
+[FLUXNET2015]
+source   = "DATA/hfls/FLUXNET2015/hfls.nc"
+weight   = 3
+
+[FLUXCOM]
+source   = "DATA/hfls/FLUXCOM/le.nc"
+land     = True
+weight   = 9
+skip_iav = True
+
+[DOLCE]
+source   = "DATA/evspsbl/DOLCE/DOLCE.nc"
+weight   = 15
+land     = True
+
+[CLASS]
+source     = "DATA/hfls/CLASS/hfls.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 25
+
+[WECANN]
+source     = "DATA/hfls/WECANN/hfls.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 9
+
+[GBAF]
+source   = "DATA/le/GBAF/le_0.5x0.5.nc"
+land     = True
+weight   = 9
+skip_iav = True
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Runoff]
+variable       = "runoff"
+alternate_vars = "mrro,QRUNOFF"
+weight         = 5
+mass_weighting = True
+
+[Dai]
+ctype          = "ConfRunoff"
+source         = "DATA/mrro/Dai/runoff.nc"
+weight         = 15
+
+[LORA]
+source         = "DATA/mrro/LORA/LORA.nc"
+table_unit     = "mm d-1"
+plot_unit      = "mm d-1"
+weight         = 15
+
+[CLASS]
+source         = "DATA/mrro/CLASS/mrro.nc"
+plot_unit      = "mm d-1"
+table_unit     = "mm d-1"
+weight         = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Sensible Heat]
+variable       = "hfss"
+alternate_vars = "sh,FSH"
+weight         = 2
+mass_weighting = True
+
+[FLUXNET2015]
+source   = "DATA/hfss/FLUXNET2015/hfss.nc"
+weight   = 9
+
+[FLUXCOM]
+source   = "DATA/hfss/FLUXCOM/sh.nc"
+weight   = 15
+skip_iav = True
+
+[CLASS]
+source     = "DATA/hfss/CLASS/hfss.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 25
+
+[WECANN]
+source     = "DATA/hfss/WECANN/hfss.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 15
+
+[GBAF]
+source   = "DATA/sh/GBAF/sh_0.5x0.5.nc"
+weight   = 15
+skip_iav = True
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Terrestrial Water Storage Anomaly]
+variable       = "twsa"
+alternate_vars = "tws,TWS"
+derived        = "RAIN+SNOW-QVEGE-QVEGT-QSOIL-QRUNOFF"
+cmap           = "Blues"
+weight         = 5
+ctype          = "ConfTWSA"
+
+[GRACE]
+source     = "DATA/twsa/GRACE/twsa_0.5x0.5.nc"
+weight     = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Snow Water Equivalent]
+variable       = "H2OSNO"
+alternate_vars = "swe,snw"
+cmap           = "Blues"
+weight         = 5
+ctype          = "ConfSWE"
+
+[CanSISE]
+source     = "DATA/swe/CanSISE/swe.nc"
+weight     = 25
+plot_unit  = "cm"
+table_unit = "cm"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Permafrost]
+variable = "TSOI"
+alternate_vars = "tsl"
+
+[Brown2002]
+ctype    = "ConfPermafrost"
+source   = "DATA/permafrost/Brown2002/Brown2002.nc"
+y0       = 1985.
+yf       = 2005.
+Teps     = 273.15
+dmax     = 3.5
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[Obu2018]
+ctype    = "ConfPermafrost"
+source   = "DATA/permafrost/Obu2018/Obu2018.nc"
+y0       = 2000.
+yf       = 2016.
+Teps     = 273.15
+dmax     = 3.5
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Soil Moisture]
+variable = "SOILWATER_10CM"
+alternate_vars = "mrsos"
+weight   = 3
+cmap     = "Blues"
+
+[WangMao]
+source   = "DATA/mrsos/WangMao/mrsos_olc.nc"
+weight   = 15
+
+###########################################################################
+
+[h1: Radiation and Energy Cycle]
+bgcolor = "#FFECE6"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Albedo]
+variable = "albedo"
+weight   = 1
+ctype    = "ConfAlbedo"
+
+[CERESed4.1]
+source   = "DATA/albedo/CERESed4.1/albedo.nc"
+weight   = 20
+
+[GEWEX.SRB]
+source   = "DATA/albedo/GEWEX.SRB/albedo_0.5x0.5.nc"
+weight   = 20
+
+#[MODIS]
+#source   = "DATA/albedo/MODIS/albedo_0.5x0.5.nc"
+#weight   = 20
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Upward SW Radiation]
+variable = "rsus"
+alternate_vars = "FSR"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rsus/CERESed4.1/rsus.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rsus/FLUXNET2015/rsus.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rsus/GEWEX.SRB/rsus_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rsus/WRMC.BSRN/rsus.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Net SW Radiation]
+variable = "rsns"
+alternate_vars = "FSA"
+derived  = "rsds-rsus"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rsns/CERESed4.1/rsns.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rsns/FLUXNET2015/rsns.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rsns/GEWEX.SRB/rsns_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rsns/WRMC.BSRN/rsns.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Upward LW Radiation]
+variable = "rlus"
+alternate_vars = "FIRE"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rlus/CERESed4.1/rlus.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rlus/FLUXNET2015/rlus.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rlus/GEWEX.SRB/rlus_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rlus/WRMC.BSRN/rlus.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Net LW Radiation]
+variable = "rlns"
+derived  = "FLDS-FIRE"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rlns/CERESed4.1/rlns.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rlns/FLUXNET2015/rlns.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rlns/GEWEX.SRB/rlns_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rlns/WRMC.BSRN/rlns.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Net Radiation]
+variable = "rns"
+derived  = "FLDS-FIRE+FSDS-FSR"
+weight = 2
+cmap     = "RdPu"
+
+[CERESed4.1]
+source   = "DATA/rns/CERESed4.1/rns.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rns/FLUXNET2015/rns.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rns/GEWEX.SRB/rns_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rns/WRMC.BSRN/rns.nc"
+weight   = 12
+
+[CLASS]
+source     = "DATA/rns/CLASS/rns.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Ground Heat Flux]
+variable       = "hfdsl"
+alternate_vars = "FGR12"
+cmap           = "Oranges"
+
+[CLASS]
+source     = "DATA/hfdsl/CLASS/hfdsl.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+
+###########################################################################
+
+[h1: Forcings]
+bgcolor = "#EDEDED"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Air Temperature]
+variable = "TSA"
+alternate_vars = "tas"
+weight   = 2
+
+[CRU4.02]
+source   = "DATA/tas/CRU4.02/tas.nc"
+weight   = 25
+
+[FLUXNET2015]
+source   = "DATA/tas/FLUXNET2015/tas.nc"
+weight   = 9
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Diurnal Max Temperature]
+variable = "TREFMXAV"
+alternate_vars = "tasmax"
+weight   = 2
+
+[CRU4.02]
+source   = "DATA/tasmax/CRU4.02/tasmax.nc"
+weight   = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Diurnal Min Temperature]
+variable = "TREFMNAV"
+alternate_vars = "tasmin"
+weight   = 2
+
+[CRU4.02]
+source   = "DATA/tasmin/CRU4.02/tasmin.nc"
+weight   = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Diurnal Temperature Range]
+variable = "dtr"
+weight   = 2
+derived  = "TREFMXAV-TREFMNAV"
+
+[CRU4.02]
+source   = "DATA/dtr/CRU4.02/dtr.nc"
+weight   = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Precipitation]
+variable       = "pr"
+derived        = "RAIN+SNOW" 
+cmap           = "Blues"
+weight         = 2
+mass_weighting = True
+
+[CMAPv1904]
+source     = "DATA/pr/CMAPv1904/pr.nc"
+land       = True
+weight     = 20
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+space_mean = True
+
+[FLUXNET2015]
+source     = "DATA/pr/FLUXNET2015/pr.nc"
+land       = True
+weight     = 9
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+
+[GPCCv2018]
+source     = "DATA/pr/GPCCv2018/pr.nc"
+land       = True
+weight     = 20
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+space_mean = True
+
+[GPCPv2.3]
+source     = "DATA/pr/GPCPv2.3/pr.nc"
+land       = True
+weight     = 20
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+space_mean = True
+
+[CLASS]
+source     = "DATA/pr/CLASS/pr.nc"
+plot_unit  = "mm d-1"
+table_unit = "mm d-1"
+weight     = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Relative Humidity]
+variable       = "RH2M"
+alternate_vars = "rhums,hurs"
+cmap           = "Blues"
+weight         = 3
+mass_weighting = True
+
+[ERA5]
+source     = "DATA/rhums/ERA5/rhums.nc"
+weight     = 10
+
+[CRU4.02]
+source     = "DATA/rhums/CRU4.02/rhums.nc"
+weight     = 10
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Downward SW Radiation]
+variable = "FSDS"
+alternate_vars = "rsds"
+weight   = 2
+
+[CERESed4.1]
+source   = "DATA/rsds/CERESed4.1/rsds.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rsds/FLUXNET2015/rsds.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rsds/GEWEX.SRB/rsds_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rsds/WRMC.BSRN/rsds.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Downward LW Radiation]
+variable = "FLDS"
+alternate_vars = "rlds"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rlds/CERESed4.1/rlds.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rlds/FLUXNET2015/rlds.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rlds/GEWEX.SRB/rlds_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rlds/WRMC.BSRN/rlds.nc"
+weight   = 12
+
diff --git a/src/ILAMB/data/ilamb_nohoff_final_CLM_SP.cfg b/src/ILAMB/data/ilamb_nohoff_final_CLM_SP.cfg
new file mode 100755
index 00000000..a07c4bf8
--- /dev/null
+++ b/src/ILAMB/data/ilamb_nohoff_final_CLM_SP.cfg
@@ -0,0 +1,610 @@
+# This configure file uses observational data which can be obtained by
+# running the following command after exporting ILAMB_ROOT to the
+# appropriate location.
+#
+#   ilamb-fetch --remote_root https://www.ilamb.org/ILAMB-Data
+#
+# Extra ilamb-run options are defined here with the '#!' symbols:
+#
+#! define_regions = DATA/regions/GlobalLand.nc
+#
+[h1: Ecosystem and Carbon Cycle]
+bgcolor = "#ECFFE6"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Gross Primary Productivity]
+variable       = "FPSN"
+alternate_vars = "gpp"
+cmap           = "Greens"
+weight         = 5
+mass_weighting = True
+
+[FLUXNET2015]
+source     = "DATA/gpp/FLUXNET2015/gpp.nc"
+weight     = 9
+table_unit = "g m-2 d-1"
+plot_unit  = "g m-2 d-1"
+
+[FLUXCOM]
+source        = "DATA/gpp/FLUXCOM/gpp.nc"
+weight        = 15
+table_unit    = "Pg yr-1"
+plot_unit     = "g m-2 d-1"
+space_mean    = False
+skip_iav      = True
+relationships = "Evapotranspiration/GLEAMv3.3a","Precipitation/GPCPv2.3","SurfaceDownwardSWRadiation/CERESed4.1","SurfaceNetSWRadiation/CERESed4.1","SurfaceAirTemperature/CRU4.02"
+
+[WECANN]
+source        = "DATA/gpp/WECANN/gpp.nc"
+weight        = 15
+table_unit    = "Pg yr-1"
+plot_unit     = "g m-2 d-1"
+space_mean    = False
+
+[GBAF]
+source        = "DATA/gpp/GBAF/gpp_0.5x0.5.nc"
+weight        = 15
+table_unit    = "Pg yr-1"
+plot_unit     = "g m-2 d-1"
+space_mean    = False
+skip_iav      = True
+relationships = "Evapotranspiration/GLEAMv3.3a","Precipitation/GPCPv2.3","SurfaceDownwardSWRadiation/CERESed4.1","SurfaceNetSWRadiation/CERESed4.1","SurfaceAirTemperature/CRU4.02"
+
+###########################################################################
+
+[h1: Hydrology Cycle]
+bgcolor = "#E6F9FF"
+
+[h2: Evapotranspiration]
+variable       = "et"
+alternate_vars = "evspsbl"
+derived        = "QVEGE+QVEGT+QSOIL"
+cmap           = "Blues"
+weight         = 5
+mass_weighting = True
+
+[GLEAMv3.3a]
+source        = "DATA/evspsbl/GLEAMv3.3a/et.nc"
+weight        = 15
+table_unit    = "mm d-1"
+plot_unit     = "mm d-1"
+relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02"
+
+[MODIS]
+source        = "DATA/evspsbl/MODIS/et_0.5x0.5.nc"
+weight        = 15
+table_unit    = "mm d-1"
+plot_unit     = "mm d-1"
+relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02"
+
+[MOD16A2]
+source        = "DATA/evspsbl/MOD16A2/et.nc"
+weight        = 15
+table_unit    = "mm d-1"
+plot_unit     = "mm d-1"
+relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Evaporative Fraction]
+variable       = "EvapFrac"
+weight         = 5
+mass_weighting = True
+ctype          = "ConfEvapFraction"
+cmap           = "BuGn"
+
+[FLUXCOM]
+hfss_source = "DATA/hfss/FLUXCOM/sh.nc"
+hfls_source = "DATA/hfls/FLUXCOM/le.nc"
+weight      = 9
+skip_rmse   = True
+
+[CLASS]
+variable    = "evapfrac"
+hfss_source = "DATA/hfss/CLASS/hfss.nc"
+hfls_source = "DATA/hfls/CLASS/hfls.nc"
+skip_rmse   = True
+weight      = 25
+
+[WECANN]
+hfss_source = "DATA/hfss/WECANN/hfss.nc"
+hfls_source = "DATA/hfls/WECANN/hfls.nc"
+weight      = 9
+skip_rmse   = True
+
+[GBAF]
+hfss_source     = "DATA/sh/GBAF/sh_0.5x0.5.nc"
+hfls_source     = "DATA/le/GBAF/le_0.5x0.5.nc"
+weight     = 9
+skip_rmse  = True
+skip_iav   = True
+limit_type = "99per"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Latent Heat]
+variable       = "hfls"
+alternate_vars = "le,EFLX_LH_TOT"
+cmap           = "Oranges"
+weight         = 5
+mass_weighting = True
+
+[FLUXNET2015]
+source   = "DATA/hfls/FLUXNET2015/hfls.nc"
+weight   = 3
+
+[FLUXCOM]
+source   = "DATA/hfls/FLUXCOM/le.nc"
+land     = True
+weight   = 9
+skip_iav = True
+
+[DOLCE]
+source   = "DATA/evspsbl/DOLCE/DOLCE.nc"
+weight   = 15
+land     = True
+
+[CLASS]
+source     = "DATA/hfls/CLASS/hfls.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 25
+
+[WECANN]
+source     = "DATA/hfls/WECANN/hfls.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 9
+
+[GBAF]
+source   = "DATA/le/GBAF/le_0.5x0.5.nc"
+land     = True
+weight   = 9
+skip_iav = True
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Runoff]
+variable       = "runoff"
+alternate_vars = "mrro,QRUNOFF"
+weight         = 5
+mass_weighting = True
+
+[Dai]
+ctype          = "ConfRunoff"
+source         = "DATA/mrro/Dai/runoff.nc"
+weight         = 15
+
+[LORA]
+source         = "DATA/mrro/LORA/LORA.nc"
+table_unit     = "mm d-1"
+plot_unit      = "mm d-1"
+weight         = 15
+
+[CLASS]
+source         = "DATA/mrro/CLASS/mrro.nc"
+plot_unit      = "mm d-1"
+table_unit     = "mm d-1"
+weight         = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Sensible Heat]
+variable       = "hfss"
+alternate_vars = "sh,FSH"
+weight         = 2
+mass_weighting = True
+
+[FLUXNET2015]
+source   = "DATA/hfss/FLUXNET2015/hfss.nc"
+weight   = 9
+
+[FLUXCOM]
+source   = "DATA/hfss/FLUXCOM/sh.nc"
+weight   = 15
+skip_iav = True
+
+[CLASS]
+source     = "DATA/hfss/CLASS/hfss.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 25
+
+[WECANN]
+source     = "DATA/hfss/WECANN/hfss.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 15
+
+[GBAF]
+source   = "DATA/sh/GBAF/sh_0.5x0.5.nc"
+weight   = 15
+skip_iav = True
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Terrestrial Water Storage Anomaly]
+variable       = "twsa"
+alternate_vars = "tws,TWS"
+derived        = "RAIN+SNOW-QVEGE-QVEGT-QSOIL-QRUNOFF"
+cmap           = "Blues"
+weight         = 5
+ctype          = "ConfTWSA"
+
+[GRACE]
+source     = "DATA/twsa/GRACE/twsa_0.5x0.5.nc"
+weight     = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Snow Water Equivalent]
+variable       = "H2OSNO"
+alternate_vars = "swe,snw"
+cmap           = "Blues"
+weight         = 5
+ctype          = "ConfSWE"
+
+[CanSISE]
+source     = "DATA/swe/CanSISE/swe.nc"
+weight     = 25
+plot_unit  = "cm"
+table_unit = "cm"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Permafrost]
+variable = "TSOI"
+alternate_vars = "tsl"
+
+[Brown2002]
+ctype    = "ConfPermafrost"
+source   = "DATA/permafrost/Brown2002/Brown2002.nc"
+y0       = 1985.
+yf       = 2005.
+Teps     = 273.15
+dmax     = 3.5
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[Obu2018]
+ctype    = "ConfPermafrost"
+source   = "DATA/permafrost/Obu2018/Obu2018.nc"
+y0       = 2000.
+yf       = 2016.
+Teps     = 273.15
+dmax     = 3.5
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Soil Moisture]
+variable = "SOILWATER_10CM"
+alternate_vars = "mrsos"
+weight   = 3
+cmap     = "Blues"
+
+[WangMao]
+source   = "DATA/mrsos/WangMao/mrsos_olc.nc"
+weight   = 15
+
+###########################################################################
+
+[h1: Radiation and Energy Cycle]
+bgcolor = "#FFECE6"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Albedo]
+variable = "albedo"
+weight   = 1
+ctype    = "ConfAlbedo"
+
+[CERESed4.1]
+source   = "DATA/albedo/CERESed4.1/albedo.nc"
+weight   = 20
+
+[GEWEX.SRB]
+source   = "DATA/albedo/GEWEX.SRB/albedo_0.5x0.5.nc"
+weight   = 20
+
+#[MODIS]
+#source   = "DATA/albedo/MODIS/albedo_0.5x0.5.nc"
+#weight   = 20
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Upward SW Radiation]
+variable = "rsus"
+alternate_vars = "FSR"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rsus/CERESed4.1/rsus.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rsus/FLUXNET2015/rsus.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rsus/GEWEX.SRB/rsus_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rsus/WRMC.BSRN/rsus.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Net SW Radiation]
+variable = "rsns"
+alternate_vars = "FSA"
+derived  = "rsds-rsus"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rsns/CERESed4.1/rsns.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rsns/FLUXNET2015/rsns.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rsns/GEWEX.SRB/rsns_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rsns/WRMC.BSRN/rsns.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Upward LW Radiation]
+variable = "rlus"
+alternate_vars = "FIRE"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rlus/CERESed4.1/rlus.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rlus/FLUXNET2015/rlus.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rlus/GEWEX.SRB/rlus_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rlus/WRMC.BSRN/rlus.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Net LW Radiation]
+variable = "rlns"
+derived  = "FLDS-FIRE"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rlns/CERESed4.1/rlns.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rlns/FLUXNET2015/rlns.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rlns/GEWEX.SRB/rlns_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rlns/WRMC.BSRN/rlns.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Net Radiation]
+variable = "rns"
+derived  = "FLDS-FIRE+FSDS-FSR"
+weight = 2
+cmap     = "RdPu"
+
+[CERESed4.1]
+source   = "DATA/rns/CERESed4.1/rns.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rns/FLUXNET2015/rns.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rns/GEWEX.SRB/rns_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rns/WRMC.BSRN/rns.nc"
+weight   = 12
+
+[CLASS]
+source     = "DATA/rns/CLASS/rns.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+weight     = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Ground Heat Flux]
+variable       = "hfdsl"
+alternate_vars = "FGR12"
+cmap           = "Oranges"
+
+[CLASS]
+source     = "DATA/hfdsl/CLASS/hfdsl.nc"
+plot_unit  = "W m-2"
+table_unit = "W m-2"
+
+###########################################################################
+
+[h1: Forcings]
+bgcolor = "#EDEDED"
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Air Temperature]
+variable = "TSA"
+alternate_vars = "tas"
+weight   = 2
+
+[CRU4.02]
+source   = "DATA/tas/CRU4.02/tas.nc"
+weight   = 25
+
+[FLUXNET2015]
+source   = "DATA/tas/FLUXNET2015/tas.nc"
+weight   = 9
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Diurnal Max Temperature]
+variable = "TREFMXAV"
+alternate_vars = "tasmax"
+weight   = 2
+
+[CRU4.02]
+source   = "DATA/tasmax/CRU4.02/tasmax.nc"
+weight   = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Diurnal Min Temperature]
+variable = "TREFMNAV"
+alternate_vars = "tasmin"
+weight   = 2
+
+[CRU4.02]
+source   = "DATA/tasmin/CRU4.02/tasmin.nc"
+weight   = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Diurnal Temperature Range]
+variable = "dtr"
+weight   = 2
+derived  = "TREFMXAV-TREFMNAV"
+
+[CRU4.02]
+source   = "DATA/dtr/CRU4.02/dtr.nc"
+weight   = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Precipitation]
+variable       = "pr"
+derived        = "RAIN+SNOW" 
+cmap           = "Blues"
+weight         = 2
+mass_weighting = True
+
+[CMAPv1904]
+source     = "DATA/pr/CMAPv1904/pr.nc"
+land       = True
+weight     = 20
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+space_mean = True
+
+[FLUXNET2015]
+source     = "DATA/pr/FLUXNET2015/pr.nc"
+land       = True
+weight     = 9
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+
+[GPCCv2018]
+source     = "DATA/pr/GPCCv2018/pr.nc"
+land       = True
+weight     = 20
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+space_mean = True
+
+[GPCPv2.3]
+source     = "DATA/pr/GPCPv2.3/pr.nc"
+land       = True
+weight     = 20
+table_unit = "mm d-1"
+plot_unit  = "mm d-1"
+space_mean = True
+
+[CLASS]
+source     = "DATA/pr/CLASS/pr.nc"
+plot_unit  = "mm d-1"
+table_unit = "mm d-1"
+weight     = 25
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Relative Humidity]
+variable       = "RH2M"
+alternate_vars = "rhums,hurs"
+cmap           = "Blues"
+weight         = 3
+mass_weighting = True
+
+[ERA5]
+source     = "DATA/rhums/ERA5/rhums.nc"
+weight     = 10
+
+[CRU4.02]
+source     = "DATA/rhums/CRU4.02/rhums.nc"
+weight     = 10
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Downward SW Radiation]
+variable = "FSDS"
+alternate_vars = "rsds"
+weight   = 2
+
+[CERESed4.1]
+source   = "DATA/rsds/CERESed4.1/rsds.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rsds/FLUXNET2015/rsds.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rsds/GEWEX.SRB/rsds_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rsds/WRMC.BSRN/rsds.nc"
+weight   = 12
+
+#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+
+[h2: Surface Downward LW Radiation]
+variable = "FLDS"
+alternate_vars = "rlds"
+weight   = 1
+
+[CERESed4.1]
+source   = "DATA/rlds/CERESed4.1/rlds.nc"
+weight   = 15
+
+[FLUXNET2015]
+source   = "DATA/rlds/FLUXNET2015/rlds.nc"
+weight   = 12
+
+[GEWEX.SRB]
+source   = "DATA/rlds/GEWEX.SRB/rlds_0.5x0.5.nc"
+weight   = 15
+
+[WRMC.BSRN]
+source   = "DATA/rlds/WRMC.BSRN/rlds.nc"
+weight   = 12
+