An automated system for designing a loading pattern (LP) for boiling water reactors (BWRs) given a reference LP and control rod (CR) sequence has been developed. This system employs the advanced nodal code SIMULATE-3 and a BWR LP optimization code FINELOAD-3, which uses a simple linear perturbation method and a modified Tabu search method to select potential optimized LP candidates. Both of these unique methods of FINELOAD-3 were developed to achieve an effective BWR LP optimization strategy and to have high computational efficiency. FINELOAD-3 also adjusts deep CR positions to compensate for the core reactivity deviation caused by fuel shuffling. The objective function is to maximize the end-of-cycle core reactivity while satisfying the specified thermal margins and cold shutdown margin constraints. This optimization system realized the practical application for real BWR LP design. Computer time needed to obtain an optimized LP for a typical BWR/5 octant core with 15 depletion steps is ~4 h using an engineering workstation. This system was extensively tested for real BWR reload core designs and showed that the developed LPs using this system are equivalent or better than the manually optimized LPs.