anchor.diff:    character(len=16), save :: castordir
     character(len=80), save :: filename_dis
   
  -!  common /grd/ myenom,mynex,mdex,myney,mdey,                        &
  -!  &nsig_tcp3,nsig_tcsg3,nsig_tcsm3,nsig_tcla3,                       &
  -!  &nsig_tcp7,nsig_tcsg7,nsig_tcsm7,nsig_tcla7,nsig_tclp,nsig_tcli,   &
  -!  &nsig_tcth1,nsig_tcth2,nsig_tcth5,nsig_tcth8,                      &
  -!  &nsig_tctv1,nsig_tctv2,nsig_tctv5,nsig_tctv8,                      &
  -!  &nsig_tcdq,nsig_tcstcdq,nsig_tdi,nsig_tcxrp,nsig_tcryo,            &
  -!  &smin_slices,smax_slices,recenter1,recenter2,                      &
  -!  &fit1_1,fit1_2,fit1_3,fit1_4,fit1_5,fit1_6,ssf1,                   &
  -!  &fit2_1,fit2_2,fit2_3,fit2_4,fit2_5,fit2_6,ssf2,                   &
  -!  &emitnx0_dist,emitny0_dist,emitnx0_collgap,emitny0_collgap,        &
  -!  &xbeat,xbeatphase,ybeat,ybeatphase,                                &
  -!  &c_rmstilt_prim,c_rmstilt_sec,c_systilt_prim,c_systilt_sec,        &
  -!  &c_rmsoffset_prim,c_rmsoffset_sec,c_sysoffset_prim,                &
  -!  &c_sysoffset_sec,c_rmserror_gap,nr,                                &
  -!  &ndr,driftsx,driftsy,pencil_offset,pencil_rmsx,pencil_rmsy,        &
  -!  &sigsecut3,sigsecut2,enerror,                                      &
  -!  &bunchlength,coll_db,name_sel,                                     &
  -!  &castordir,filename_dis,nloop,rnd_seed,c_offsettilt_seed,          &
  -!  &ibeam,jobnumber,do_thisdis,n_slices,pencil_distr,                 &
  -!  &do_coll,                                                          &
  -!  &do_select,do_nominal,dowrite_dist,do_oneside,dowrite_impact,      &
  -!  &dowrite_secondary,dowrite_amplitude,radial,systilt_antisymm,      &
  -!  &dowritetracks,cern,cdb_doNSig,do_mingap
  -!+cd info
     integer, save :: ie, iturn, nabs_total
  -!  common  /info/ ie,iturn,nabs_total
   
    !-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-
   !+cd dbcommon
  
positive.code:   integer, save :: pencil_distr
  
    real(kind=fPrec), save :: myenom,mynex,mdex,myney,mdey,                    &
    &nsig_tcp3,nsig_tcsg3,nsig_tcsm3,nsig_tcla3,                       &
    &nsig_tcp7,nsig_tcsg7,nsig_tcsm7,nsig_tcla7,nsig_tclp,nsig_tcli,   &
    &nsig_tcth1,nsig_tcth2,nsig_tcth5,nsig_tcth8,                      &
    &nsig_tctv1,nsig_tctv2,nsig_tctv5,nsig_tctv8,                      &
    &nsig_tcdq,nsig_tcstcdq,nsig_tdi,nsig_tcxrp,nsig_tcryo,            &
  !SEPT2005 add these lines for the slicing procedure
    &smin_slices,smax_slices,recenter1,recenter2,                      &
    &fit1_1,fit1_2,fit1_3,fit1_4,fit1_5,fit1_6,ssf1,                   &
  
negative.code:   &driftsx,driftsy,pencil_offset,pencil_rmsx,pencil_rmsy,            &
    &sigsecut3,sigsecut2,enerror,bunchlength
  
    real(kind=fPrec), private, save :: nsig_tcp3    = zero
    real(kind=fPrec), private, save :: nsig_tcsg3   = zero
    real(kind=fPrec), private, save :: nsig_tcsm3   = zero
    real(kind=fPrec), private, save :: nsig_tcla3   = zero
    real(kind=fPrec), private, save :: nsig_tcp7    = zero
    real(kind=fPrec), private, save :: nsig_tcsg7   = zero
    real(kind=fPrec), private, save :: nsig_tcsm7   = zero
    real(kind=fPrec), private, save :: nsig_tcla7   = zero
    real(kind=fPrec), private, save :: nsig_tclp    = zero
    real(kind=fPrec), private, save :: nsig_tcli    = zero
    real(kind=fPrec), private, save :: nsig_tcth1   = zero
    real(kind=fPrec), private, save :: nsig_tcth2   = zero
    real(kind=fPrec), private, save :: nsig_tcth5   = zero
    real(kind=fPrec), private, save :: nsig_tcth8   = zero
    real(kind=fPrec), private, save :: nsig_tctv1   = zero
    real(kind=fPrec), private, save :: nsig_tctv2   = zero
    real(kind=fPrec), private, save :: nsig_tctv5   = zero
    real(kind=fPrec), private, save :: nsig_tctv8   = zero
    real(kind=fPrec), private, save :: nsig_tcdq    = zero
    real(kind=fPrec), private, save :: nsig_tcstcdq = zero
    real(kind=fPrec), private, save :: nsig_tdi     = zero
    real(kind=fPrec), private, save :: nsig_tcxrp   = zero
    real(kind=fPrec), private, save :: nsig_tcryo   = zero
  
    ! From collimation_comnul
    real(kind=fPrec), public,  save :: emitnx0_dist    = zero
    real(kind=fPrec), public,  save :: emitny0_dist    = zero
    real(kind=fPrec), public,  save :: emitnx0_collgap = zero
    real(kind=fPrec), public,  save :: emitny0_collgap = zero
  
    real(kind=fPrec), private, save :: nr
  
  