code: def density_water(temp):
      
      ut.check_range([temp, ">0", "Temperature in Kelvin"])
      rhointerpolated = interpolate.CubicSpline(WATER_DENSITY_TABLE[0],
                                                      WATER_DENSITY_TABLE[1])
      return np.asscalar(rhointerpolated(temp))
docs: Return the density of water at a given temperature.
  
      If given units, the function will automatically convert to Kelvin.
      If not given units, the function will assume Kelvin.