Dear Paolo Can you please explain what the variable "zdotc" in elphon.f90 exactly does? I think it calculates the derivative of the potential. It is declared as an external variable but i couldnot figure out the location it has been calculated. Can you show me the the path as well where it is calculated ? Thanks Ajit