add expokit.f90