timei = (ti + 1) * data_time_increment
if 0.0 <= timei <= 1.05:
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')
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)
fdata_si = []
wdata_si = []
for ti in range(start_t, end_t + 1):
timei = (ti + 1) * data_time_increment
if 0.7 <= timei <= 0.8:
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_si.append(grid_vz)
wdata_si.append(wgeo_data)
#----------------------------------------
#------read data for each pa-------------
for pai in pa:
if pai == 45.0:
pf = 0.125
else:
pf = 0.25