fix p gradient sign

python-hydro · Oct 21, 2024 · 0bde429 · 0bde429
1 parent 26ffe6f
commit 0bde429
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/pyro/compressible/simulation.py b/pyro/compressible/simulation.py
@@ -347,9 +347,9 @@ def evolve(self):
         # fluxes
 
         if myg.coord_type == 1:
-            xmom.v()[:, :] += self.dt * (qx.ip(1, n=self.ivars.ip) -
+            xmom.v()[:, :] -= self.dt * (qx.ip(1, n=self.ivars.ip) -
                                          qx.v(n=self.ivars.ip)) / myg.Lx.v()
-            ymom.v()[:, :] += self.dt * (qy.jp(1, n=self.ivars.ip) -
+            ymom.v()[:, :] -= self.dt * (qy.jp(1, n=self.ivars.ip) -
                                          qy.v(n=self.ivars.ip)) / myg.Ly.v()
 
         # now the external sources (including gravity).  We are going