Skip to content

Commit

Permalink
use plot_deltaY_dir instead of project_dir
Browse files Browse the repository at this point in the history
  • Loading branch information
ajnonaka committed Dec 5, 2024
1 parent dd17133 commit 9f1ccd2
Show file tree
Hide file tree
Showing 3 changed files with 8 additions and 8 deletions.
6 changes: 3 additions & 3 deletions src_common/common_functions.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -240,7 +240,7 @@ int common::plot_means;
int common::plot_vars;
int common::plot_covars;
int common::plot_cross;
int common::plot_deltaY;
int common::plot_deltaY_dir;
int common::particle_motion;

AMREX_GPU_MANAGED amrex::Real common::turb_a;
Expand Down Expand Up @@ -609,7 +609,7 @@ void InitializeCommonNamespace() {
plot_vars = 0;
plot_covars = 0;
plot_cross = 0;
plot_deltaY = 0;
plot_deltaY_dir = -1;
particle_motion = 0;

// turblent forcing parameters
Expand Down Expand Up @@ -1137,7 +1137,7 @@ void InitializeCommonNamespace() {
pp.query("plot_vars",plot_vars);
pp.query("plot_covars",plot_covars);
pp.query("plot_cross",plot_cross);
pp.query("plot_deltaY",plot_deltaY);
pp.query("plot_deltaY_dir",plot_deltaY_dir);
pp.query("particle_motion",particle_motion);
pp.query("turb_a",turb_a);
pp.query("turb_b",turb_b);
Expand Down
2 changes: 1 addition & 1 deletion src_common/common_namespace.H
Original file line number Diff line number Diff line change
Expand Up @@ -314,7 +314,7 @@ namespace common {
extern int plot_vars;
extern int plot_covars;
extern int plot_cross;
extern int plot_deltaY;
extern int plot_deltaY_dir;
extern int particle_motion;

// parameters for turbulent forcing example
Expand Down
8 changes: 4 additions & 4 deletions src_compressible_stag/writePlotFileStag.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -88,7 +88,7 @@ void WritePlotFileStag(int step,
if (nspec_surfcov>0) nplot += nspec_surfcov*6;
}

if (plot_deltaY == 1) {
if (plot_deltaY_dir != -1) {
nplot += nspecies;
}

Expand Down Expand Up @@ -261,10 +261,10 @@ void WritePlotFileStag(int step,
}
}

if (plot_deltaY == 1) {
if (plot_deltaY_dir != -1) {
MultiFab Ybar(ba, dmap, nspecies, 0);
// Yk is component 6: in prim
WriteHorizontalAverageToMF(prim,Ybar,project_dir,6,nspecies,0);
WriteHorizontalAverageToMF(prim,Ybar,plot_deltaY_dir,6,nspecies,0);
Ybar.mult(-1.);
amrex::MultiFab::Add(Ybar,prim,6,0,nspecies,0);
amrex::MultiFab::Copy(plotfile,Ybar,0,cnt,nspecies,0);
Expand Down Expand Up @@ -459,7 +459,7 @@ void WritePlotFileStag(int step,

}

if (plot_deltaY == 1) {
if (plot_deltaY_dir != 1) {
x = "deltaYk_";
for (i=0; i<nspecies; i++) {
varNames[cnt] = x;
Expand Down

0 comments on commit 9f1ccd2

Please sign in to comment.