  We briefly discuss some issues concerning including the pseudoscalar glueball
interpolating operators into the variational basis for computing the masses of
the eta and etaprime mesons. As a start we present some preliminary results for
the correlators from pseudoscalar glueball operators on nf = 2+1+1 twisted mass
gauge configurations. Preliminary results for the effect of open boundary
conditions on the masses of three glueballs computed in quenched QCD are
presented. The statistics of glueball correlators are briefly discussed.
