A simplified nonlinear numerical model for the development of incompressible magnetohydrodynamics in the presence of a strong magnetic field B and stratification, nicknamed "Shell-Atm," is presented. In planes orthogonal to the mean field, the nonlinear incompressible dynamics is replaced by two-dimensional shell models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long integrations to obtain good statistics at moderate computational cost. The shell models of different planes are coupled by Alfvén waves propagating along B. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated by means of turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfvén waves interact nonlinearly and form turbulent spectra in the directions perpendicular and, through propagation, also parallel to the mean field. A heating function is obtained and shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona and is proportional to the aspect ratio of the loop to the -1.5 power, and characteristic properties of heating events are distributed as power laws. Cross-correlations show a delay of dissipation compared with energy content.