An optimized and automated protocol for determining the location of guest sorbate molecules in highly siliceous zeolites from Si-29 INADEQUATE and H-1/Si-29 cross polarization (CP) magic-angle spinning (MAS) NMR experiments is described. With the peaks in the Si-29 MAS NMR spectrum assigned to the unique Si sites in the zeolite framework by a 2D(29)Si INADEQUATE experiment, the location of the sorbate molecule is found by systematically searching for sorbate locations for which the measured rates of H-1/Si-29 cross polarization of the different Si sites correlate linearly with H-1/Si-29 second moments calculated from H-Si distances. Due to the H-1/Si-29 cross polarization being in the "slow CP regime" for many zeolite-sorbate complexes, it is proposed that the CP rate constants are best measured by H-1/Si-29 cross polarization drain experiments, if possible, to avoid complications that may arise from fast H-1 and Si-29 T-1 rho relaxations. An algorithm for determining the sorbate molecule location is described in detail. A number of ways to effectively summarize and display the large number of solutions which typically result from a prediction of the structure from the CP MAS NMR data are presented, including estimates of the errors involved in the structure determinations. As a working example throughout this paper, the structure of the low loaded rho-dichlorobenzene/ZSM-5 complex is determined under different conditions from solid-state H-1/Si-29 CP MAS NMR data, and the solutions are shown to be in excellent agreement with the known single-crystal X-ray diffraction structure. This structure determination approach is shown to be quite insensitive to the use of relative rate constants rather than absolute values, to the detailed structure of the zeolite framework, and relatively insensitive to temperature and motions.