Open JeffreyEarly opened 5 months ago
Seems like this should maybe be implement as part of the forcing modules. The nonlinear flux could 'own' the various forcing components. Quite possibility we want the NonlinearFlux to know how to compute the spatial and spectral forcing effects, as well as the inertial flux.
This would allow us to dump all operations associated with a particular nonlinear flux when it has been removed.
Still need to solve the problem that u_g gets redefined, so what do we do?