From 57f0b7cb11a78c6ac7d469171d3974213220fd68 Mon Sep 17 00:00:00 2001 From: chengzhuzhang Date: Fri, 28 Jun 2024 15:53:09 -0500 Subject: [PATCH] update use of ilamb --- zppy/templates/default.ini | 2 +- zppy/templates/ilamb/cmip.cfg | 122 ++++++--- zppy/templates/ilamb/ilamb.cfg | 445 +++++++++++++++++++++++++-------- zppy/templates/ts.bash | 2 +- 4 files changed, 422 insertions(+), 149 deletions(-) diff --git a/zppy/templates/default.ini b/zppy/templates/default.ini index 8ccc8a58..f543c9c0 100644 --- a/zppy/templates/default.ini +++ b/zppy/templates/default.ini @@ -334,7 +334,7 @@ ts_years = string_list(default=list("")) # `years = "1-100",` would plot years 1 to 100 on the graphs. [ilamb] -cfg = string(default="ilamb/cmip.cfg") +cfg = string(default="ilamb/ilamb.cfg") # observational data path, default uses mache's "diagnostics_base_path" ilamb_obs = string(default="") # for land_only run diff --git a/zppy/templates/ilamb/cmip.cfg b/zppy/templates/ilamb/cmip.cfg index c50ba0bb..96ca0a3c 100644 --- a/zppy/templates/ilamb/cmip.cfg +++ b/zppy/templates/ilamb/cmip.cfg @@ -104,7 +104,7 @@ 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" +relationships = "Evapotranspiration/GLEAMv3.3a","Precipitation/GPCPv2.3","SurfaceDownwardSWRadiation/CERESed4.2","SurfaceNetSWRadiation/CERESed4.2","SurfaceAirTemperature/CRU4.02" #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -191,7 +191,7 @@ weight = 5 mass_weighting = True [HWSD] -source = "DATA/soilc/HWSD/soilc_0.5x0.5.nc" +source = "DATA/cSoil/HWSD/soilc_0.5x0.5.nc" weight = 15 table_unit = "Pg" plot_unit = "kg m-2" @@ -199,7 +199,7 @@ space_mean = False skip_rmse = True [NCSCDV22] -source = "DATA/soilc/NCSCDV22/soilc_0.5x0.5.nc" +source = "DATA/cSoil/NCSCDV22/soilc_0.5x0.5.nc" weight = 12 table_unit = "Pg" plot_unit = "kg m-2" @@ -208,14 +208,14 @@ skip_rmse = True [Koven] ctype = "ConfSoilCarbon" -source = "DATA/soilc/NCSCDV22/soilc_0.5x0.5.nc" +source = "DATA/cSoil/NCSCDV22/soilc_0.5x0.5.nc" weight = 15 -soilc_source = "DATA/soilc/NCSCDV22/soilc_0.5x0.5.nc, DATA/soilc/HWSD/soilc_0.5x0.5.nc" +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/soilc/Koven/npp_0.5x0.5.nc" -pet_source = "DATA/soilc/Koven/pet_0.5x0.5.nc" -fracpeat_source = "DATA/soilc/Koven/fracpeat_0.5x0.5.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" ########################################################################### @@ -230,21 +230,21 @@ weight = 5 mass_weighting = True [GLEAMv3.3a] -source = "DATA/et/GLEAMv3.3a/et.nc" +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/et/MODIS/et_0.5x0.5.nc" +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/et/MOD16A2/et.nc" +source = "DATA/evspsbl/MOD16A2/et.nc" weight = 15 table_unit = "mm d-1" plot_unit = "mm d-1" @@ -259,8 +259,14 @@ mass_weighting = True ctype = "ConfEvapFraction" [FLUXCOM] -hfss_source = "DATA/sh/FLUXCOM/sh.nc" -hfls_source = "DATA/le/FLUXCOM/le.nc" +hfss_source = "DATA/hfss/FLUXCOM/sh.nc" +hfls_source = "DATA/hfls/FLUXCOM/le.nc" +weight = 9 +skip_rmse = True + +[CLASS] +hfss_source = "DATA/hfss/CLASS/hfss.nc" +hfls_source = "DATA/hfls/CLASS/hfls.nc" weight = 9 skip_rmse = True @@ -274,20 +280,24 @@ weight = 5 mass_weighting = True [FLUXNET2015] -source = "DATA/le/FLUXNET2015/hfls.nc" +source = "DATA/hfls/FLUXNET2015/hfls.nc" weight = 3 [FLUXCOM] -source = "DATA/le/FLUXCOM/le.nc" +source = "DATA/hfls/FLUXCOM/le.nc" land = True weight = 9 skip_iav = True [DOLCE] -source = "DATA/et/DOLCE/DOLCE.nc" +source = "DATA/evspsbl/DOLCE/DOLCE.nc" weight = 15 land = True +[CLASS] +source = "DATA/hfls/CLASS/hfls.nc" +weight = 15 + #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Runoff] @@ -298,11 +308,17 @@ mass_weighting = True [Dai] ctype = "ConfRunoff" -source = "DATA/runoff/Dai/runoff.nc" +source = "DATA/mrro/Dai/runoff.nc" weight = 15 [LORA] -source = "DATA/runoff/LORA/LORA.nc" +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" table_unit = "mm d-1" plot_unit = "mm d-1" weight = 15 @@ -316,14 +332,18 @@ weight = 2 mass_weighting = True [FLUXNET2015] -source = "DATA/sh/FLUXNET2015/hfss.nc" +source = "DATA/hfss/FLUXNET2015/hfss.nc" weight = 9 [FLUXCOM] -source = "DATA/sh/FLUXCOM/sh.nc" +source = "DATA/hfss/FLUXCOM/sh.nc" weight = 15 skip_iav = True +[CLASS] +source = "DATA/hfss/CLASS/hfss.nc" +weight = 15 + #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Terrestrial Water Storage Anomaly] @@ -343,11 +363,21 @@ weight = 25 [h2: Permafrost] variable = "tsl" -[NSIDC] +[Brown2002] ctype = "ConfPermafrost" -source = "DATA/permafrost/NSIDC/NSIDC_0.5x0.5.nc" -y0 = 1970. -yf = 2000. +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 @@ -363,8 +393,8 @@ variable = "albedo" weight = 1 ctype = "ConfAlbedo" -[CERESed4.1] -source = "DATA/albedo/CERESed4.1/albedo.nc" +[CERESed4.2] +source = "DATA/albedo/CERESed4.2/albedo.nc" weight = 20 [GEWEX.SRB] @@ -377,8 +407,8 @@ weight = 20 variable = "rsus" weight = 1 -[CERESed4.1] -source = "DATA/rsus/CERESed4.1/rsus.nc" +[CERESed4.2] +source = "DATA/rsus/CERESed4.2/rsus.nc" weight = 15 [FLUXNET2015] @@ -400,8 +430,8 @@ variable = "rsns" derived = "rsds-rsus" weight = 1 -[CERESed4.1] -source = "DATA/rsns/CERESed4.1/rsns.nc" +[CERESed4.2] +source = "DATA/rsns/CERESed4.2/rsns.nc" weight = 15 [FLUXNET2015] @@ -422,8 +452,8 @@ weight = 12 variable = "rlus" weight = 1 -[CERESed4.1] -source = "DATA/rlus/CERESed4.1/rlus.nc" +[CERESed4.2] +source = "DATA/rlus/CERESed4.2/rlus.nc" weight = 15 [FLUXNET2015] @@ -445,8 +475,8 @@ variable = "rlns" derived = "rlds-rlus" weight = 1 -[CERESed4.1] -source = "DATA/rlns/CERESed4.1/rlns.nc" +[CERESed4.2] +source = "DATA/rlns/CERESed4.2/rlns.nc" weight = 15 [FLUXNET2015] @@ -468,8 +498,8 @@ variable = "rns" derived = "rlds-rlus+rsds-rsus" weight = 2 -[CERESed4.1] -source = "DATA/rns/CERESed4.1/rns.nc" +[CERESed4.2] +source = "DATA/rns/CERESed4.2/rns.nc" weight = 15 [FLUXNET2015] @@ -484,6 +514,10 @@ weight = 15 source = "DATA/rns/WRMC.BSRN/rns.nc" weight = 12 +[CLASS] +source = "DATA/rns/CLASS/rns.nc" +weight = 15 + ########################################################################### [h1: Forcings] @@ -573,6 +607,14 @@ table_unit = "mm d-1" plot_unit = "mm d-1" space_mean = True +[CLASS] +source = "DATA/pr/CLASS/pr.nc" +land = True +weight = 15 +table_unit = "mm d-1" +plot_unit = "mm d-1" +space_mean = True + #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Surface Relative Humidity] @@ -596,8 +638,8 @@ weight = 10 variable = "rsds" weight = 2 -[CERESed4.1] -source = "DATA/rsds/CERESed4.1/rsds.nc" +[CERESed4.2] +source = "DATA/rsds/CERESed4.2/rsds.nc" weight = 15 [FLUXNET2015] @@ -618,8 +660,8 @@ weight = 12 variable = "rlds" weight = 1 -[CERESed4.1] -source = "DATA/rlds/CERESed4.1/rlds.nc" +[CERESed4.2] +source = "DATA/rlds/CERESed4.2/rlds.nc" weight = 15 [FLUXNET2015] diff --git a/zppy/templates/ilamb/ilamb.cfg b/zppy/templates/ilamb/ilamb.cfg index 770240b0..5311d3f6 100644 --- a/zppy/templates/ilamb/ilamb.cfg +++ b/zppy/templates/ilamb/ilamb.cfg @@ -4,6 +4,10 @@ # # 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" @@ -16,66 +20,91 @@ weight = 5 skip_rmse = True mass_weighting = True -[GEOCARBON] -source = "DATA/biomass/GEOCARBON/biomass_0.5x0.5.nc" +[ESACCI] +source = "DATA/biomass/ESACCI/biomass.nc" +scale_factor = 0.714 # x 0.5 (total mass to carbon) / 0.7 (above only to total) weight = 20 table_unit = "Pg" plot_unit = "kg m-2" space_mean = False -[Tropical] -source = "DATA/biomass/Tropical/biomass_0.5x0.5.nc" -weight = 20 +[GEOCARBON] +source = "DATA/biomass/GEOCARBON/biomass.nc" +scale_factor = 0.714 # x 0.5 (total mass to carbon) / 0.7 (above only to total) +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 -[GlobalCarbon] -source = "DATA/biomass/GLOBAL.CARBON/biomass_0.5x0.5.nc" +[Saatchi2011] +source = "DATA/biomass/Saatchi2011/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 +[Thurner] +source = "DATA/biomass/Thurner/biomass_0.5x0.5.nc" +weight = 16 table_unit = "Pg" plot_unit = "kg m-2" space_mean = False [USForest] -source = "DATA/biomass/US.FOREST/biomass_0.5x0.5.nc" +source = "DATA/biomass/USForest/biomass_0.5x0.5.nc" +scale_factor = 0.714 # x 0.5 (total mass to carbon) / 0.7 (above only to total) weight = 8 table_unit = "Pg" plot_unit = "kg m-2" space_mean = False +[XuSaatchi2021] +source = "DATA/biomass/XuSaatchi2021/XuSaatchi.nc" +weight = 20 +table_unit = "Pg" +plot_unit = "kg m-2" +space_mean = False + #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Burned Area] variable = "burntArea" +alternate_vars = "burntFractionAll" weight = 4 cmap = "OrRd" mass_weighting = True -[GFED4S] -source = "DATA/burntArea/GFED4S/burntArea_0.5x0.5.nc" +[GFED4.1S] +source = "DATA/burntArea/GFED4.1S/burntArea.nc" weight = 20 -relationships = "Precipitation/GPCP2","SurfaceAirTemperature/CRU" +relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02" #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Carbon Dioxide] variable = "co2" -ctype = "ConfCO2" weight = 5 -[NOAA] +[NOAA.Emulated] +ctype = "ConfCO2" 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/HIPPO_AToM/HIPPO_AToM.nc" +model_flux = "ra+rh-gpp" #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -85,20 +114,27 @@ cmap = "Greens" weight = 5 mass_weighting = True -[Fluxnet] -source = "DATA/gpp/FLUXNET/gpp.nc" +[FLUXNET2015] +source = "DATA/gpp/FLUXNET2015/gpp.nc" weight = 9 table_unit = "g m-2 d-1" plot_unit = "g m-2 d-1" -[GBAF] -source = "DATA/gpp/GBAF/gpp_0.5x0.5.nc" +[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/GLEAM","Precipitation/GPCP2","SurfaceDownwardSWRadiation/CERES","SurfaceNetSWRadiation/CERES","SurfaceAirTemperature/CRU" +relationships = "Evapotranspiration/GLEAMv3.3a","Precipitation/GPCPv2.3","SurfaceDownwardSWRadiation/CERESed4.2","SurfaceNetSWRadiation/CERESed4.2","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 #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -111,12 +147,17 @@ mass_weighting = True [AVHRR] source = "DATA/lai/AVHRR/lai_0.5x0.5.nc" weight = 15 -relationships = "Precipitation/GPCP2" +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/GPCP2" +relationships = "Precipitation/GPCPv2.3" skip_iav = True #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -127,7 +168,7 @@ weight = 5 ctype = "ConfNBP" [GCP] -source = "DATA/nbp/GCP/nbp_1959-2012.nc" +source = "DATA/nbp/GCP/nbp_1959-2016.nc" weight = 20 [Hoffman] @@ -139,24 +180,16 @@ skip_taylor = True [h2: Net Ecosystem Exchange] variable = "nee" -derived = "gpp-ra-rh" +derived = "ra+rh-gpp" weight = 5 mass_weighting = True -[Fluxnet] -source = "DATA/nee/FLUXNET/nee.nc" +[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] @@ -165,14 +198,14 @@ derived = "ra+rh" weight = 4 mass_weighting = True -[Fluxnet] -source = "DATA/reco/FLUXNET/reco.nc" +[FLUXNET2015] +source = "DATA/reco/FLUXNET2015/reco.nc" weight = 6 table_unit = "g m-2 d-1" plot_unit = "g m-2 d-1" -[GBAF] -source = "DATA/reco/GBAF/reco_0.5x0.5.nc" +[FLUXCOM] +source = "DATA/reco/FLUXCOM/reco.nc" weight = 4 table_unit = "Pg yr-1" plot_unit = "g m-2 d-1" @@ -182,13 +215,13 @@ skip_iav = True #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Soil Carbon] -variable = "cSoil" -alternate_vars = "soilc" +variable = "cSoilAbove1m" +alternate_vars = "cSoil" weight = 5 mass_weighting = True [HWSD] -source = "DATA/soilc/HWSD/soilc_0.5x0.5.nc" +source = "DATA/cSoil/HWSD/soilc_0.5x0.5.nc" weight = 15 table_unit = "Pg" plot_unit = "kg m-2" @@ -196,13 +229,51 @@ space_mean = False skip_rmse = True [NCSCDV22] -source = "DATA/soilc/NCSCDV22/soilc_0.5x0.5.nc" +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 = "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 + +#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +[h2: Methane] +variable = "FCH4" +alternate_vars = "ch4" +cmap = "Greens" +mass_weighting = True + +[FluxnetANN] +source = "DATA/ch4/FluxnetANN/FCH4_F_ANN_monthly_wetland_tier1.nc" +table_unit = "g m-2 d-1" +plot_unit = "g m-2 d-1" + ########################################################################### [h1: Hydrology Cycle] @@ -215,19 +286,26 @@ cmap = "Blues" weight = 5 mass_weighting = True -[GLEAM] -source = "DATA/et/GLEAM/et_0.5x0.5.nc" +[GLEAMv3.3a] +source = "DATA/evspsbl/GLEAMv3.3a/et.nc" weight = 15 table_unit = "mm d-1" plot_unit = "mm d-1" -relationships = "Precipitation/GPCP2","SurfaceAirTemperature/CRU" +relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02" [MODIS] -source = "DATA/et/MODIS/et_0.5x0.5.nc" +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/GPCP2","SurfaceAirTemperature/CRU" +relationships = "Precipitation/GPCPv2.3","SurfaceAirTemperature/CRU4.02" #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -236,13 +314,25 @@ variable = "EvapFrac" weight = 5 mass_weighting = True ctype = "ConfEvapFraction" - -[GBAF] -source = "DATA/EvapFrac/GBAF/EvapFrac_0.5x0.5.nc" -weight = 9 -skip_rmse = True -skip_iav = True -limit_type = "99per" +cmap = "BuGn" + +[FLUXCOM] +hfss_source = "DATA/hfss/FLUXCOM/sh.nc" +hfls_source = "DATA/hfls/FLUXCOM/le.nc" +weight = 9 +skip_rmse = True + +[CLASS] +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 #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -253,28 +343,58 @@ cmap = "Oranges" weight = 5 mass_weighting = True -[Fluxnet] -source = "DATA/le/FLUXNET/le.nc" +[FLUXNET2015] +source = "DATA/hfls/FLUXNET2015/hfls.nc" weight = 3 -[GBAF] -source = "DATA/le/GBAF/le_0.5x0.5.nc" +[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 + #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Runoff] variable = "runoff" alternate_vars = "mrro" weight = 5 +mass_weighting = True [Dai] ctype = "ConfRunoff" -source = "DATA/runoff/Dai/runoff.nc" +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] @@ -283,15 +403,27 @@ alternate_vars = "sh" weight = 2 mass_weighting = True -[Fluxnet] -source = "DATA/sh/FLUXNET/sh.nc" +[FLUXNET2015] +source = "DATA/hfss/FLUXNET2015/hfss.nc" weight = 9 -[GBAF] -source = "DATA/sh/GBAF/sh_0.5x0.5.nc" +[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 + #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Terrestrial Water Storage Anomaly] @@ -326,14 +458,35 @@ table_unit = "cm" [h2: Permafrost] variable = "tsl" -[NSIDC] +[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/NSIDC/NSIDC_0.5x0.5.nc" -y0 = 1970. -yf = 2000. +source = "DATA/permafrost/Obu2018/Obu2018.nc" +y0 = 2000. +yf = 2016. Teps = 273.15 dmax = 3.5 +#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +[h2: Surface Soil Moisture] +variable = "mrsos" +weight = 3 +cmap = "Blues" + +[WangMao] +source = "DATA/mrsos/WangMao/mrsos_olc.nc" +weight = 15 + ########################################################################### [h1: Radiation and Energy Cycle] @@ -346,28 +499,28 @@ variable = "albedo" weight = 1 ctype = "ConfAlbedo" -[CERES] -source = "DATA/albedo/CERES/albedo_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/albedo/CERESed4.2/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" weight = 1 -[CERES] -source = "DATA/rsus/CERES/rsus_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/rsus/CERESed4.2/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 @@ -383,10 +536,14 @@ variable = "rsns" derived = "rsds-rsus" weight = 1 -[CERES] -source = "DATA/rsns/CERES/rsns_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/rsns/CERESed4.2/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 @@ -401,10 +558,14 @@ weight = 12 variable = "rlus" weight = 1 -[CERES] -source = "DATA/rlus/CERES/rlus_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/rlus/CERESed4.2/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 @@ -420,10 +581,14 @@ variable = "rlns" derived = "rlds-rlus" weight = 1 -[CERES] -source = "DATA/rlns/CERES/rlns_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/rlns/CERESed4.2/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 @@ -437,14 +602,15 @@ weight = 12 [h2: Surface Net Radiation] variable = "rns" derived = "rlds-rlus+rsds-rsus" -weight = 2 +weight = 2 +cmap = "RdPu" -[CERES] -source = "DATA/rns/CERES/rns_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/rns/CERESed4.2/rns.nc" weight = 15 -[Fluxnet] -source = "DATA/rns/FLUXNET/rns.nc" +[FLUXNET2015] +source = "DATA/rns/FLUXNET2015/rns.nc" weight = 12 [GEWEX.SRB] @@ -455,6 +621,24 @@ weight = 15 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] @@ -466,53 +650,92 @@ bgcolor = "#EDEDED" variable = "tas" weight = 2 -[CRU] -source = "DATA/tas/CRU/tas_0.5x0.5.nc" +[CRU4.02] +source = "DATA/tas/CRU4.02/tas.nc" weight = 25 -[Fluxnet] -source = "DATA/tas/FLUXNET/tas.nc" +[FLUXNET2015] +source = "DATA/tas/FLUXNET2015/tas.nc" weight = 9 #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ +[h2: Diurnal Max Temperature] +variable = "tasmax" +weight = 2 + +[CRU4.02] +source = "DATA/tasmax/CRU4.02/tasmax.nc" +weight = 25 + +#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +[h2: Diurnal Min Temperature] +variable = "tasmin" +weight = 2 + +[CRU4.02] +source = "DATA/tasmin/CRU4.02/tasmin.nc" +weight = 25 + +#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +[h2: Diurnal Temperature Range] +variable = "dtr" +weight = 2 +derived = "tasmax-tasmin" + +[CRU4.02] +source = "DATA/dtr/CRU4.02/dtr.nc" +weight = 25 + +#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + [h2: Precipitation] variable = "pr" cmap = "Blues" weight = 2 mass_weighting = True -[CMAP] -source = "DATA/pr/CMAP/pr_0.5x0.5.nc" +[CMAPv1904] +source = "DATA/pr/CMAPv1904/pr.nc" land = True weight = 20 table_unit = "mm d-1" plot_unit = "mm d-1" space_mean = True -[Fluxnet] -source = "DATA/pr/FLUXNET/pr.nc" +[FLUXNET2015] +source = "DATA/pr/FLUXNET2015/pr.nc" land = True weight = 9 table_unit = "mm d-1" plot_unit = "mm d-1" -[GPCC] -source = "DATA/pr/GPCC/pr_0.5x0.5.nc" +[GPCCv2018] +source = "DATA/pr/GPCCv2018/pr.nc" land = True weight = 20 table_unit = "mm d-1" plot_unit = "mm d-1" space_mean = True -[GPCP2] -source = "DATA/pr/GPCP2/pr_0.5x0.5.nc" +[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" +land = True +weight = 15 +table_unit = "mm d-1" +plot_unit = "mm d-1" +space_mean = True + #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [h2: Surface Relative Humidity] @@ -522,8 +745,12 @@ cmap = "Blues" weight = 3 mass_weighting = True -[ERA] -source = "DATA/rhums/ERA/rhums_0.5x0.5.nc" +[ERA5] +source = "DATA/rhums/ERA5/rhums.nc" +weight = 10 + +[CRU4.02] +source = "DATA/rhums/CRU4.02/rhums.nc" weight = 10 #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ @@ -532,12 +759,12 @@ weight = 10 variable = "rsds" weight = 2 -[CERES] -source = "DATA/rsds/CERES/rsds_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/rsds/CERESed4.2/rsds.nc" weight = 15 -[Fluxnet] -source = "DATA/rsds/FLUXNET/rsds.nc" +[FLUXNET2015] +source = "DATA/rsds/FLUXNET2015/rsds.nc" weight = 12 [GEWEX.SRB] @@ -554,10 +781,14 @@ weight = 12 variable = "rlds" weight = 1 -[CERES] -source = "DATA/rlds/CERES/rlds_0.5x0.5.nc" +[CERESed4.2] +source = "DATA/rlds/CERESed4.2/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 diff --git a/zppy/templates/ts.bash b/zppy/templates/ts.bash index 7141e60a..91dceccf 100644 --- a/zppy/templates/ts.bash +++ b/zppy/templates/ts.bash @@ -159,7 +159,7 @@ EOF {% endif -%} {% if input_files.split(".")[0] == 'cam' or input_files.split(".")[0] == 'eam' -%} --var-list \ - 'pr, tas, rsds, rlds, rsus' \ + 'pfull, phalf, tas, ts' \ --realm \ atm \ {% endif -%}