Open wwieder opened 3 years ago
I believe what's happening is that by moving tauMod outside of the sqrt, it's applied consistently for all points, rather than having a weighted effect based on ANPP.
To test if the change in the equation is truly beneficial may be hard, since we'll likely need to find a new value for tauMod. My parameterization script could be modified to help find that ~value.
@piersond found a better parameterization for tauMod at Reynolds Creek with the following function
SQRT(ANPP)/tauModDenom
.https://github.com/wwieder/biogeochem_testbed_1.1/blob/52d460fd39720c3f7dc7b47265f20f4c6149796d/SOURCE_CODE/mimics_cycle_CN.f90#L1602
This is something we can consider trying in global parameterizations?