The new streamflow-routing package (SFR1) (Prudic and others, 2004) has five options for computing the stage from stream flow in MODFLOW-2000 simulations. The USGS collected more than 10,000 field measurements of width, depth, and velocity on the Gila River and selected tributaries in east-central Arizona from the early 1900s to the present. Inspection of graphs of these field measurements at numerous gages indicates that the depth- and width-flow relationships change across the range of flow rates. Field inspection of the stream channel morphology confirms that distinct channel configurations will dominate discrete ranges of flow. The three newly available options (Manning Equation–8-PointChannel Configuration, Power Function, and Table of Values) provide for distinct depth- and width- flow relationships within distinct ranges of flow. The Power Function option was found to be less than adequately descriptive of the relationships observed at the gages considered here. The data compiled from thousands of USGS field measurements provide a substantial foundation for application of the remaining two newly available options. In this case, the Table of Values option was most convenient to apply and appears to be adequate in its characterization of the depth- and width- flow relationships. Suitable SFR1 input is identified with the Table of Values option using the extensive USGS data sets available for the gaging stations selected for this evaluation.