Earth System Models (ESM) are essential tools to investigate projected changes in precipitation patterns due to different climate scenarios. However, the atmospheric component of traditional ESMs, an atmospheric Global Circulation Model (GCM) with resolution of around 100 km, cannot directly simulate precipitation-generating convective processes, as they occur on scales of a few kilometers, so on much smaller length scales than the grid resolution (Bony et al., 2015;Randall et al., 2003). Therefore, GCMs rely on parameterizations to represent the effect of convective sub-grid-scale processes on the large-scale resolved state (Randall, 2013;Randall et al., 2003). However, the models exhibit large persisting systematic biases such as the presence of a Double Inter Tropical Convergence Zone (ITCZ, a zonal band of strong precipitation in the tropics that forms the ascending