I think I know what went wrong for the SO2 binning. You probably got the following error message
ERROR: variables don't have the same number of dimensions (SO2_column_number_density_weight)
This is a bug in HARP that I just fixed. You can work around it for the moment by removing the SO2_column_number_density_weight
variable before binning (this will create slightly different values for the daily/monthly averages, but this should not be significant):
product = harp.import_product(r"/Users/jpoetzscher/Downloads/SO2 OCT 1/*.nc",
operations='SO2_column_number_density_validity>50;keep(latitude_bounds,longitude_bounds,SO2_column_number_density);bin_spatial(1801,-90,0.1,3601,-180,0.1);derive(SO2_column_number_density [DU]);exclude(SO2_column_number_density_weight)',
post_operations='bin();squash(time, (latitude_bounds,longitude_bounds));derive(latitude {latitude});derive(longitude {longitude});exclude(latitude_bounds,longitude_bounds,latitude_bounds_weight,longitude_bounds_weight,count,weight)')