If you know already at which k-point the band gap is, and esp. if it's direct, then it's fairly easy. For instance, assuming it's direct at the Gamma point, you can do
import sys
b = nlread(sys.argv[1], BulkConfiguration)[0]
t = Bandstructure(b, route=['G','G'], points_per_segment=2)
vbmax_ix = numpy.where(t.evaluate()[0]<=0.*eV)[0][-1]
cbmin_ix = numpy.where(t.evaluate()[0]>=0.*eV)[0][0]
bandgap = t.evaluate()[0][cbmin_ix]-t.evaluate()[0][vbmax_ix]
print "Direct band gap at Gamma point", bandgap
Save this as "bandgap.py", then run as "atkpython bandgap.py file.nc" where file.nc contains the converged BulkConfiguration.