A four-dimensional (4D) intermolecular potential energy surface (IPES) for the water trimer was calculated in the subspace of the three torsional coordinates and the symmetric intermolecular stretch coordinate, employing high-level ab initio theory. Torsionally adiabatic potential energy curves for the intermolecular symmetric stretching vibration were constructed based on this IPES. They were used to calculate the symmetric stretch fundamentals and stretching-averaged interoxygen R(O ... O) distances [R-H] and [R-D], for the n = 0-6 torsional levels of (H2O)(3) and (D2O)(3). [R-H] and [R-D] increase with n up to n = 5, and decrease for n = 6. Torsionally averaged rotational constants A, B, and C of all 20 isotopomers of water trimer, for the torsional levels n =0, 5 and 6, were obtained by averaging the inverse inertia tensor over the 3D torsional wave functions. Two approaches were examined: (i) setting the interoxygen R(O ... O) distance to a fixed value, independent of torsional excitation; (ii) effectively incorporating the vibrational averaging due to the intermolecular symmetric stretching mode by using the appropriate [R-H] and [R-D] values to define the R(O ... O) distances for the isotopomers in the torsional state n. Both approaches yielded n=0, 5 rotational constants in good agreement with experiment. However, only approach (ii) reproduced the experimentally observed decrease in the rotational constants A and B upon 5 <--0 torsional excitation. Fixing the R(O ... O) distances to R-H and R-D values obtained by fitting the 3D torsionally averaged rotational constants to the experimental values for (H2O)(3) and (D2O)(3), only marginally improved the agreement with experiment for other isotopomers.