Code based on venus-levis for computation of delta-f. Basic venus-levis routines are kept (e.g. reading of equilibrium, guiding-centre pushing), but higher-level modules have been removed and substituted by new ones.