some conversion to single precision. Now I have trouble with rvavg being a vector of NaNs in some locations. This most likely has something to do with the shift to interp1q in interpolate.m . I might change back but there is code to handle that corner case now.
updates to interpolate.m - interp1q, can handle multiple timesteps at one go