(Both options are implemented in the same .F file and you can try to fix the bug, by comparing the file to siesta 2)
1. SlabDipoleCorrection .true.
2. %block ExternalElectricField
0.0 0.0 1.0 V/Ang
%endblock ExternalElectricField
