These data are a realistic simulation of the radiation
damage situation faced with a lysozyme-sized protein
growing ~5 micron crystals and shot with a 6 micron beam.

The exposure time was adjusted to get decent resolution
on the first image, but unfortunately, you don't get 
very many shots before the crystal dies!  And once you
start trying to scale and merge the data, it gets even
worse.  The rad dam creates "non-isomorphism" that
rapidly becomes unmanageable, despite the fact that
the damage model used here is actually a very simple
equation (described by Holton and Frankel, Acta D 2010).

The real trick with this dataset, however, is the
fact that the a and b axes are the same length, but
the space group is P212121.  This means that autoindexing
will get a and b swapped for about half of the wedges 
and you will need to check each one of them and "flip"
the ones that don't agree with the rest.  Is this a common
problem when mergeing data from many crystals?  Yes.
Is there a program for doing this automatically?  No.

Good luck.