Should be final version - very slow! probably ca. 6 h to run - could do with parallization, but won't bother now