anchor.diff:  ! 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
  
positive.code: ! Calculate stress Green's function for each fault-fault pair
  do i = 1,fault%nrows
  
      ! Station coordinates
      if (gf_model.eq.'okada_rect') then
          sta(1) = fault%array(i,1)
          sta(2) = fault%array(i,2)
  
negative.code:     elseif (trim(tag).eq.'-los') then
          i = i + 1
          call get_command_argument(i,los%file)
      elseif (trim(tag).eq.'-prests'.or.tag.eq.'-prestress') then
          i = i + 1
          call get_command_argument(i,prestress%file)
      elseif (trim(tag).eq.'-cov') then
  