Commit 2a7736a9 authored by Christopher Subich's avatar Christopher Subich
Browse files

wave_reader, tracer file output

Previously for the wave_reader case file, "tracer.0" and
"tracer_reader.m" files were being output regardless of whether the
configuration specified that there was a passive tracer.  This arose
because of a scoping error during tracer initalization; tracer.0
contained a copy of rho.0.

This has been fixed in this commit.  No numerical results were in error
because of this bug.
parent 1d616e81
......@@ -397,9 +397,9 @@ class userControl : public BaseCase {
read_array(rho,tracer_filename.c_str(),Nx,Ny,Nz);
break;
}
write_array(rho,"tracer",plotnum);
write_reader(rho,"tracer",true);
}
write_array(rho,"tracer",plotnum);
write_reader(rho,"tracer",true);
}
void forcing(double t, DTArray & u, DTArray & u_f,
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment