poor resolution due to output compression
The data range of the 2D output variable total column integrated cloud water (tqc_dia
) has been discretised to only a few values, as can be seen via this MWE:
import eurec4a
import pandas as pd
import matplotlib.pyplot as plt
cat = eurec4a.get_intake_catalog()
timeslice = slice(pd.to_datetime('2020-01-11-10:00'),pd.to_datetime('2020-01-11-12:00'))
data = cat['simulations']['ICON']['LES_CampaignDomain_control']['surface_DOM01'].to_dask().sel(time=timeslice)
plt.hist(data['tqc_dia'].values.ravel(), range=(.01, .1), bins=100)
# or :
# plt.hist(np.log(data['tqc_dia'].values.ravel()), range=(-12, 2), bins=100)
plt.show()
Edited by Hernan Campos