mid-pressure calculation in gsw_distance.m

[chrisdane]

Hi

I was wondering if

GSW-Matlab/Toolbox/gsw_distance.m

Line 128 in 08a4b3a

p_mid = 0.5*(p(:,1:nla-1) + p(:,1:nla-1));

is correct since the same index 1:nla-1 is used to calculate the average pressure of different locations. So I think

>> p = [200 1000];
>> p_mid = 0.5*(p(:,1:nla-1) + p(:,1:nla-1))

p_mid =

   200

is wrong, as it should be rather

>> p_mid = 0.5*(p(:,1:nla-1) + p(:,2:nla))  

p_mid =

   600

as it is also used for latitude:

GSW-Matlab/Toolbox/gsw_distance.m

Line 129 in 08a4b3a

lat_mid = 0.5*(lat(:,1:nla-1) + lat(:,2:nla));

[PaulMBarker]

I have corrected this and the change will appear in version 3.06.12