Routines for surface and volume averaging

Translation from the corresponding IDL libraries
Small correction in order to make it work with stop animation feature
update to elm-pb manual

More verbose algebra included. Still some issues to work out. Please
check for errors/typos
post processing python routines