calIsenFlux1d.m

function flux1d  = calIsenFlux1d(vcomp1d,PT1d,pfull,isenCoordEdge,pInterp)
% unit is not converted into standard using G

vcomp1dInterp = interp1(pfull,vcomp1d,    pInterp,'spline','extrap');
PT1dInterp    = interp1(pfull,PT1d,       pInterp,'spline','extrap');

[~, bin] = histc(PT1dInterp,isenCoordEdge);

flux1d = zeros(length(isenCoordEdge)-1,1);

for i = 1:length(bin)
    if bin(i) ~= 0
        flux1d(bin(i)) = flux1d(bin(i))+vcomp1dInterp(i);
    end
end