讀取輸入并設置要計算的量的標志
1.1) 讀取 pw.x 寫入的所有量
1.2) 讀取贗勢數據

確定必須計算的內容。
2.1) 如果還沒有在disk上，計算 q 個點的網格和
所有 q 點的所有modes并保存在disk上 (SD means save on disk)
2.2) 如果要求image 并行化，則在images之間劃分工作

進行恢復運行檢查，檢查 .xml 文件中已有的內容，并將以完成的部分設置為TRUE

開始對于q點的主循環：do_phonon.f90
4.1) 計算所有不依賴于系統響應的量 initialize_ph()
4.2）檢查是否需要進行能帶計算并進行。
注意：以下幾點僅在 q 為 Gamma 時執行。
4.3) 在頻率上開始循環
4.3.1) 計算極化率作為 iu SD 的函數
4.4) 計算對電場的響應
4.5) 計算 epsilon 和 SD
4.6）計算zeu和SD
4.7) 計算電光系數和SD
4.8) 計算 E 的二階響應
4.9) 計算拉曼張量和 SD
結束注

在不可約表示上開始循環
5.1) 計算對不可約表示的響應
代碼解析