Added Rutger's anisitropy code