anchor.diff:  ! Local variables
   integer :: i, j, ios, ierr
   character(len=32) :: fmt_str
  +character(len=512) :: input_line
   
   
   ierr = 0
  
positive.code: ! Read the file, in free format
  open(unit=21,file=val%file,status='old')
  do i = 1,val%nrows
      read(21,*,iostat=ios) (val%array(i,j),j=1,val%ncols)
      if (ios.ne.0) then
          write(stderr,*) 'read_fltinv_data_file: read error on file ',trim(val%file),' at line ',i
          ierr = 1
          return
      endif
  
negative.code:                 nn = ndisp_dof
                  mm = ndisp_dof
              else
                  if (index(disp_components,nchar).le.0) then
                      write(stderr,*) 'read_inputs: first component (',nchar,') is not used in inversion'
                      write(stderr,*) 'It is not being loaded into the covariance matrix.'
                      cycle
                  endif
                  if (index(disp_components,mchar).le.0) then
                      write(stderr,*) 'read_inputs: second component (',mchar,') is not used in inversion'
                      write(stderr,*) 'It is not being loaded into the covariance matrix.'
                      cycle
                  endif
                  nn = (index(disp_components,nchar)-1)*displacement%nrows
                  mm = (index(disp_components,mchar)-1)*displacement%nrows
              endif
  