C units: length in m; energy,momentum,mass in GeV; angles in radians
      common /cophys/pi,me,alpha,hbarc
      real me,alpha,hbarc
      double precision pi
      integer nsites
      parameter (nsites=8)
      common /cotarg/Z,a,radlen,Aphonon,mospread,betaFF,ucell(3,nsites)
      real Z,a,radlen,Aphonon,mospread,betaFF,ucell
      common /cosetup/thx,thy,rotate(3,3),E,emit,spot,D,t,collim
      double precision thx,thy,rotate
      real E,emit,spot,D,t,collim