Data and code associated with "Land surface conductance linked to precipitation: co-evolution of vegetation and climate in Earth system models"

This record contains four sets of files related to the paper "Land surface conductance linked to precipitation: co-evolution of vegetation and climate in Earth system models";

1. Two netCDF files containing high resolution global values of g1m and g1b, calculated using the MAP method (see paper for details).
2. Post-processed netCDF files based on output from the NCAR Community Land Model 5 (CLM), for the 5 variables analysed in the paper (see table 1 in paper).
3. Code modifications for CLM that implement the g1 modifications used in the simulations described in the paper.
4. Two surface data files ("surfdat" files) for CLM that include the 2D g1 values needed for the "BB MAP" and "MED MAP" simulations (see Table 1 in the paper). These are the same g1 values in item 1 above but interpolated to our model's resolution and named so that code modifications in item 3 can access them.

These four sets of files are described further in the sections below.

NOTE ON CESM SET UP: CESM2.1.3 was used to access CLM. The compset "FHIST" was used with the grid "f09_f09_mg17".

1. Global g1m and g1b files

g1m_MAP-method.nc and g1b_MAP-method.nc contain g1m and g1b, respectively, as shown and discussed in the paper (Figure 1). Note this data is at ~5 km resolution and NULL values over the ocean have been filled using the poisson_grid_fill function from the NCL programming language. However, open ocean values should not be trusted (ocean values are provided in case coastal cells that are partially land but registered as ocean require values).

2. CLM output netCDF files

These files contain post-processed output from the CLM. The naming convention for netCDF files is: 

<variable_name>.<simulation_name>.<time_period>.<frequency>.nc

<variable_name> can be:

<simulation_name> can be:

<time_period> is either "1981-2010" or "2071-2100".

<frequency> is "monthly", "yearmax" or "h3yearmax" and refers to the frequency of the variable in the file. Only 3 of the 5 variables have "yearmax" ("h3yearmax" in the case of TSA) since these were the only variables with corresponding daily observations. There are eight netCDF files for each simulation.

3. CLM code modifications

Two zip files containing modified CLM code (in the form of FORTRAN files .f90 files) are included in this record. These zip files must be extracted into the SourceMods/src.clm/ subdirectory of the CESM simulation prior to building the model. These files work for CLM5, specifically the version of CLM5 released with CESM2.1.3. 

CLM5_g1MAP_PFT.zip contains files for the "Default" and "MAP PFT" simulations (see table above).

CLM5_g1MAP_static2D.zip contains files for the "MAP" simulations (see table above).

4. CLM surface data files

surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr1850_c190214_g1m-MAP_06082021.nc

The above file is the input surface data file for CLM used in the "MED MAP" simulation. It contains a new 2D variable named "g1" that contains g1m calculated from the MAP method described in the manuscript. The code modifications in CLM5_g1MAP_static2D.zip (described in the section above) make CLM read these g1 values instead of using the default PFT-specific g1 values.

surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr1850_c190214_g1b-MAP_07072021.nc

The above file is the same except that it contains g1b instead of g1m.