fix sign error

This commit is contained in:
James McClure
2019-10-25 10:53:47 +02:00
parent 48b3f3b019
commit ac5ec4492e

View File

@@ -288,7 +288,7 @@ void SubPhase::Basic(){
// compute the pressure drop
double pressure_drop = (Pressure(Nx*Ny + Nx + 1) - 1.0) / 3.0;
double length = ((Nz-2)*Dm->nprocz());
force_mag += pressure_drop/length;
force_mag -= pressure_drop/length;
}
if (force_mag == 0.0){
// default to z direction