f: left side numerical reconstruction at (i+1/2) interface
function wcL(v1,v2,v3,v4,v5)
eps = 1.0e-6