anchor.diff:  ! M. Fitterer, FNAL, A. Mereghtti, CERN
  -! last modified: 11-04-2019
  -! Common block for electron lens definition
  +! last modified: 25-04-2019
  +! Common block for electron lens
   module elens
   
     use parpro
  
positive.code:         ! UNIFORM: eLens with uniform profile
          ! formula: (r^2-r1^2)/(r2^2-r1^2)
  !FOX    FRR=((RR+ELE_R1)*(RR-ELE_R1))/ELENOR ;
           
        case (2)
          ! GAUSSIAN: eLens with Gaussian profile
          ! formula: (exp(-r1^2/2sig^2)-exp(-r^2/2sig^2))/(exp(-r1^2/2sig^2)-exp(-r2^2/2sig^2))
  
negative.code: !FOX    LINLIN=TMPCOF ;
  !FOX    TMPTMP=RR;
          do kk=2,mPoints
            TMPCOF=COF(kk)
            RKK=real(kk-1,fPrec)
  !FOX      LINLIN=LINLIN+(TMPTMP*TMPCOF)/RKK ;
  !FOX      TMPTMP=TMPTMP*RR ;
          end do
  !FOX    FRR=(LINLIN-ELEFR1)/ELENOR ;
          call dealloc(cof,'cof')
  