fixed_charge_continuationでTH=0.0になる

お世話になります。

PHASE/0 2021.01にて、fixed_charge_continuationで計算を継続すると、時々出力が下記の様になります。

ENERGY EIGEN VALUE SUM. FOR     1 -TH ITER (K=  2182:  2182) =       0.00000000 <eko_old-eko_new>  NaN  NaN 0.00000D+00)

K=2182の値はこの限りでなく、どの時点で継続計算したかによります。

やろうとしている計算はBoltzTraP用の固定電荷計算で、入力ファイルはユーザーマニュアルの5.1.10節のサンプルファイルなど。

共通点を探すと、iteration_electronic=0で終えた計算から継続する時に起きるみたいです。

何か回避方法はないでしょうか。

取り急ぎ、本問題を回避するためのソースコード変更点（二次元版）をお知らせします。

src_phase/Initial_Electronic_Structure.F90の1081, 1082行を削除してください。

もしくは、添付の修正パッチを添付いたしましたので、ご利用ください。
添付ファイルをsrc_phaseディレクトリに配置してから、以下のコマンドを実行していただきますと
ソースプログラムが修正されます。

$ patch < ek_cont.diff

いずれの場合もソースプログラム修正後にmakeしてください。
修正版のekcalが生成されます。

なお、三次元版には他にも問題があるようです。
こちらは次期公開をお待ちくださいますようお願いします。

本件詳しい報告を下さったおかげで、問題を特定しその解決に至ることができました。
今後ともPHASE/0をご利用いただき、お気づきの点ございましたらコメントくださいますよう、
よろしくお願い申し上げます。