Updated my branch to be based on master

Hi Adam, thanks for catching this bug and fixing it. Your fix looks good and works well, and I appreciate the updates to make it more readable as well.

One question I have is for the addDefaultPlotVars.m, is there a point to just adding the plot variables to the last dataset? I don't see an instance where the user would pass in multiple hydro structs, because the readX.m scripts already apply the default plotting variables. Your implementation here works fine with just 1 so no issues with it, I'm just curious.

Hi @jtgrasb Thanks for the review! On addDefaultPlotVars -- This will break when hydro is an array of structures. For example if I wanted to read multiple output files to compare or combine them:

hydro = readCapytaine(hydro, 'file1.nc', ...)
hydro = readCapytaine(hydro, 'file2.nc', ...)
hydro = combineBEM(hydro)
plotBEMIO(hydro)

Then the second readCapytaine call will fail because addDefaultPlotVars tries to add plotDofs, etc to an array instead of a structure. I was doing this because I was processing a lot of hydro files at once and comparing them.

I think we previously had a set-up like this in compareBEMIO, but it looks like it was changed to save computational time:

@akeeste Thanks for the explanation, makes sense. Everything looks good then, I will merge this today.