Before change to bbar based code
Added numerical estimation of stress