startcount = 1 r = pd.DataFrame allmatrix = [] for j in range(len(pathstring)): for p in sorted( Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')): print(p) try: r_single = mu.getAvgR( str(p)) #, 'zero', 'cylinder_shell.out', p.parts[-2][2]) inseg = mu.gettubeILTT(p, 'adsorbed') # choose itube/otube/iotube outseg = mu.gettubeOLTT( p, 'adsorbed' ) # [[# adsorbed chain, fluc], [# adsorbed seg, fluc]] except ValueError: continue # print(inseg) # print('#adsorbed chain', inseg[0][0]) # print('#adsorbed chain fluc', inseg[0][1]) # print('#adsorbed seg', inseg[1][0]) # print('#adsorbed seg fluc', inseg[1][1]) # print(outseg) labelstring = 'ph = {:d}'.format(int(p.parts[-2][2])) print(p.parts) print(p.parts[-2]) print(p.parts[-2][2])
seg_io_loop_errup = [] num_io_loop = [] num_io_loop_errdown = [] num_io_loop_errup = [] ph = [] for p in sorted( Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')): # zip(pathstring, linecolor) print(p) #io_loops = mu.gettubeIOLTT(p, 'loops') #choose itube/otube/iotube #io_tails = mu.gettubeIOLTT(p, 'tails') #io_trains = mu.gettubeIOLTT(p, 'trains') i_loops = mu.gettubeILTT(p, 'loops') #choose itube/otube/iotube o_loops = mu.gettubeOLTT(p, 'loops') #choose itube/otube/iotube io_loops = i_loops + o_loops ## Input: fin - file in ## restype - type of result. Can be 'loops', 'tails', 'trains' or 'adsorbed' ## Output: 2D array of results on the form ## [number of x] [variation (fluctuation)], ## [length of x] [variation], ## [number of segments in x] [variation] ## in the case of x = 'loops', 'tails' or 'trains', or on the form ## [number of adsorbed chains] [variation], ## [number of adsorbed segments] [variation] ## in the case of 'adsorbed' ## Note: loops, tails and trains are average number per (adsorbed?) chain! if j == 0:
for p in sorted( Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')): # zip(pathstring, linecolor) print(p) io_loops = mu.gettubeIOLTT(p, 'loops') #choose itube/otube/iotube io_tails = mu.gettubeIOLTT(p, 'tails') io_trains = mu.gettubeIOLTT(p, 'trains') i_loops = mu.gettubeILTT(p, 'loops') #choose itube/otube/iotube i_tails = mu.gettubeILTT(p, 'tails') i_trains = mu.gettubeILTT(p, 'trains') o_loops = mu.gettubeOLTT(p, 'loops') #choose itube/otube/iotube o_tails = mu.gettubeOLTT(p, 'tails') o_trains = mu.gettubeOLTT(p, 'trains') ## Input: fin - file in ## restype - type of result. Can be 'loops', 'tails', 'trains' or 'adsorbed' ## Output: 2D array of results on the form ## [number of x] [variation (fluctuation)], ## [length of x] [variation], ## [number of segments in x] [variation] ## in the case of x = 'loops', 'tails' or 'trains', or on the form ## [number of adsorbed chains] [variation], ## [number of adsorbed segments] [variation] ## in the case of 'adsorbed' ## Note: loops, tails and trains are average number per (adsorbed?) chain!
seg_io_loop_errup = [] num_io_loop = [] num_io_loop_errdown = [] num_io_loop_errup = [] ph = [] for p in sorted(Path(str(pathstring[j])).glob('zero/ph*/cylinder_shell.out')): # zip(pathstring, linecolor) print(p) #io_loops = mu.gettubeIOLTT(p, 'loops') #choose itube/otube/iotube #io_tails = mu.gettubeIOLTT(p, 'tails') #io_trains = mu.gettubeIOLTT(p, 'trains') i_loops = mu.gettubeILTT(p, 'adsorbed') #choose itube/otube/iotube o_loops = mu.gettubeOLTT(p, 'adsorbed') #choose itube/otube/iotube io_loops = i_loops + o_loops if j == 0: titlestring = 'neg, out' elif j == 1: titlestring = 'pos, out' elif j == 2: titlestring = 'neg, in' elif j == 3: titlestring = 'pos, in' #titlestring = '{:s}'.format(str(p.parts[-4])) labelstring = 'pH = {:d}'.format(int(p.parts[-2][2]))