anchor.diff:  if (auto_mode.eq.'h'.or.auto_mode.eq.'v') then
       call update_auto_stations()
       call calc_deformation()
  +    if (auto_mode.eq.'v') then
  +        call project_deformation()
  +    endif
   endif
   
   end
  
positive.code: if (verbosity.ge.1) then
      write(stdout,*) 'update_auto_stations: finished'
      write(stdout,*)
  endif
  
  return
  
negative.code: integer :: i, j, ierr
  double precision :: disp(2002,6), disp_mag, xmin, xmax, ymin, ymax
  double precision :: d, lon, lat
  character(len=512) :: proj_sta_file
  
  
  if (verbosity.ge.1) then
  