return xyz,direction,speed (xyz,direction,speed) = load_input_file("input.txt") mesh = Mesh(mesh="../geom/conv_pipe_mesh.h5m") # load dag geometry load("../geom/conv_pipe_geom.h5m") # discretise geom onto mesh cell_fracs = discretize_geom(mesh) mesh.vol_id = IMeshTag(1,int) mesh.source_strength = IMeshTag(1,float) id_list = [] for element in cell_fracs: id_list.append(element["cell"]) mesh.vol_id[:]=id_list # loop along path, use point in volume to determine the vol id volumes = get_all_cells(xyz,direction,speed) # loop along path, use point in volume to determine the vol id # but tag each tet in the range with right source strength dump = get_all_cells_detailed(xyz,direction,speed,mesh) sys.exit()