Added shear calculation to script
Added removal of empty interface cells for shear plotting in f_contactline.py