anchor.diff:  ! parameter(fnavo=6.02214129e23_fPrec)
     real(kind=fPrec), save :: cgen(200,nmat)
     character(4), save :: mname(nmat)
  -
  -!  common/mater/anuc(nmat),zatom(nmat),rho(nmat),emr(nmat)
  -!  common/coul/tlcut,hcut(nmat),cgen(200,nmat),mcurr
  -!  common/scat/cs(0:5,nmat),csref(0:5,nmat),bnref(nmat),freep(nmat)
  -!  common/scatu/cprob(0:5,nmat),bn(nmat),bpp,xln15s,ecmsq
  -!  common/scatu2/xintl(nmat),radl(nmat),mname
  -!  common/scatpp/pptot,ppel,ppsd
  -!  common/sppref/pptref,pperef,pref,pptco,ppeco,sdcoe,freeco
  -!! real(kind=fPrec) exenergy
  -!! common/meanexen/exenergy(nmat)
  -!  common/cmcs1/zlm1
  -!  common/sindif/xpsd,zpsd,psd
  -!  common/cdpodx/dpodx
  -!  common/cions/edens(nmat),pleng(nmat)
  -
  -
  -!  common/materia/mat
  -!  common/phase/x,xp,z,zp,dpop
  -!  common/nommom/p0
  -!  common/cjaw1/zlm
   
   !>
   !! block data scdata
  
positive.code:   real(kind=fPrec), allocatable, save :: tdispx(:)  !(nblz)
    real(kind=fPrec), allocatable, save :: tdispy(:)  !(nblz)
  
  !  common /rtwiss/ tbetax,tbetay,talphax,talphay,torbx,torbxp,torby,torbyp,tdispx,tdispy
  !
  !-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-
  !
  ! Variables for finding the collimator with the smallest gap
  
negative.code:   real(kind=fPrec), save :: beamsize1, beamsize2,betax1,betax2,betay1,betay2, alphax1, alphax2,alphay1,alphay2,minAmpl
  !SEPT2005
  
    integer, private, save :: nev
  !  real(kind=fPrec), private, save :: c_xmin,c_xmax,c_xpmin,c_xpmax,c_zmin,c_zmax,c_zpmin,c_zpmax,length
  !
  !  real(kind=fPrec), private, save :: c_mybetax,c_mybetaz,mymux,mymuz,atdi
  
    real(kind=fPrec), private, save :: length
  
  