This paper develops an edge-based finite element method for 3D CSEM modeling. The algorithm is capable of dealing with anisotropic conductivity. We adopt a total field formulation and propose a new advanced boundary condition. We use a parallelized multifrontal method to solve the system of equations. The developed method is effective in modeling complex geometry such as bathymetry.