I have version compute h/f with long-range coulomb and short range
pair interactions. I can sent it.
It works with Ewald method only.
The script I posted works with all type of interactions but not with
long range coulomb.
I checked it for MEAM, SW and Tersoff potentials for SiC system.
It's not a big deal to include k-space but again it will work with
Ewald method only.
German