This article presents a novel ramp rate control and active power smoothing and shifting methodology for net-load profiles in large power distribution networks where high amount of photovoltaic (PV) penetration levels exist. The novelty is that the proposed methodology uses dynamic state-of-charge (SoC) management, energy storage optimal use of any given size and fast-Fourier transform-based reference curves that adequately fit energy storage system (ESS) active and reactive power capabilities. Then, a least square minimization-based optimization method is applied taking both ancillary applications into consideration and to maintain SoC within limits based on ESS capacity. The methodology is applied on a 8500 nodes US power grid model using CYMEDist simulation platform and applied on a real US power distribution system. Results show significant level of net-load ramp rate reduction as well as smoothing of fast fluctuations (up to 60%) due to both PV intermittency and load changes. The approach is also capable of breaking large systems into virtual local net-load locations where ESS's are installed to mitigate the severity of the anticipated shape of system net-loads.