def cs_water_Tp(T=298, p=1):
    """
    Calculates chemical shift of water based on Temparature and pressure (according to Maciej).
    """
#   according to Maciej
    cs = 5.945e-7*p*T  -1.612e-4*p  -1.025e-2*T  + 7.866
    return(cs)