Gradient fusion

I am using Gauss-legendre with p=5 and quad. Anti-aliasing is disabled. I am doing packing since I want to use baseflow gradients in the future time steps and this baseflow gradient graph is only run once when starting. This is following our previous discussion in Save view matrices - #9 by Zhenyang

Regards,
Zhenyang