Dear Stefano de Gironcoli, Thank you very much. It really does work! By reading your detailed README file,I understand this method further.But I also found a little bug in the espresso-2.1. When plus u is used in the calculation, the projwfc.x can not run successfully.An error "array has already allocated" appears.If plus u is not used,projwfc.x can run without any interruption. I just insert a line "deallocate (swfcatom)" in the projwfc.f90 to slove it.But I think you may have a better way to fix this problem. Thanks a lot.
Best regards, Sincerely, Ruijuan Xiao Institute of Physics, Chinese Academy of Sciences
