more fixes

python-hydro · Jan 24, 2025 · 0fad800 · 0fad800
1 parent 23e859c
commit 0fad800
Show file tree

Hide file tree

Showing 5 changed files with 5 additions and 2 deletions.
diff --git a/pyro/compressible_fv4/_defaults b/pyro/compressible_fv4/_defaults
@@ -23,6 +23,7 @@ grav = 0.0                ; gravitational acceleration (in y-direction)
 
 riemann = CGF
 
+small_dens = -1.e200      ; minimum allowed density
 
 [sponge]
 do_sponge = 0             ; do we include a sponge source term

diff --git a/pyro/compressible_fv4/simulation.py b/pyro/compressible_fv4/simulation.py
@@ -19,7 +19,7 @@ def substep(self, myd):
         compute the advective source term for the given state
         """
 
-        self.clean_state(myd)
+        self.clean_state(myd.data)
 
         myg = myd.grid
 

diff --git a/pyro/compressible_rk/_defaults b/pyro/compressible_rk/_defaults
@@ -25,6 +25,7 @@ riemann = HLLC            ; HLLC or CGF
 
 well_balanced = 0         ; use a well-balanced scheme to keep the model in hydrostatic equilibrium
 
+small_dens = -1.e200      ; minimum allowed density
 
 [sponge]
 do_sponge = 0             ; do we include a sponge source term

diff --git a/pyro/compressible_rk/simulation.py b/pyro/compressible_rk/simulation.py
@@ -15,7 +15,7 @@ def substep(self, myd):
         conservative state defined as part of myd
         """
 
-        self.clean_state(myd)
+        self.clean_state(myd.data)
 
         myg = myd.grid
 

diff --git a/pyro/compressible_sdc/_defaults b/pyro/compressible_sdc/_defaults
@@ -23,6 +23,7 @@ grav = 0.0                ; gravitational acceleration (in y-direction)
 
 riemann = CGF
 
+small_dens = -1.e200      ; minimum allowed density
 
 [sponge]
 do_sponge = 0             ; do we include a sponge source term