Development and qualification results are described for a three-dimensional, time-domain core dynamics analysis program for commercial boiling water reactors (BWRs). The program allows analysis of the reactor core with a detailed mesh division, which eliminates cal-culational ambiguity in the nuclear-thermal-hydraulic stability analysis caused by reactor core regional division. During development, emphasis was placed on high calculational speed and large memory size as attained by the latest supercomputer technology. The program consists of six major modules, namely a core neutronics module, a fuel heat conduction/transfer module, a fuel channel thermal-hydraulic module, an upper plenum/separator module, a feedwater/recirculation flow module, and a control system module. Its core neutronics module is based on the modified one-group neutron kinetics equation with the prompt jump approximation and with six delayed neutron precursor groups. The module is used to analyze one fuel bun dle of the reactor core with one mesh (region). The fuel heat conduction/transfer module solves the onedimensional heat conduction equation in the radial direction with ten nodes in the fuel pin. The fuel channel thermal-hydraulic module is based on separated three-equation, two-phase flow equations with the drift flux correlation, and it analyzes one fuel bundle of the reactor core with one channel to evaluate flow redistribution between channels precisely. Thermal margin is evaluated by using the GEXL correlation, for example, in the module. In the upper plenum/separator module, the upper plenum is modeled as a single volume in the thermal-equilibrium state and water spiraling in the separator is modeled by an effective length in the momentum equation. In the feedwater/recirculation flow module, the single-phase flow model is solved with the assumption of incompressive flow. Finally, the control system module includes the recirculation flow control minimodule, the pressure control minimodule, and the feedwater control minimodule, as well as the interlock functions, which work during a transient to allow analysis of general transient phenomena. The program was verified to provide satisfactory results within reasonable computational time based on application analysis of stability and scram phenomena in a BWR-5 type plant.