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