function h = muw_ow(p, vwi, prw, cvw)
h = vwi + cvw * (p - prw);