Removing useless terms from the equations in the 2 layer QG example

Author: jodemaey

examples/atmospheric/two_layer_quasi_geostrophic.py

@@ -121,9 +121,6 @@
 psi2_advec = Jacobian(psi1, psi2, atmospheric_basis.coordinate_system, sign=-1, prefactors=(alphap_1, alphap_1))
 psi1_equation.add_rhs_terms(psi2_advec)
 
-psi1_advec = Jacobian(psi1, psi1, atmospheric_basis.coordinate_system, sign=-1, prefactors=(dalpha_1, dalpha_1))
-psi1_equation.add_rhs_terms(psi1_advec)
-
 # adding the beta term
 beta_term = OperatorTerm(psi1, D, x, prefactor=beta_nondim, sign=-1)
 psi1_equation.add_rhs_term(beta_term)
@@ -144,9 +141,6 @@
 advection_term = vorticity_advection(psi2, psi2, atmospheric_basis.coordinate_system, sign=-1)
 psi2_equation.add_rhs_terms(advection_term)
 
-psi2_advec = Jacobian(psi2, psi2, atmospheric_basis.coordinate_system, sign=1, prefactors=(alphap_2, alphap_2))
-psi2_equation.add_rhs_terms(psi2_advec)
-
 psi1_advec = Jacobian(psi2, psi1, atmospheric_basis.coordinate_system, sign=-1, prefactors=(alphap_2, alphap_2))
 psi2_equation.add_rhs_terms(psi1_advec)