I have a large netcdf file with several variables. I need to do discrete integration along a dimension to a variable say temperature of shape (80,100,300000) with dimensions (time, depth, nodes). So, I tried with dividing the large dataset into chunks with xarray and then tried to apply the function scipy.integrate.simps, but failed. I have a large netcdf file with several variabl