3D electromagnetic modeling using staggered finite differences