anchor.diff:  write(stderr,*) '                       MORVEL56 (Argus et al., 2011)'
   write(stderr,*) '                       NUVEL1A  (DeMets et al., 1994)'
   write(stderr,*) '-list MODEL        List plates in MODEL (use ALL or leave blank to see all models)'
  +write(stderr,*) '-geo2xyz           Convert Euler pole in geographic coordinates to Cartesian'
  +write(stderr,*) '-xyz2geo           Convert Euler pole in Cartesian coordinates to geographic (assumes LON/LAT/VEL is X/Y/Z)'
   write(stderr,*) '-f IFILE           Input file (default: stdin)'
   write(stderr,*) '-o OFILE           Output file (default: stdout)'
   write(stderr,*)
  
positive.code: 
  ! Parse command line
  call gcmdln()
  
  
  ! Calculate pole of rotation between plate pair
  
negative.code:                        output_file, &
                         plates, &
                         pole, &
                         model_name, &
                         convert_mode
  
  implicit none
  
  