Bug in rain 2M+PEDMF terminal velocity calculation

[haakon-e]

In the updrafts, we pass ᶜρʲ * ᶜn_raiʲ to CM2.rain_terminal_velocity, see e.g.:

ClimaAtmos.jl/src/cache/prognostic_edmf_precomputed_quantities.jl

Lines 750 to 759 in 13875ac

	@. ᶜwᵣʲs.:($$j) = getindex(
	CM2.rain_terminal_velocity(
	cm2p.sb,
	cm2p.rtv,
	max(zero(Y.c.ρ), Y.c.sgsʲs.:($$j).q_rai),
	ᶜρʲs.:($$j),
	max(zero(Y.c.ρ), ᶜρʲs.:($$j) * Y.c.sgsʲs.:($$j).n_rai),
	),
	2,
	)

while in the environment, we only pass ᶜn_rai⁰ to CM2.rain_terminal_velocity, see e.g.:

ClimaAtmos.jl/src/cache/precipitation_precomputed_quantities.jl

Lines 444 to 454 in 13875ac

	# Rain (mass)
	@. ᶜw⁰ = getindex(
	CM2.rain_terminal_velocity(
	cm2p.sb,
	cm2p.rtv,
	max(zero(Y.c.ρ), ᶜq_rai⁰),
	ᶜρ⁰,
	max(zero(Y.c.ρ), ᶜn_rai⁰),
	),
	2,
	)

This applies to all condensate species.

The correct approach is to pass ρ * n