pg_utils.pg_model.forcing.Le_p_cg

pg_utils.pg_model.forcing.Le_p_cg = -sqrt(2)*I*(B_{+, z}^e(s, \phi, t) - B_{-, z}^e(s, \phi, t))*B_{z}^e(s, \phi, t)/2 + sqrt(2)*(B_{+}^e(s, \phi, t) + B_{-}^e(s, \phi, t))*Derivative(-sqrt(2)*I*(B_{+}^e(s, \phi, t) - B_{-}^e(s, \phi, t))/2, s)/2 - I*(B_{+}^e(s, \phi, t) - B_{-}^e(s, \phi, t))*(B_{+}^e(s, \phi, t) + B_{-}^e(s, \phi, t))/(2*s) - sqrt(2)*I*(B_{+}^e(s, \phi, t) - B_{-}^e(s, \phi, t))*Derivative(-sqrt(2)*I*(B_{+}^e(s, \phi, t) - B_{-}^e(s, \phi, t))/2, \phi)/(2*s)

Expression: the azimuthal Lorentz force in the equatorial plane in conjugate vars