fix errors not filling mixing coeffs

hklion · Dec 19, 2023 · 9927344 · 9927344
1 parent 50a5e13
commit 9927344
Show file tree

Hide file tree

Showing 4 changed files with 8 additions and 6 deletions.
diff --git a/Exec/Advection/inputs b/Exec/Advection/inputs
@@ -2,9 +2,9 @@
 max_step = 10
 stop_time = 30000.0
 
-amrex.fpe_trap_invalid = 0
+amrex.fpe_trap_invalid = 1
 
-fabarray.mfiter_tile_size = 16 16 1024
+fabarray.mfiter_tile_size = 1024 1024 1024
 
 # PROBLEM SIZE & GEOMETRY
 geometry.prob_lo     =      0.     0.    -150.

diff --git a/Exec/DoublyPeriodic/prob.cpp b/Exec/DoublyPeriodic/prob.cpp
@@ -213,7 +213,9 @@ init_custom_vmix(const Geometry& /*geom*/, MultiFab& mf_Akv, MultiFab& mf_Akt,
       [=] AMREX_GPU_DEVICE (int i, int j, int k)
       {
         Akv(i,j,k) = 2.0e-03+8.0e-03*std::exp(z_w(i,j,k)/150.0);
-        Akt(i,j,k) = 1.0e-6_rt;
+        Akt(i,j,k,Temp_comp) = 1.0e-6_rt;
+        Akt(i,j,k,Salt_comp) = 1.0e-6_rt;
+        Akt(i,j,k,Scalar_comp) = 0.0_rt;
       });
     }
 }

diff --git a/Exec/Seamount/prob.cpp b/Exec/Seamount/prob.cpp
@@ -186,7 +186,7 @@ init_custom_hmix(const Geometry& /*geom*/, MultiFab& mf_visc2_p, MultiFab& mf_vi
       Box bx = mfi.tilebox();
       bx.grow(IntVect(NGROW,NGROW,0));
 
-      int ncomp = 2; // temperature and salt
+      int ncomp = mf_diff2.nComp(); // temperature and salt and scalar
       Gpu::streamSynchronize();
 
       amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k)

diff --git a/Tests/test_files/Advection/Advection.i b/Tests/test_files/Advection/Advection.i
@@ -2,9 +2,9 @@
 max_step = 10
 stop_time = 30000.0
 
-amrex.fpe_trap_invalid = 0
+amrex.fpe_trap_invalid = 1
 
-fabarray.mfiter_tile_size = 16 16 1024
+fabarray.mfiter_tile_size = 1024 1024 1024
 
 # PROBLEM SIZE & GEOMETRY
 geometry.prob_lo     =      0.     0.    -150.