  We describe Zhu recursion for a vertex operator algebra (VOA) and its modules
on a genus $g$ Riemann surface in the Schottky uniformisation. We show that
$n$-point (intertwiner) correlation functions are written as linear
combinations of $(n-1)$-point functions with universal coefficients given by
derivatives of the differential of the third kind, the Bers quasiform and
certain holomorphic forms. We use this formula to describe conformal Ward
identities framed in terms of a canonical differential operator which acts with
respect to the Schottky moduli and to the insertion points of the $n$-point
function. We consider the generalised Heisenberg VOA and determine all its
correlation functions by Zhu recursion. We also use Zhu recursion to derive
linear partial differential equations for the Heisenberg VOA partition function
and various structures such as the bidifferential of the second kind,
holomorphic $1$-forms, the prime form and the period matrix. Finally, we
compute the genus $g$ partition function for any rational Euclidean lattice
generalised VOA.
