HARP does not support direct ingestion of MODIS data. But you can try to read the data using some other library in Python and then construct a HARP Product + Variables in Python.
You can then pass all products of a day to a harp.execute_operations invocation and have this perform a spatial bin and merge to produce a daily average.