code: def dens_floc(ConcAl, ConcClay, DIM_FRACTAL, DiamTarget, coag, material, Temp):
      
      WaterDensity = pc.density_water(Temp).magnitude
      return ((dens_floc_init(ConcAl, ConcClay, coag, material).magnitude
               - WaterDensity
               )
              * (material.Diameter / DiamTarget)**(3 - DIM_FRACTAL)
              + WaterDensity
              )
docs: Calculate floc density as a function of size.