function g = MultiD_times(f,a)
[Nx,Nv,Nz]= size(f);
f = reshape(f,Nx*Nv, Nz);
g = f*diag(a);
g = reshape(g,Nx,Nv,Nz);
