All singularities on the solvent excluded surface are computed analytically. The data structures of all components of the molecular surfaces are pre-computed. The advancing-front method is used to mesh each spherical SES patch. Each vertex of the mesh lies exactly on the molecular surface. An efficient algorithm for filling molecular inner holes is proposed.