1)，處理總電荷密度，取出某壹表面電荷密度時pp.x的輸入文件格式:
& ampinputpp
前綴= 'si '
outdir = '。/tmp '
filplot = 'sicharge '
plot_num= 0
自旋分量=0
/
& amp情節
nfile = 1
filepp(1) = 'sicharge '
重量(1) = 1.0
iflag = 2
output_format = 2
fileout = 'si.rho.dat '
e1(1) =1.0，e1(2)=1.0，e1(3) = 0.0，
e2(1) =0.0，e2(2)=0.0，e2(3) = 1.0，
nx=56，ny=40
/
prefix在前面的步驟中計算電荷密度文件時設置系統的名稱；
Outdir設置上壹次自洽計算的輸出文件所在的目錄；
Filplot設置要處理的文件的名稱，這裏我們處理的是總電荷密度，是上壹步自洽計算得到的；
Plot_num設置要處理的數據類型。這裏，因為要處理電荷密度，所以將其設置為0；
Spin_component設置是處理總電荷密度還是向上或向下旋轉。總計設置為0；1是自旋電荷密度；2是自旋向下電荷密度。
nfile設置中有幾個電荷密度文件需要處理。這裏只有壹個，所以後面只設置filepp(1)和weight(1)。如果有兩個以上，則需要指定每個的文件名和對應的權重。
Filepp(1)設置要處理的第壹個文件的文件名；權重(1)是這個文件對應的權重。
Iflag設置繪制什麽樣的地圖，因為是繪制2D等高線地圖，所以設置為2，意思是繪制2D地圖。
Output_format設置處理後數據的輸出方式，這裏設置為2，表示以plotrho.x要求的格式輸出表面的電荷密度值。
Fileout是用於輸出表面電荷密度的文件名。
E1 (1)、E1 (2)和E1 (3)是用於確定表面的第壹個向量。
E2 (1)、E2 (2)和E2 (3)是用於確定曲面的第二個向量，它們的單位是alat。
這個平面的原點由x0 (1)，x0 (2)，x0 (3)決定。
Nx和ny用於設置平面上網格的大小，nx表示沿第壹向量方向的劃分數，ny表示沿第二向量方向的劃分數。