3D finite difference modeling of controlled-source electromagnetic response in frequency domain based on a modified curl-curl equation