Boussinesq equation =================== .. only:: solution Reference solution ------------------ .. literalinclude:: nse_cylinder.py :start-after: # Begin code