From 150f9c3f5eb11d9183ac0952427fd2507f7bf8f8 Mon Sep 17 00:00:00 2001
From: mjburton <mjburton@mit.edu>
Date: Thu, 26 Oct 2017 09:44:59 -0400
Subject: [PATCH] add tail testing script

---
 gpkitmodels/GP/aircraft/tail/tail_tests.py | 46 ++++++++++++++++++++++
 1 file changed, 46 insertions(+)
 create mode 100644 gpkitmodels/GP/aircraft/tail/tail_tests.py

diff --git a/gpkitmodels/GP/aircraft/tail/tail_tests.py b/gpkitmodels/GP/aircraft/tail/tail_tests.py
new file mode 100644
index 00000000..c2cc5ba7
--- /dev/null
+++ b/gpkitmodels/GP/aircraft/tail/tail_tests.py
@@ -0,0 +1,46 @@
+" test tail models "
+from gpkitmodels.GP.aircraft.tail.horizontal_tail import HorizontalTail
+from gpkitmodels.GP.aircraft.tail.vertical_tail import VerticalTail
+from gpkitmodels.GP.aircraft.tail.empennage import Empennage
+from gpkitmodels.GP.aircraft.wing.wing_test import FlightState
+from gpkit import Model
+
+def test_htail():
+
+    ht = HorizontalTail()
+    fs = FlightState()
+    ht.substitutions.update({ht.topvar("W"): 5, "m_h": 0.01, "AR": 4})
+    perf = ht.flight_model(ht, fs)
+
+    m = Model(perf["C_d"], [ht, perf])
+    m.solve("mosek")
+
+def test_vtail():
+
+    ht = VerticalTail()
+    fs = FlightState()
+    ht.substitutions.update({ht.topvar("W"): 5, "AR": 3})
+    perf = ht.flight_model(ht, fs)
+
+    m = Model(perf["C_d"], [ht, perf])
+    m.solve("mosek")
+
+def test_emp():
+
+    emp = Empennage()
+    fs = FlightState()
+    emp.substitutions.update({emp.topvar("W"): 10, "l": 5,
+                              emp.htail.planform["AR"]: 4,
+                              emp.vtail.planform["AR"]: 4})
+    htperf = emp.htail.flight_model(emp.htail, fs)
+    vtperf = emp.vtail.flight_model(emp.vtail, fs)
+    tbperf = emp.tailboom.flight_model(fs)
+
+    m = Model(htperf["C_d"] + vtperf["C_d"] + tbperf["C_f"],
+              [emp, fs, htperf, vtperf, tbperf])
+    m.solve("mosek")
+
+if __name__ == "__main__":
+    test_htail()
+    test_vtail()
+    test_emp()