]
#---------------------------------------
cwd = os.getcwd()
data_dir = os.path.join(os.path.dirname(os.path.dirname(os.path.dirname(cwd))),
'wake_convection_data/3_SIM_RESULTS')
image_out_path = cwd
#---------------------------------------
mcf_array = []
for cfi in cfd_data_list:
cfd_datai = os.path.join(data_dir, cfi)
mcf_arrayi = read_cfd_data(cfd_datai)
mcf_array.append(mcf_arrayi)
#---------------------------------------
with open('meancf_all.dat', 'w') as f:
for data_name, mvalues in zip(cfd_data_list, mcf_array):
f.write("%s:\n" % data_name)
f.write(
"mcl_s = %s, mcl_w = %s, ratio_l = %s, mcd_s = %s, mcd_w = %s, ratio_d = %s, ratio_ld = %s\n"
% (
'{0:.8g}'.format(mvalues[0]),
'{0:.8g}'.format(mvalues[1]),
'{0:.8g}'.format(mvalues[2]),
'{0:.8g}'.format(mvalues[3]),
'{0:.8g}'.format(mvalues[4]),
'{0:.8g}'.format(mvalues[5]),
'{0:.8g}'.format(mvalues[6]),
))
cf_plotter(x_data, mcf_array, markc, markEffects, legends, x_range, y_range,
y_label, image_out_path)
# x_range = 'all'
# y_range = 'all'
x_range = [1, 8]
cl_range = [0.9, 1.7]
cd_range = [1.0, 2.8]
pf_range = [0.4, 1.0]
y_range = [cl_range, cd_range, pf_range]
y_label = [r'$\bar{C_L}$', r'$\bar{C_D}$', r'$\frac{1}{P^\ast}$']
#---------------------------------------
cwd = os.getcwd()
kinematics_dir = os.path.join(
os.path.dirname(os.path.dirname(os.path.dirname(cwd))),
'geometry_effect_fw_data/1_kinematic_cases/3dbm_kinematic_cases')
data_dir = os.path.join(os.path.dirname(os.path.dirname(os.path.dirname(cwd))),
'geometry_effect_fw_data/5_SIM_RESULTS')
image_out_path = cwd
#---------------------------------------
CF_file_names = [f.name for f in os.scandir(data_dir) if f.is_file()]
mcf_array = []
for cfi in cfd_data_list:
for CF_name in CF_file_names:
if CF_name.startswith(cfi):
kinematics_datai = os.path.join(kinematics_dir, CF_name)
cfd_datai = os.path.join(data_dir, CF_name)
u2, karr, dkarr, ddkarr = read_kinematics_data(kinematics_datai)
mcf_arrayi = read_cfd_data(cfd_datai, u2, karr, dkarr)
mcf_array.append(mcf_arrayi)
#---------------------------------------
cf_plotter(x_data, mcf_array, legends, x_range, y_range, y_label,
image_out_path, title, x_label)
cl_range = [0.9, 1.7]
cd_range = [1.0, 2.8]
pf_range = [0.4, 1.0]
y_range = [cl_range, cd_range, pf_range]
y_label = [r'$\bar{C_L}$', r'$\bar{C_D}$', r'$\frac{1}{P^\ast}$']
#---------------------------------------
cwd = os.getcwd()
kinematics_dir = os.path.join(
os.path.dirname(os.path.dirname(os.path.dirname(cwd))),
'geometry_effect_fw_data/1_kinematic_cases/3dbm_kinematic_cases')
data_dir = os.path.join(os.path.dirname(os.path.dirname(os.path.dirname(cwd))),
'geometry_effect_fw_data/5_SIM_RESULTS')
image_out_path = cwd
#---------------------------------------
CF_file_names = [f.name for f in os.scandir(data_dir) if f.is_file()]
mcf_array = []
Ro = []
for cfi in cfd_data_list:
for CF_name in CF_file_names:
if CF_name.startswith(cfi):
Roi = float(CF_name.split('Ro')[-1])
Ro.append(Roi)
kinematics_datai = os.path.join(kinematics_dir, CF_name)
cfd_datai = os.path.join(data_dir, CF_name)
u2, karr, dkarr, ddkarr = read_kinematics_data(kinematics_datai)
mcf_arrayi = read_cfd_data(cfd_datai, u2, karr, dkarr)
mcf_array.append(mcf_arrayi)
#---------------------------------------
cf_plotter(Ro, no_x, mcf_array, legends, x_range, y_range, y_label,
image_out_path)