Timesteps with multiple airflow conditions (temps) were creating flattened XC shapes (ceiling) during phreatic flow conditions with large erosional rates. This was resulting from change in position of the ceiling without recalculation of flow (causing ceiling to not get eroded). I have fixed this by adding flow calculations in for each temperature sub-step. However, there may be a more efficient way to refactor this such that flow is only calculated once and all erosion is applied at the end of the timestep.
Timesteps with multiple airflow conditions (temps) were creating flattened XC shapes (ceiling) during phreatic flow conditions with large erosional rates. This was resulting from change in position of the ceiling without recalculation of flow (causing ceiling to not get eroded). I have fixed this by adding flow calculations in for each temperature sub-step. However, there may be a more efficient way to refactor this such that flow is only calculated once and all erosion is applied at the end of the timestep.