diff --git a/demos/quasigeostrophy_1layer/qg_1layer_wave.py.rst b/demos/quasigeostrophy_1layer/qg_1layer_wave.py.rst
index d0d90c0383..25138f82ff 100644
--- a/demos/quasigeostrophy_1layer/qg_1layer_wave.py.rst
+++ b/demos/quasigeostrophy_1layer/qg_1layer_wave.py.rst
@@ -147,8 +147,8 @@ along with the physical parameters of the model. ::
 Next, we define the variational problems.  First the elliptic problem
 for the stream function. ::
 
-  psi = TrialFunction(Vcg)  # Test function
-  phi = TestFunction(Vcg)   # Trial function
+  psi = TrialFunction(Vcg)
+  phi = TestFunction(Vcg)
 
   # Build the weak form for the inversion
   Apsi = (inner(grad(psi),grad(phi)) + F*psi*phi)*dx