Если я правильно понимаю ваш вопрос, у вас есть массивы dask SX и VX с формами ( i , j ) и ( l , m , n ) соответственно, и требуется вычислить своего рода "внешнюю" разницу между ними, в результате чего получается массив DX формы ( i , j , l , m , n )?В этом случае, то, что вы ищете, это
DX = SX[:,:,None,None,None] - VX[None,None,:,:,:]
Это приведет SX и VX к формам ( i , j , 1,1,1) и (1,1, l , m , n ) соответственно, что позволило бывычитание для правильной трансляции по размерам.Получившийся массив dask DX будет иметь форму ( i , j , l , m , n ) и будет ленивым по своей природе (значения не будут фактически рассчитываться до тех пор, пока они не понадобятся).