vorz_data_file = os.path.join(vorz_folder,
'vorz_' + time_instance + '.csv')
q_data_file = os.path.join(q_folder,
'q_' + time_instance + '.csv')
wgeo_data_file = os.path.join(wgeo_folder,
'wgeo_' + time_instance + '.csv')
vorz_data = read_sfield(vorz_data_file)
# q_data = read_sfield(q_data_file)
wgeo_data, w_centroid = read_wgeo(wgeo_data_file)
#-------------------------------------
grid_x, grid_y, grid_vz = grid_vorz(window, resolution,
vorz_data)
# grid_x, grid_y, grid_q = grid_vorz(window, resolution, q_data)
#-----write wing centroid history-----
centroid_ti = [
str(timei),
str(w_centroid[0]),
str(w_centroid[1])
]
fdata_si.append(grid_vz)
wdata_si.append(wgeo_data)
#----------------------------------------
fdata_all.append(fdata_si)
wdata_all.append(wdata_si)
# -------------plot all vortices-----------
field_plot(windows, fdata_all, wdata_all, marks, oimage_file, 'save')
plt.close()
# ----------------------------------------------
ufield_data_file = os.path.join(ufield_folder,
'ufield_' + time_instance + '.csv')
geop_data_file = os.path.join(wgeo_folder,
'geop_' + time_instance + '.csv')
vorz_data = read_sfield(vorz_data_file)
q_data = read_sfield(q_data_file)
ufield_data = read_vfield(ufield_data_file)
wgeo_data, w_centroid = read_wgeo(geop_data_file)
#-------------------------------------
grid_x, grid_y, grid_vz = grid_vorz(window, resolution, vorz_data)
grid_x, grid_y, grid_q = grid_vorz(window, resolution, q_data)
grid_x, grid_y, grid_ux, grid_uy = grid_ufield(window,
resolution_vector,
ufield_data)
#-----write wing centroid history-----
centroid_ti = [str(timei), str(w_centroid[0]), str(w_centroid[1])]
ufield_gridx.append(grid_x.T)
ufield_gridy.append(grid_y.T)
sImgdata.append(grid_vz.T)
sCtrdata.append(grid_q.T)
vdata.append([grid_ux.T, grid_uy.T])
wdata.append(wgeo_data)
markt.append(timei)
# -------------plot all vortices-----------
field_plot(window, ufield_gridx, ufield_gridy, sImgdata, sCtrdata, vdata,
wdata, markt, oimage_file, 'save')
plt.close()
# ----------------------------------------------
fdata_all = []
wdata_all = []
for ti in range(start_t, end_t + 1):
timei = (ti + 1) * data_time_increment
if 0.8 <= timei <= 1.1:
time_instance = str(ti).zfill(4)
#--------setting up file dirs-----------
vorz_data_file = os.path.join(vorz_folder,
'vorz_' + time_instance + '.csv')
q_data_file = os.path.join(q_folder, 'q_' + time_instance + '.csv')
wgeo_data_file = os.path.join(wgeo_folder,
'wgeo_' + time_instance + '.csv')
vorz_data = read_sfield(vorz_data_file)
q_data = read_sfield(q_data_file)
wgeo_data, w_centroid = read_wgeo(wgeo_data_file)
#-------------------------------------
grid_x, grid_y, grid_vz = grid_vorz(window, resolution, vorz_data)
grid_x, grid_y, grid_q = grid_vorz(window, resolution, q_data)
#-----write wing centroid history-----
centroid_ti = [str(timei), str(w_centroid[0]), str(w_centroid[1])]
fdata_all.append(grid_vz)
wdata_all.append(wgeo_data)
# -------------plot all vortices-----------
field_plot(window, fdata_all, wdata_all, oimage_file, 'save')
plt.close()
# ----------------------------------------------