anchor.diff:           real(kind=fPrec)            p0
            real(kind=fPrec)            zlm
            real(kind=fPrec)            x,xp
  -         real(kind=fPrec)            shift
  -         real(kind=fPrec)            x_shift, xp_shift,s_shift !coordinates after shift/rotation
  -         real(kind=fPrec)            x_rot, xp_rot,s_rot
  -         real(kind=fPrec)            x_temp, xp_temp,s_temp !!all the _temp variables are used when you hit the cry from below
  -         real(kind=fPrec)            tilt_int, x_int,xp_int,s_int       !all the _int variables are used when you hit the cry from below (int=interaction point)
            real(kind=fPrec)            x00
            real(kind=fPrec)            z
            real(kind=fPrec)            z00
  
positive.code:          real(kind=fPrec)            p
           real(kind=fPrec)            dpop
           real(kind=fPrec)            s
           real(kind=fPrec)            a_eq,b_eq,c_eq,Delta
           real(kind=fPrec)            ENOM
  !
           real(kind=fPrec)            x_PRINT,xp_PRINT,y_PRINT,yp_PRINT
  
negative.code: 
    integer, parameter :: max_ncoll = 99
  
    logical,           save :: bool_create
  
    integer,           save :: iProc
    integer,           save :: n_chan
    integer,           save :: n_VR
    integer,           save :: n_amorphous
    logical,           save :: changed_tilt1(max_ncoll)
    logical,           save :: changed_tilt2(max_ncoll)
  
  