Get init_interp working correctly with hillslope initial conditions files #1806
Labels
enhancement
new capability or improved behavior of existing capability
priority: low
Background task that doesn't need to be done right away.
I'm pretty sure that init_interp won't currently work correctly with initial conditions files from a hillslope run. (I'm also not positive it's working ideally for interpolation from a non-hillslope run to a hillslope run, but it seems likely that it's at least doing something reasonable in that situation.) One need – maybe the main one – is to label each column with appropriate metadata so that, in the interpolation, we can copy data from the source lowland column to the destination lowland column, etc. This can be done with discrete labels giving the type of each column and/or a continuous (real-valued) variable that lets us find the closest match. Currently, all vegetated columns are given the same type index and init_interp doesn't have a way to handle multiple columns in the vegetated landunit. (In contrast, for some other landunits, like crop and glacier, we have mechanisms to ensure that the appropriate source column from a landunit is associated with a given destination column.)
A good first-order check of whether this is working correctly is whether interpolating an initial conditions file onto an identical configuration will result in a replica of the original initial conditions file (this is what the LII test checks). Currently I believe this would not be the case. But eventually - at least if hillslope becomes default - we also want this to be doing something reasonable when interpolating from one resolution / hillslope configuration to another.
@swensosc
The text was updated successfully, but these errors were encountered: