import sys
import numpy
import dakota_utils
from build_mesh import build_mesh
from apame_utils import run_case
x,varnames,=dakota_utils.read_input(sys.argv[1])
filename = 'wing.vtp'
build_mesh(filename,
semi_span=x[0],
root_chord=x[1],
tip_chord=x[2],
root_m=x[3],
root_p=x[4],
tip_m=x[5],
tip_p=x[6],
t=0.1,
nu=50,
nv=21)
Patm = 101325.#Pa
rho = 1.225#kg.m-3
c = numpy.sqrt(1.4*Patm/rho)
Mach = 0.17
airspeed = c*Mach
run_case(filename,
wake_length=20.,
alpha=[x[7]],
tip_chord=1./1.1
root_m=0.01
root_p=0.4
tip_m=0.01
tip_p=0.4
t = 0.11
filename = 'wing.vtp'
build_mesh(filename,
semi_span,
root_chord,
tip_chord,
root_m,
root_p,
tip_m,
tip_p,
t=t,
nu=50,
nv=41)
Patm = 101325.#Pa
rho = 1.225#kg.m-3
c = numpy.sqrt(1.4*Patm/rho)
Mach = 0.15
airspeed = c*Mach
mu=1.5e-5
Re = rho*airspeed*root_chord/mu
print "Re=",Re
print "airspeed=",airspeed*3.6