Consolidate aerosol number concentration arrays

[climbfuji]

Originally raised as https://github.com/NCAR/FV3/issues/99:

With the updates to the CPT suite, three arrays holding aerosol information are added to the Tbd DDT by EMC:

Tbd%in_nm(:,:) ! IN number concentration
Tbd%ccn_nm(:,:) ! CCN number concentration
Tbd%aer_nm(:,:) ! GOCART aerosol climo
These are used by Morrison-Gettelman microphysics at the moment, Tbd%aer_nm is mutually exclusive with Tbd%in_nm + Tbd%ccn_nm. On the other hand, we have Statein%qgrs(:,:,ntwa) and Statein%qgrs(:,:,ntia) (and corresponding Stateout%gq0(:,:,ntwa) and Stateout%gq0(:,:,nita)) that are used by the aerosol-aware Thompson microphysics scheme, reading in its own version of the aerosol climatology. Thompson also uses the Coupling%nwfa2d(:) and Coupling%nifa2d(:) arrays to hold surface emission rates for the two species.

Questions to be answered:

Are the values in EMC's aerosol climatology from the same source as those used by Thompson MP?
If the values derived from the climatology and used by the MP schemes compatible?
Can we use the same arrays for both Thompson and MG at least? How does MG use the Statein and Stateout arrays?