make long wave orbital velocity a function of z, and evaluate it at z = eta in the integrand calculation
make long wave orbital velocity a function of z, and evaluate it at z = eta in the integrand calculation

milancurcic/stokes-doppler,milancurcic/stokes-doppler