anchor.diff:  C----
   
         IMPLICIT none
  -      REAL olon,olat,ilon,ilat,plon,plat,angle
  -      REAL ix,iy,iz,px,py,pz,ox,oy,oz
  +      REAL*8 olon,olat,ilon,ilat,plon,plat,angle
  +      REAL*8 ix,iy,iz,px,py,pz,ox,oy,oz
         CHARACTER*30 ifile,ofile
   
         call gcmdln(ifile,ofile)
  
positive.code:   102     continue
        endif
  
   9999 format(2F8.2)
  
        END
  
  
negative.code: 
        SUBROUTINE rotate(olon,olat,ilon,ilat,plon,plat,angle)
        IMPLICIT NONE
        REAL*8 olon,olat,ilon,ilat,plon,plat,angle
        REAL*8 ix,iy,iz,px,py,pz,ox,oy,oz
        REAL*8 pi,d2r,r2d
        PARAMETER (pi=4.0d0*atan(1.0d0),d2r=pi/1.8d2,r2d=1.8d2/pi)
  
        call lonlat2xyz(ix,iy,iz,ilon,ilat)
        call lonlat2xyz(px,py,pz,plon,plat)
  