Add time-vary & spatially varying body force for MMC simulations

[lawrenceccheung]

Is this feature request related to a problem? Please describe.

For MMC coupled wind farm simulations, we need to add a body force which varies in both z and in time. The body forces are derived from the ABLMesoscaleForcing and are horizontally homogeneous. This feature was alluded to in #1022 and #1025, and would generalize the capabilities currently available.

Describe the preferred solution

One possibility is to have the body forces be read from the abl stats netcdf file directly, which would avoid any reformatting of data (the abl stats file contains abl_meso_forcing_mom_x and abl_meso_forcing_mom_y which is what we want).

The same body forcing term would be necessary for the temperature equation here: https://github.com/Exawind/amr-wind/blob/main/amr-wind/equation_systems/temperature/source_terms/BodyForce.cpp.

Additional context

This capability is/will be necessary for FLOWMAS, AWAKEN, and other such projects.