Implicit buoyancy?

[pdziekan]

Buoyancy is supposed to be added implicitly. If so, shouldn't it be included in normalize_vip?

TODO:

• add implicit velocity forcings through vip_rhs, not rhs
• add correct vip_rhs_impl_fnlz
• add correct normalize_vip
• make velocity absorber implicit? (vip_vab == impl)

Examples: mpdata_rhs_vip_common.hpp mpdata_rhs_vip_prs_common.hpp boussinesq_impl.hpp supercell.hpp (from mwarusz/supercell branch)

[pdziekan]

No need to change anything. boussinesq_impl and supercell model perturbation of tht and not full tht, therefore there's a tht forcing proportional to w*grad(tht_e) which, together with trapezoidal integration of absorbers, leads to vip normalization.