anchor.diff:                    x(i,2) = 65.5
                     x(i,3) = 0.968
                 elseif (plate(i).eq.'ITRF') then
  -                  x(i,1) = -68.2+180.0
  -                  x(i,2) = -63.4
  +                  x(i,1) = -68.2
  +                  x(i,2) = 63.4
                     x(i,3) = 0.677
                 elseif (plate(i).eq.'PA') then
                     x(i,1) = 0.0
  
positive.code: 
        ! Compute Cartesian Euler pole
        if (i.gt.0) then
            if (verbos.eq.1) then
                write(0,*) 'getpole says: Using plate circuit model'
            endif
            call model2pole(pole,polec,model,verbos)
        elseif (i.lt.0) then
  
negative.code: 
        SUBROUTINE getpole(pole,polec,model,verbos)
        IMPLICIT none
        CHARACTER*80 polec,model_name
        REAL*8 pole(3),x(3)
        INTEGER i,isave,model,verbos
  
  