Error when trying to fit data to gamma distribution

Here is my code:

data("airquality")
head(airquality)

Solar <- airquality$Solar.R

library(fitdistrplus)
fit.solar <- fitdist(Solar, distr = "gamma", method = "mle")

Any thoughts? Thank you!