Generating reference evapotranspiration surfaces from the Hargreaves equation at watershed scale