computer code is working (though h5py mpio didn't work too well).

Some thought needs to be done on the physics.
For instance,
* we have the two phasepols, but they are different.
Which one should we use.
* it could happen that we average both fh's in time and frequency
  which would reduce the resolution of the data. A warning is printed
  but no code has been implemented for interpolation.

* I need to double-check the order of the frequency channels