def test_distribute_lwp_to_liquid_clouds():
lwc = np.array([[1, 1, 1],
[2, 2, 2]])
lwp = np.array([2, 8])
result = atmos.distribute_lwp_to_liquid_clouds(lwc, lwp)
correct = [[2/3, 2/3, 2/3], [16/6, 16/6, 16/6]]
assert_array_equal(result, correct)
def _adiabatic_lwc_to_lwc(self) -> np.ndarray:
"""Initialises liquid water content (g/m3).
Calculates LWC for ALL profiles (rain, lwp > theoretical, etc.),
"""
lwc_scaled = atmos.distribute_lwp_to_liquid_clouds(self.lwc_adiabatic, self.lwc_source.lwp)
return lwc_scaled / self.dheight