anchor.diff:                         output_file, &
                          plates, &
                          pole, &
  -                       model_name
  +                       model_name, &
  +                       convert_mode
   
   implicit none
   
  
positive.code: plates = ''
  pole = 0.0d0
  model_name = 'MORVEL56'
  
  ! Number of arguments
  narg = command_argument_count()
  
negative.code: 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,*)
  