def TOC_schmoker1979(rho_b, A=157.0, B=58.1):
    """
    Schmoker (1979) method of TOC caluculation from bulk density to estimate 
    TOC in devonian shales.
    
    bulk density units: g/cc
    """
    TOC = (A / rho_b) - B
    return TOC