We present a numerical differentiation technique for the calculation of directional derivatives of electronic wave functions evaluated using a multireference configuration interaction method. The technique has been implemented in the set of programs MELD, and its application is illustrated by calculations of nonadiabatic couplings for the ArH2+ quasimolecule.