Commit 74357af5 authored by Geoffroy Lesur's avatar Geoffroy Lesur 👻
Browse files

fix the typo regarding shearing box BCs

parent 25a32295
......@@ -107,7 +107,9 @@ void Setup::InitFlow(DataBlock &data) {
#else
d.Vc(RHO,k,j,i) = 1.0;
#endif
d.Vc(PRS,k,j,i) = d.Vc(RHO,k,j,i)/cs2*gammaIdeal;
#ifdef ISOTHERMAL
d.Vc(PRS,k,j,i) = d.Vc(RHO,k,j,i)*cs2/gammaIdeal;
#endif
d.Vc(VX1,k,j,i) = 1e-5*sin(2.0*M_PI*(y+4*z));
d.Vc(VX2,k,j,i) = shear*x;
d.Vc(VX3,k,j,i) = 0.0;
......
......@@ -111,7 +111,9 @@ void Setup::InitFlow(DataBlock &data) {
#else
d.Vc(RHO,k,j,i) = 1.0;
#endif
d.Vc(PRS,k,j,i) = d.Vc(RHO,k,j,i)/cs2*gammaIdeal;
#ifndef ISOTHERMAL
d.Vc(PRS,k,j,i) = d.Vc(RHO,k,j,i)*cs2/gammaIdeal;
#endif
d.Vc(VX1,k,j,i) = 1e-5*sin(2.0*M_PI*(y+2*z));
d.Vc(VX2,k,j,i) = shear*x;
d.Vc(VX3,k,j,i) = 0.0;
......
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