Two things:

1. Remove the totally unnecessary loop iterating through the vertical levels on the vertical interpolation secion. No idea what that was doing there.

2. Change the interpolation of the vertical profiles to scale the POLCOMS depths to the current FVCOM node's depth range. This means the full profile structure is maintained in the interpolated profiles but that it is squashed or stretched to fit the scaled depths. This also means the surface and seabed value are fully captured.
Big changes here:

1. Do the vertical interpolation of the u and v components separately then at the end calculate the depth averaged velocity as well as the depth averaged u and v components into separate fields in Mobj. This is all necessary for outputting the ASCII files for mean flow and the open boundaries.

2. Make the vertical interpolation scale the POLCOMS-ERSEM depth range to match the current FVCOM node's depth range, thus squeezing or stretching the vertical velocity profiles into the FVCOM depth.

3. Add a commented out section which can be used to plot the vertical profiles to make sure they have interpolated correctly.