anchor.diff:  ! Calculate stress Green's function for each fault-fault pair
   do i = 1,fault%nrows
   
  -    ! Station coordinates
  +    ! Station (target fault) coordinates
       if (gf_model.eq.'okada_rect') then
           sta(1) = fault%array(i,1)
           sta(2) = fault%array(i,2)
  
positive.code:         if (coord_type.eq.'geographic') then
              call tri_geo2cart(pt1,pt2,pt3,fault%array(i,1:3),fault%array(i,4:6), &
                                fault%array(i,7:9),'m')
          endif
          call tri_geometry(nvec,tvec1,tvec2,pt1,pt2,pt3)
      endif
  
      do j = 1,fault%nrows
  
negative.code: ! enddo
  
  ! Calculate chi-squared
  if (verbosity.ge.3) then
      write(stdout,*) 'anneal_psc_objective: starting misfit_chi2'
  endif
  call misfit_chi2(obs,pre,cov_matrix,nobs,anneal_psc_objective)
  anneal_psc_objective = -0.5d0*anneal_psc_objective
  if (verbosity.ge.3) then
      write(stdout,*) 'anneal_psc_objective: finished misfit_chi2'
  endif
  
  if (verbosity.ge.3) then
      write(stdout,*) 'anneal_psc_objective: objective=',anneal_psc_objective
  