real function xsgama(energy)
      real energy
      write(6,*) 'this should never be called!'
      end