Implemented AdvectionSrcForMom_EB.

[skang67]

This PR implements RHS calculation for momentum variables (AdvectionSrcForMom_EB).

1. The highest upwind order is determined using flag.isCovered(). If the specified advection scheme is not applicable, it is reduced to CENTERED2 or UPWIND3.
2. Fluxes for staggered grids are computed using EB datasets (eb_aux_).

Validation:

• When the area and volume fractions are commented out, the function AdvectionSrcForMom_EB is backward compatible with AdvectionSrcForMom_ConstantDz, except for WENO schemes.
• AdvectionSrcForMom_EB has not been validated for the EB test. We need SRD for momentum variables to have a stable simulation.