Modelling of compressional alfvén eigencmodes in axisymmetric toroidal fusion plasmas