else: # nothing given print("read_parts error: Invalid commandline arguments.") print("Usage: ") print(" ./read_parts.py <./path/to/sim/output> <start_time>") print(" or") print(" ./read_parts.py <./path/to/sim/output>") sys.exit() # initialize the reader times = bbparts.init(data_dir) # visit all outputted time values for time in times: # open the CGNS file for this particular output time bbparts.open(time) # read the CGNS file t = bbparts.read_time() n = bbparts.read_nparts() (x, y, z) = bbparts.read_part_position() (u, v, w) = bbparts.read_part_velocity() print("time = ", time, "t =", t, "n =", n) print(u) sys.exit() # close the CGNS file bbparts.close()
# number of time outputs nt = 0; for realization in ensemble: timeseries = bb.init(realization + "/output")[int(timestart/DT_out):] # open final output in timeseries and read time for comparison bb.open(timeseries[-1]) t_tmp = bb.read_time() if t_tmp < t_end: t_end = t_tmp nt = len(timeseries) # close this output bb.close() # overwrite number of time outputs to read (for testing) #nt = 300 #t_end = 3 #################################################################### # store particle position data at initial time t_init #################################################################### timeseries = bb.init(ensemble[0] + "/output")[int(timestart/DT_out):] bb.open(timeseries[0]) # using time step 0 for now (X_init, Y_init, Z_init) = bb.read_part_position() T_init = bb.read_time() np = len(X_init) # also store particle number bb.close()