A partial-wave method is developed to deal with small molecules dominated by a central atom as an extension of earlier single-center methods. In particular, a model potential for the water molecule is expanded over a basis of spherical harmonics. A finite element method is employed to generate local polynomial functions in subintervals over a finite range for the radial variable.The angular parts of the system are represented by spherical harmonics. The problem of Stark resonances is treated with the exterior complex scaling method which incorporates a wavefunction discontinuity at the scaling radius. The resultant non-hermitian matrix eigenvalue problem yields resonance positions and widths (decay rates). We present these DC Stark shifts and exponential decay rates for the valence orbitals 1b 1 , 3a 1 , and the bonding orbital 1b 2 . Furthermore, comparison is made with total molecular decay rates and DC shifts obtained recently within the Hartree-Fock and coupled-cluster approaches.