This dataset is used to understand the pressure drag of shallow cumulus clouds, the relationship between the pressure drag of individual clouds and the cloud ensemble. We performed large eddy simulations of shallow cumulus clouds based on the Barbados Oceanographic and Meteorological Experiment (BOMEX) using the Met Office-NERC (National Environment Research Council) Cloud model (MONC). The grid spacing is 25 m in all directions and the domain size is (15X15X3) km^3. As the whole simulation has a huge dataset, we present a subset of the results that can be used to reproduce the figures in 'Pressure drag for shallow cumulus clouds: from thermals to cloud ensemble' (submitted to Geophysical Research Letters; 2020). Access to the all simulation data can be provided on request. The full dataset will be deposited in the CEDA Archive in due course. The dataset contains four sub-directories with generated data and another sub-directory with scripts for plotting the figures in the article. The generated data include: [1]. The conditionally averaged budget terms in vertical velocity equation for different parts of cloud ensemble; [2]. The conditionally averaged budget terms in vertical velocity equation for all tracked clouds; [3]. The conditionally averaged budget terms in vertical velocity equation for tracked clouds with different lifetime; [4]. The conditionally averaged budget terms in vertical velocity equation for each individually tracked cloud over its lifetime; [5]. The vertical cross section data of cloud structure; [6]. The positions of tracked clouds; [7]. The composited structure of clouds at 600m, 1000m and 1800m.