Frozen orbit is an attractive option for orbital design owing to its characteristics (its argument of pericenter and eccentricity are kept constant on an average). Solar sails are attractive solutions for massive and expensive missions. However, the solar radiation pressure effect represents an additional force on the solar sail that may greatly affect its orbital behavior in the long run. Thus, this force must be included as a perturbation force in the dynamical model for more accuracy. This study shows the calculations of initial conditions for a lunar solar sail frozen orbit. The disturbing function of the problem was developed to include the lunar gravitational field that is characterized by uneven mass distribution, third body perturbation, and the effect of solar radiation. An averaging technique was used to reduce the dynamical problem to a long period system. Lagrange planetary equations were utilized to formulate the rate of change of the argument of pericenter and eccentricity. Using the reduced system, frozen orbits for the Moon sail orbiter were constructed. The resulting frozen orbits are shown by two 3Dsurface (semimajor, eccentricity, inclination) figures. To simplify the analysis, we showed inclination–eccentricity contours for different values of semi-major axis, argument of pericenter, and values of sail lightness number.