A minor problem with this: vertical advection is the only our regression that use temporaries in cuda; our perftests use power of two domain sizes. Hence we need yet another test to track the performance improvement of that change.