default to z direction for unsteady

This commit is contained in:
James McClure 2019-10-22 14:27:30 +02:00
parent 64e7b1211a
commit 48b3f3b019

View File

@ -278,6 +278,12 @@ void SubPhase::Basic(){
dir_y = Fy/force_mag; dir_y = Fy/force_mag;
dir_z = Fz/force_mag; dir_z = Fz/force_mag;
} }
else {
// default to z direction
dir_x = 0.0;
dir_y = 0.0;
dir_z = 1.0;
}
if (Dm->BoundaryCondition > 0 ){ if (Dm->BoundaryCondition > 0 ){
// compute the pressure drop // compute the pressure drop
double pressure_drop = (Pressure(Nx*Ny + Nx + 1) - 1.0) / 3.0; double pressure_drop = (Pressure(Nx*Ny + Nx + 1) - 1.0) / 3.0;