def pres_from_hybrid(psfc, hya, hyb, p0=100000.):
    """Return pressure field on hybrid-sigma coordinates, 
       assuming formula is 
       p = a(k)*p0 + b(k)*ps.
    """
    return hya*p0 + hyb*psfc