anchor.diff:  use io, only: stdout, stdin
   use trig, only: d2r
   use algebra, only: cross_product, normalize
  -use earth, only: get_pole, pole_geo2xyz, deg2km, semimajor_grs80, eccentricity_grs80, ellipsoid_geo2xyz
  +use earth, only: get_pole, pole_geo2xyz, pole_xyz2geo, deg2km
   
   use platemotion, only: input_file, &
                          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: character(len=16) :: plates(2)
  double precision :: pole(3)
  character(len=32) :: model_name
  character(len=32) :: convert_mode
  
  end module
  
  