Dear Petr Khomyakov,
Thank you very much for your prompt reply. Attached are the .py files with doping, n-doping and p-doping.
Strangely, the EDP plots obtained are almost the same, with no influence of the doping type and doping level.
BTW, I am using ATK2019.03.
Best regards,
lknife_ATK