Question

我有几个清单：

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.dat和121.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

我无法写出来自Xs，Ys，Zs和ATOMIC_NUMBERS的信息，因为这4个列表的长度与P0和{{}不同{1}}。

我设法将P2，Xs，Ys和Zs重写为一个ATOMIC_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)):列表的长度约为50个左右。

Answer 1

使用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。

Answer 2

以下内容将创建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

同时写入几个文件，不同长度的列表元素

2 个答案: