• Christopher Subich's avatar
    Gravity_current: compute v^2 KE if f != 0 · ae19aedf
    Christopher Subich authored
    This change allows the gravity_current case to compute the v^2 component
    of the kinetic energy term even if NY=1 when rot_f != 0; in that case
    the solution can still develop a 2D geostrophic balance where v-velocity
    balances a pressure gradient.
    
    Additionally, the respective component sums are explicitly initialized to
    0.
    ae19aedf
gravity_current.cpp 22.8 KB