From 48b3f3b0190b2f17cf1cf2b626769a5df7328914 Mon Sep 17 00:00:00 2001
From: James McClure <mcclurej@vt.edu>
Date: Tue, 22 Oct 2019 14:27:30 +0200
Subject: [PATCH] default to z direction for unsteady

---
 analysis/SubPhase.cpp | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/analysis/SubPhase.cpp b/analysis/SubPhase.cpp
index 18e67aa9..d17187de 100644
--- a/analysis/SubPhase.cpp
+++ b/analysis/SubPhase.cpp
@@ -278,6 +278,12 @@ void SubPhase::Basic(){
 			dir_y = Fy/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 ){
 			// compute the pressure drop
 			double pressure_drop = (Pressure(Nx*Ny + Nx + 1) - 1.0) / 3.0;