I found a different data source that seems to work well with @Patrick's original solution. I skipped the subset portion and just used the profile function, because it doesn't hurt me to download a larger time frame and that way I don't have to worry about a bounding box.

This is my new source: https://opendap.cr.usgs.gov/opendap/hyrax/MOD13C1.061/MOD13C1.061.ncml