同时写入多个文件,不同长度的列表元素

Writing simultaneously, into several files, the elements of lists of different length

我有几个列表:

VOLUMES =  ['119.823364', '121.143469']
P0 =  ['4.97568007', '4.98494429']
P2 =  ['16.76591397', '16.88768068']
Xs =  ['0.000000000000E+00', '3.333333333333E-01', '-4.090760942850E-01', '0.000000000000E+00', '3.333333333333E-01', '-4.093755657782E-01']
Ys =  ['0.000000000000E+00', '-3.333333333333E-01', '-3.333333333333E-01', '0.000000000000E+00', '-3.333333333333E-01', '-3.333333333333E-01']
Zs =  ['0.000000000000E+00', '-8.333333333333E-02', '-8.333333333333E-02', '0.000000000000E+00', '-8.333333333333E-02', '-8.333333333333E-02']
ATOMIC_NUMBERS =  ['20', '6', '8', '20', '6', '8']

我想生成2个文件,以VOLUMES列表中的项目命名:119.823364.dat121.143469.dat,例如每个文件包含以下内容:

119.823364.dat 文件:

some stuff
other stuff
4.97568007   16.76591397
3
20 0.000000000000E+00    0.000000000000E+00   0.000000000000E+00
6  3.333333333333E-01   -3.333333333333E-01  -8.333333333333E-02
8 -4.090760942850E-01   -3.333333333333E-01  -8.333333333333E-02
other stuff
some other stuff

121.143469.dat 文件:

some stuff
other stuff
4.98494429  16.88768068
3
20 0.000000000000E+00    0.000000000000E+00   0.000000000000E+00
6  3.333333333333E-01   -3.333333333333E-01  -8.333333333333E-02
8 -4.093755657782E-01   -3.333333333333E-01  -8.333333333333E-02
other stuff
some other stuff

存在以下问题:

len(VOLUMES) = len(P0) = len(P2) = 2

但是:

len(Xs) = len(Ys) = len(Zs) = 6

我已经完成了第一部分:

# Remove *.dat files, to clean first: 
for f in glob.glob("*.dat"):
    os.remove(f)

# Create the files:
filenames = []
for V in VOLUMES:
    filename = "{}.dat".format(V)
    print 'filename = ', filename
    filenames.append(filename)
print filenames

# Write to files:
for i in xrange(len(P0)):
       with open(filenames[i],'w') as f:
        f.write("""some stuff
other stuff\n""")
        f.write("{} {}\n".format(P0[i], P2[i]))
        f.write("{}\n".format(N_atom_irreducible_unit))

创建以下内容:

119.823364.dat 文件:

some stuff
other stuff
4.97568007 16.76591397
3

121.143469.dat 文件:

some stuff
other stuff
4.98494429 16.88768068
3

我无法写出来自 XsYsZsATOMIC_NUMBERS 的信息,因为这 4 个列表的长度不同于 P0P2.

我设法将 XsYsZsATOMIC_NUMBERS 重写成一个 list of list of lists:

for index_vol in range(len(VOLUMES)):
  for index in range(len(ATOMIC_NUMBERS)):
    atoms_per_frame = [ATOMIC_NUMBERS[index], Xs[index], Ys[index], Zs[index]]
    atoms_all_frames[index_vol].append(atoms_per_frame)

print atoms_all_frames  

打印以下内容:

[[['20', '0.000000000000E+00', '0.000000000000E+00', '0.000000000000E+00'], ['6', '3.333333333333E-01', '-3.333333333333E-01', '-8.333333333333E-02'], ['8', '-4.090760942850E-01', '-3.333333333333E-01', '-8.333333333333E-02'], ['20', '0.000000000000E+00', '0.000000000000E+00', '0.000000000000E+00'], ['6', '3.333333333333E-01', '-3.333333333333E-01', '-8.333333333333E-02'], ['8', '-4.093755657782E-01', '-3.333333333333E-01', '-8.333333333333E-02']], [['20', '0.000000000000E+00', '0.000000000000E+00', '0.000000000000E+00'], ['6', '3.333333333333E-01', '-3.333333333333E-01', '-8.333333333333E-02'], ['8', '-4.090760942850E-01', '-3.333333333333E-01', '-8.333333333333E-02'], ['20', '0.000000000000E+00', '0.000000000000E+00', '0.000000000000E+00'], ['6', '3.333333333333E-01', '-3.333333333333E-01', '-8.333333333333E-02'], ['8', '-4.093755657782E-01', '-3.333333333333E-01', '-8.333333333333E-02']]]

我不知道创建这个 list of list of lists 是否是能够与 for i in xrange(len(P0)): 一起循环的解决方案,但我无法做到这一点。

实际上,VOLUMES 列表的长度约为 50 项。

以下将创建 119.823364.dat:

VOLUMES = ['119.823364', '121.143469']
P0 = ['4.97568007', '4.98494429']
P2 = ['16.76591397', '16.88768068']
Xs = ['0.000000000000E+00', '3.333333333333E-01', '-4.090760942850E-01', 
      '0.000000000000E+00', '3.333333333333E-01', '-4.093755657782E-01']
Ys = ['0.000000000000E+00', '-3.333333333333E-01', '-3.333333333333E-01', 
      '0.000000000000E+00', '-3.333333333333E-01', '-3.333333333333E-01']
Zs = ['0.000000000000E+00', '-8.333333333333E-02', '-8.333333333333E-02', 
      '0.000000000000E+00', '-8.333333333333E-02', '-8.333333333333E-02']
ATOMIC_NUMBERS = ['20', '6', '8', '20', '6', '8']

with open('{}.dat'.format(VOLUMES[0]), 'w+') as file:
    file.write('some stuff \n')
    file.write('other stuff \n')
    for item in P0:
        file.write('{}   '.format(item))
    file.write('\n')
    file.write('3 \n')
    for i in range(3):
        file.write('{} {} {} {} \n'.format(ATOMIC_NUMBERS[i], 
                                           Xs[i], Ys[i], Zs[i]))
    file.write('other stuff')

您将遵循类似的过程来创建 121.143469.dat

使用 zip 构建您的线条,并将结果拆分为 len(rows) / len(volumes) 大小的块。然后,将每个块写入其各自的文件。

headers = list(zip(P0, P2))
rows = [row for row in zip(ATOMIC_NUMBERS, Xs, Ys, Zs)]
interval = int(len(rows) / len(VOLUMES))

for block_i, vol_i in zip(range(0, len(rows), interval), range(len(VOLUMES))):
    # Create the lines for the file
    lines = [' '.join(headers[vol_i]), '3']
    lines += [' '.join(row) for row in rows[block_i : block_i + interval]]
    # Write the file
    with open(VOLUMES[vol_i] + '.dat', 'w') as f:
        # Preceding lines
        f.write('some stuff\nother stuff')
        # Lines of data
        for line in lines:
            f.write(line + '\n')
        # Trailing lines
        f.write('other stuff\nsome other stuff')

文件 119.823364.dat 将包含:

some stuff
some other stuff
4.97568007 16.76591397
3
20 0.000000000000E+00 0.000000000000E+00 0.000000000000E+00
6 3.333333333333E-01 -3.333333333333E-01 -8.333333333333E-02
8 -4.090760942850E-01 -3.333333333333E-01 -8.333333333333E-02
other stuff
some other stuff

请注意,此方法是动态的,适用于任何长度的 VOLUMES