MATLAB: How to calculate total precipitation per day using hourly data ? (netcdf)

I have hourly data from ECMWF ERA5 for each day in a specific year. I want to convert that data from hourly to daily. Copernicus has a Python code for this here https://confluence.ecmwf.int/display/CKB/ERA5%3A+How+to+calculate+daily+total+precipitation.
I want to know what is the matlab code to do this? I was upload the netcdf file in my google drive here: https://drive.google.com/open?id=1qm5AGj5zRC3ifD1_V-ne2nDT1ch_Khik time steps of each day are:

0:00
1:00
2:00
3:00
4:00
5:00
6:00
7:00
8:00
9:00
10:00
11:00
12:00
13:00
14:00
15:00
16:00
17:00
18:00
19:00
20:00
21:00
22:00
23:00

Notice to cover total precipitation for 1st January 2017 for example, we need two days of data: 1st January 2017 time = 01 - 23 will give you total precipitation data to cover 00 - 23 UTC for 1st January 2017 2nd January 2017 time = 00 will give you total precipitation data to cover 23 - 24 UTC for 1st January 2017 here is ncdisp():

>> ncdisp(filename)
Source:
           C:\Users\Behzad\Desktop\download.nc
Format:
           64bit
Global Attributes:
           Conventions = 'CF-1.6'
           history     = '2019-11-01 07:36:15 GMT by grib_to_netcdf-2.14.0: /opt/ecmwf/eccodes/bin/grib_to_netcdf -o /cache/data6/adaptor.mars.internal-1572593007.3569295-19224-27-449cad76-bcd6-4cfa-9767-8a3c1219c0bb.nc /cache/tmp/449cad76-bcd6-4cfa-9767-8a3c1219c0bb-adaptor.mars.internal-1572593007.35751-19224-4-tmp.grib'
Dimensions:
           longitude = 49
           latitude  = 41
           time      = 8760
Variables:
    longitude
           Size:       49x1
           Dimensions: longitude
           Datatype:   single
           Attributes:
                       units     = 'degrees_east'
                       long_name = 'longitude'
    latitude 
           Size:       41x1
           Dimensions: latitude
           Datatype:   single
           Attributes:
                       units     = 'degrees_north'
                       long_name = 'latitude'
    time     
           Size:       8760x1
           Dimensions: time
           Datatype:   int32
           Attributes:
                       units     = 'hours since 1900-01-01 00:00:00.0'
                       long_name = 'time'
                       calendar  = 'gregorian'
    tp       
           Size:       49x41x8760
           Dimensions: longitude,latitude,time
           Datatype:   int16
           Attributes:
                       scale_factor  = 3.0792e-07
                       add_offset    = 0.010089
                       _FillValue    = -32767
                       missing_value = -32767
                       units         = 'm'
                       long_name     = 'Total precipitation'

tp is my variable which have 3 dimensions (lon*lat*time) = 49*41*8760 I want it in the 49*41*365 for a non-leap year.
The result should be the daily values for the whole year.

While some vectorized versions may exist that reshape your vector into 4 dimensions, a simple for loop will do the job.

tp_daily=zeros(size(tp,1),size(tp,2),365);
for ii=0:364
    day=tp(:,:,ii*24+1:(ii+1)*24); %grab an entire day 
    tp_daily(:,:,ii+1)=sum(day,3); % add the third dimension
end