anchor.diff:                 endif
                  fixed_rad = .true.
                  iarg = iarg + 2
  +            case('-v')
  +               value = .true.
  +               iarg = iarg + 1
               case('-h')
                  call usage
               case default
  
positive.code:          stop 1
        endif
        call geog2cart(lon, lat, r, x, y, z, degrees=.true.)
        write(*,*) x, y, z
     enddo
  
  contains
  
negative.code:    iostat = 0
     do while (iostat == 0)
        if (fixed_rad) then
           if (value) then
              read(*,*,iostat=iostat) lon, lat, val
           else
              read(*,*,iostat=iostat) lon, lat
           endif
        else
           if (value) then
              read(*,*,iostat=iostat) lon, lat, r, val
           else
              read(*,*,iostat=iostat) lon, lat, r
           endif
        endif
        if (iostat < 0) exit
        if (iostat > 0) then
  