Originally noted in #753

In compute_u_with_rbf there is a loop around ens_size, for get_reconstruct but uval(1:ens_size) is set to set to the single value ureconstructzonal/ureconstructmeridional so you are getting the whole ensemble set each step of the loop

whole ensemble == last value from the loop.

DART/models/mpas_atm/model_mod.f90

Lines 6214 to 6226 in e218864

	do e = 1, ens_size
	call get_reconstruct(nedges, lat*deg2rad, lon*deg2rad, &
	coeffs_reconstruct, on_a_sphere, veldata(:, e), &
	ureconstructx, ureconstructy, ureconstructz, &
	ureconstructzonal, ureconstructmeridional)
	if (zonal) then
	uval = ureconstructzonal
	else
	uval = ureconstructmeridional
	endif
	enddo

Originally posted by @hkershaw-brown in #753