function [JV] = CalcJV(Wvl,BinWidth,FluxSingle,XSQY);

[nRxn,nWvl] = size(XSQY); 
JV = zeros(1,nRxn);
for iRxn = 1:nRxn
    
    jvFrac = zeros(1,nWvl);
    
    for iWvl = 1:nWvl
        jvFrac(iWvl) = XSQY(iRxn,iWvl)*FluxSingle(iWvl)*BinWidth(iWvl);
    end
    
     
    
    
    
    JV(iRxn) = sum(jvFrac);
    
    
    
   
end