求助：异质结构电荷密度的计算及后处理所对应的输入文件怎么写呢？

ZhangBigPower

学习过程中看到知乎的一个不详细的教程，设置了石墨4个C原子吸附一个O原子的示例。
他提供了自洽计算输入文件和一个pp后处理的输入文件。
最后需要导出：总体系、碳原子、氧原子三者的电荷密度xsf文件，教程只能用来得到总体系的xsf文件，请问怎么改编输入文件再分别得到两个分项的电荷密度xsf文件呢？
教程提供输入文件见下：
&CONTROL
  pseudo_dir = './'
  calculation = 'scf'
  prefix = 'grap-o'/&SYSTEM
  nat= 5
  ntyp= 2
  ibrav= 0
  ecutwfc= 20.0
/
&ELECTRONS
  mixing_beta = 0.7
  conv_thr =  1.0d-8
/
&IONS
/
&CELL
/
ATOMIC_SPECIES
C   12.011  C.pbe-n-kjpaw_psl.1.0.0.UPF
O   15.999  O.pbe-n-kjpaw_psl.1.0.0.UPF
CELL_PARAMETERS angstrom
      2.46400000        0.00000000        0.00000000
      0.00000000        4.26780000        0.00000000
      0.00000000        0.00000000       12.00000000
ATOMIC_POSITIONS angstrom
C         0.00332640        0.08834346        4.92792000
C         1.23453792        2.27618845        5.26944000
C         1.23520320        0.82829462        5.01600000
C         0.00280896        3.01511534        4.98420000
O         1.23103904        1.62522092        6.50928000
K_POINTS automatic
7 5 1  0 0 0
pp后处理输入文件
&inputpp
    prefix = 'grap-o'
    outdir = '.'
    filplot = 'grap-o'
    plot_num= 0
/
&plot
    nfile = 1
    filepp(1) = 'grap-o'
    weight(1) = 1.0
    iflag = 3
    output_format = 5
    fileout = 'grap-o.xsf'
/

ZhangBigPower

已解决