An algorithm for simulating multisite streamflow or water quality data for hydrologic management and planning studies is demonstrated. Key features are: (1) Straightforward consideration of cross-linkage and Markovian correlations between gaging sites, and (2) use of time-independent residual data having low average skewness for parameter estimation. The FORTRAN IV code is available upon request.