anchor.diff:  
     double precision aa, bb, cc, dd, ee, ff, gg, hh, ii, jj, kk, ll, t, beta, gamma, denom;
     double precision eihf, gfdi, dheg, akjb, jcal, blkc;
  -  double precision P(1:3), N1(1:3), normal(1:3)
  +  double precision P(1:3), normal(1:3)
     double precision A(1:3) , B(1:3) , C(1:3)
     double precision e(1:3), d(1:3)
     double precision hither,  yon
  
positive.code: 
    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
    logical result 
    integer elmnew(1:3,1:nbElements), i  , iRays , iReflection, nearestElement , s , k
    double precision nearestIntersection
    double precision omega_zero, epsilono ,z0, RADX, disp,opticalPower, lightSpeed
    double precision x , y , z, my_cap_center(1:3), H
  
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
  
  