anchor.diff:  
     double precision nm, np , stress(1:4)
     double precision  XFN(1:3,1:nbNodes) , Selements(1:4,1:nbElements)
  -  double precision A(1:3)  , B(1:3),C(1:3),FA(1:4)   , mag
  -  double precision rays(1:9,1:numberOfRays), P(1:3),N1(1:3), nearestPoint(1:3), normal(1:3) , nearestNormal(1:3)
  -  double precision e(1:3) , d(1:3),r(1:3), hither,  yon, t , rAngle , iAngle
  +  double precision A(1:3)  , B(1:3),C(1:3), mag
  +  double precision rays(1:9,1:numberOfRays), P(1:3), nearestPoint(1:3), normal(1:3) , nearestNormal(1:3)
  +  double precision e(1:3) , d(1:3), hither,  yon, t
     logical result 
  -  integer elmnew(1:3,1:nbElements), i  , iRays , iReflection, nearestElement , s , k
  +  integer elmnew(1:3,1:nbElements), i  , iRays , iReflection, nearestElement
     double precision nearestIntersection
     double precision omega_zero, epsilono ,z0, RADX, disp,opticalPower, lightSpeed
     double precision x , y , z, my_cap_center(1:3), H
  
positive.code: 
    double precision N(1:3,1:nbElements) , zc, PI
  
    double precision A(1:3),B(1:3),C(1:3),center(1:3),tail(1:3) , t
    double precision XFN(1:3,1:nbNodes)
    integer elmnew(1:3,1:nbElements)
  
  
negative.code:   integer elmnew(1:3,1:nbElements), ii
    double precision  XFN(1:3,1:nbNodes) , Fnodes(1:3,1:nbNodes)
  
    double precision  shpfs(1:7,1:nbElements)
  
    double precision Felements(1:3,1:nbElements),Selements(1:4,1:nbElements)
    double precision z
    double precision total, Fphi, FOSTAR
    REAL (8) FX_Back,FZ_Back, FX_FRONT,FZ_FRONT, zc
    double precision rays(1:9,1:numberOfRays)
    integer i,j ,k
  
  