A charge-scaling implementation of the variational electrostatic projection method