restore energy computation in integrals
add mean flow forcing skeleton