Add linear component 
loop: 1 <= j <= N 
 Calculate Mlocal
 loop: |s| <= Mlocal
 if (mode J(j,s) exists) then 
 Add linear contribution (including wave dissipation) 
 end 
 end 
end 