anchor.diff:    102     continue
         endif
   
  - 9999 format(2F8.2)
  + 9999 format(2F14.6)
   
         END
   
  
positive.code: 
        SUBROUTINE rotate(olon,olat,ilon,ilat,plon,plat,angle)
        IMPLICIT NONE
        REAL olon,olat,ilon,ilat,plon,plat,angle
        REAL ix,iy,iz,px,py,pz,ox,oy,oz
        REAL pi,d2r,r2d
        PARAMETER (pi=3.14159265,d2r=pi/180.,r2d=180./pi)
  
        call lonlat2xyz(ix,iy,iz,ilon,ilat)
        call lonlat2xyz(px,py,pz,plon,plat)
  
negative.code:       SUBROUTINE lonlat2xyz(x,y,z,lon,lat)
  
        IMPLICIT none
        REAL*8 pi,d2r
        PARAMETER (pi=4.0d0*atan(1.0d0),d2r=pi/1.8d2)
        REAL*8 lon,lat,x,y,z
  
        x = cos(lon*d2r)*cos(lat*d2r)
        y = sin(lon*d2r)*cos(lat*d2r)
  