function h = Bg_gw(p, BG, Ppr)
h = interptable(Ppr, BG, p);
end