  A new next-to-leading order Monte Carlo program for the calculation of fully
differential heavy quark cross sections in electroproduction is described. A
comparison between the theoretical predictions and the latest charm production
data from H1 and ZEUS at HERA is presented.
