function h = Bw_ow(p, Bwi, prw, cw)
h = Bwi * (1 - cw * (p - prw));