This work employs a Landau-Ginzburg-type nonlocal rheology model to account for shear localization in a wall-bounded dense granular flow. The configuration is a 3D shear cell in which the bottom bumpy wall moves at a constant speed, while a load pressure is applied at the top bumpy wall, with flat but frictional lateral walls. At a fixed pressure, shear zones transit from the top to the bottom when increasing lateral wall friction coefficient. With a quasi-2D model simplification, asymptotic solutions for fluidization order parameters near the top and bottom boundaries are sought separately. Both solutions are the Airy function in terms of a depth coordinate scaled by a characteristic length which measures the width of the corresponding shear zone. The theoretical predictions for the shear zone widths against lateral wall friction coefficient and load pressure agree well with data extracted from particle-based simulation for the flow.