Top row: time profiles for the internal, magnetic, kinetic, and gravitational energy density variations relative to the initial state, total (internal + magnetic + kinetic + gravitational) energy density difference, and energy density provided by the driver. All the quantities are volume averaged for the whole computational domain. From left to right: ColdI, ColdR, and ColdV models. The contributions from the Poynting flux and energy flux due to the plasma displacement through the side boundaries, have been incorporated to the magnetic and internal energy densities. Bottom row: time averaged 1D power spectra of kinetic energy, magnetic energy, and pressure at the apex, averaged over the last oscillation period, for the ColdI, ColdR, ColdV, and ColdV2 models.