The function calc_adaptive_dtau() is not only incredibly slow, it is also not completely finished in the sense that it does not include all the robustness control capabilities included in GEMS. As I recall, the most expensive component of the adaptive dtau is calculating the mixture dynamic viscosity, which is not previously computed (it's only calculated at the face average). As such, the adaptive dtau should be completely fleshed out and some means of making the calculation less expensive should be sussed out.
The function
calc_adaptive_dtau()is not only incredibly slow, it is also not completely finished in the sense that it does not include all the robustness control capabilities included in GEMS. As I recall, the most expensive component of the adaptive dtau is calculating the mixture dynamic viscosity, which is not previously computed (it's only calculated at the face average). As such, the adaptive dtau should be completely fleshed out and some means of making the calculation less expensive should be sussed out.