Updated Init.cpp to incorporate crit_nucleus simulations

AMReX-FHD · Jan 24, 2025 · 5128d99 · 5128d99
1 parent be88b5c
commit 5128d99
Showing 1 changed file with 3 additions and 4 deletions.
diff --git a/exec/multispec/Init.cpp b/exec/multispec/Init.cpp
@@ -44,8 +44,7 @@ void InitRhoUmac(std::array< MultiFab, AMREX_SPACEDIM >& umac,
               c=c_init_1(:) inside, c=c_init_2(:) outside
               can be discontinous or smooth depending on smoothing_width
             */
-            Real rad = L[0] / 4.;
-	    rad = radius_cyl;
+            Real rad = radius_cyl;
 
             amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
             {
@@ -61,11 +60,11 @@ void InitRhoUmac(std::array< MultiFab, AMREX_SPACEDIM >& umac,
                     // discontinuous interface
                     if (r < rad) {
                         for (int n=0; n<nspecies; ++n) {
-                            c(i,j,k,n) = c_init_1[n];
+                            c(i,j,k,n) = c_init_2[n];
                         }
                     } else {
                         for (int n=0; n<nspecies; ++n) {
-                            c(i,j,k,n) = c_init_2[n];
+                            c(i,j,k,n) = c_init_1[n];
                         }
                     }