New: (1) The hydrological parameters: mcs, mcf, mcw (as well as the _lin coefficients derived from mcs/mcf/mcw), are now SOC-dependent, therefore varies for each pixel. No vertical variation at the moment.
To be specific:
(i) read SOC map (in unit of gC/m2) in hydrol_initialize
(ii) change mcs as a linear combination of mineral soil and organic soil, the same as how we treat soil porosity
(iii) use the Van Genuchten relationship to calculate mcf and mcw for each USDA soil texture and for each pixel.
To activate this option, set "use_refSOC_hydrol=y" in run.def. (for now only works for the USDA soil texture)
(2) nroot (the weighting factor to calculate water stress) is now dynamic for PFT9 (boreal deciduous needleleaf) in order to optimize its usage of water: soil layers with higher relative moisture have larger nroot (similar to Beer et al., 2007).
To activate this option, set "dyn_nroot_larix=y" in run.def
Done by Dan
