Given a surface with isolated ordinary double points in P^3 over F_q and an equisingular liftof its equation to Z_q, where q = p^a for p > 3, we give an algorithm that computes its zeta function via p-adic cohomology, using results of Dimca and Saito on singular hypersurfaces in characteristic zero.