Skip to content

Commit

Permalink
bugs in 1D SF
Browse files Browse the repository at this point in the history
  • Loading branch information
ajnonaka committed Jan 30, 2025
1 parent 00f1f9d commit cffcc5c
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions src_compressible_stag/main_driver.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -849,7 +849,7 @@ void main_driver(const char* argv)
}
}
}

if (n_ads_spec > 0 && (n_cells[0] != 1 || n_cells[1] != 1) ) {

MultiFab Flattened; // flattened multifab defined below
Expand Down Expand Up @@ -899,8 +899,8 @@ void main_driver(const char* argv)
structFactConsPencil.resize(n_cells[1]*n_cells[2]);

for (int i = 0; i < n_cells[1]*n_cells[2]; ++i) {
structFactPrimPencil[i].define(ba_flat,dmap_flat,prim_var_names,var_scaling_prim);
structFactConsPencil[i].define(ba_flat,dmap_flat,cons_var_names,var_scaling_cons);
structFactPrimPencil[i].define(ba_pencil,dmap_pencil,prim_var_names,var_scaling_prim);
structFactConsPencil[i].define(ba_pencil,dmap_pencil,cons_var_names,var_scaling_cons);
}

}
Expand Down Expand Up @@ -1444,14 +1444,14 @@ void main_driver(const char* argv)
{
MultiFab pencil;

ExtractSlice(structFactPrimMF, pencil, i/n_cells[1], i%n_cells[1], 0, structVarsPrim);
ExtractXPencil(structFactPrimMF, pencil, i/n_cells[1], i%n_cells[1], 0, structVarsPrim);
structFactPrimPencil[i].FortStructure(pencil);
}

{
MultiFab pencil;

ExtractSlice(structFactConsMF, pencil, i/n_cells[1], i%n_cells[1], 0, structVarsCons);
ExtractXPencil(structFactConsMF, pencil, i/n_cells[1], i%n_cells[1], 0, structVarsCons);
structFactConsPencil[i].FortStructure(pencil);
}

Expand Down

0 comments on commit cffcc5c

Please sign in to comment.