def JRC_SDF(water_depth):
    """Replication of the SDF provided by the JRC for South America"""

    dmg = (1.006049 + (0.001284178 - 1.006049) /
           (1 + (water_depth / 797962.3) ** 0.9707786) ** 681383)
    return(dmg)