Allow N tracers to be advected

[FlorianDeconinck]

GEOS-FP advects 28 tracers. DSL port does 9 hard coded tracers.

Parent: https://github.com/GEOS-ESM/SMT-Nebulae/issues/55

---

• [ ] ndsl dynamical core can advect N tracers

[FlorianDeconinck]

GEOS breakdown:

• Bundle of tracers, some are flagged friendly to dynamics, the flagged are all advected
  • Within the dynamics tracers, FV_State mode singles out the water/cloud species
  • AdvCore is a standalone advection, which borrows some (how much?) code from FV and is used to
  • Do offline advection for Chemistry
  • Do advection of the non-water species on a separate DT
• nwat is used. ncnst is not

[FlorianDeconinck]

pyFV3 breakdown:

• Code is split between TracerAdvection and Remapping:mapn_tracers.
• The "advectable" tracers are denoted by NQ in constants.py, which means order of tracers in the array matters (sic)
• TracerAdvection actually names the tracers in its code, NQ is passed but unused
• Mapn_tracers uses NQ to loop over the tracers and remap them

[FlorianDeconinck]

Potential solutions:

• Replicate the "bundle" of tracers via DataClass/Register combo, with a way to show water species and non-water species on top of names
• Allow the initialization of stencils by name of species so that we can do a simple for q_name in advectable_tracers: stencil[q_name](). That way on orchestration we will be able to use the nounroll and have a clean and faster SDFG