anchor.diff:    516         continue
                 obj0(j) = obj
             endif
  -C         If solution is better than previous, save it
  +!         If solution is better than previous, save it
             if (obj.lt.objmin) then
                 do 519 j = 1,nflt
                     best(j,1) = soln(j,1)
                     best(j,2) = soln(j,2)
                     objmin = obj
     519         continue
  -C         When using Metropolis-Hastings, allow solution chance to get worse
  +!         When using Metropolis-Hastings, allow solution chance to get worse
             elseif (algrthm.eq.'metropolis') then
                 test = dexp((objmin-obj)/temp)
                 random = ran2(idum)
  
positive.code:       PROGRAM main
  C----
  C Invert displacement observations for fault slip, assuming slip in an
  C elastic half-space. Based on equations of Okada (1992).
  C
  C Inversion options:
  C     - Least squares with LAPACK tools
  C     - Simulated annealing using Gibbs sampler (THIS DOES NOT WORK YET!)
  C----
        IMPLICIT none
  C Command line inputs
        CHARACTER*80 fltf,obsf,smoof,haff,ofile,annf,lkfile
        REAL*8 smooth,damp,damp0,fact,stscon
        INTEGER geo,p,invert,vrb,ostyle,misfit,getdmp
  C Inversion array variables
        INTEGER nobs,OBSMAX,nflt,FLTMAX
        PARAMETER (OBSMAX=1000,FLTMAX=1000)
        REAL*8 obs(OBSMAX,6),flt(FLTMAX,7),gf(OBSMAX,FLTMAX,6),
       1       sgf(FLTMAX,FLTMAX,4)
        REAL*8 soln(FLTMAX,2)
        REAL*8 pre(OBSMAX,3),tot
        INTEGER isThisFaultLocked(FLTMAX),nsts
  C Other variables
        REAL*8 haf(3),mxdisp,disp
        INTEGER i,j,k
        COMMON /VERBOSE/ vrb
  
  C----
  C Parse command line and test inputs
  C----
        call gcmdln(fltf,obsf,smoof,smooth,damp,ofile,p,geo,haff,invert,
       1            annf,ostyle,misfit,getdmp,fact,stscon,lkfile)
        if (vrb.ge.1) then
            write(0,'("FLTINV VERBOSE OUTPUT ON")')
        endif
  
negative.code:           nsoln = 1
            return
        endif
  ! Check if annealing control file exists
        inquire(file=annf,exist=ex)
        if (.not.ex) then
            write(*,*) '!! Warning: no annealing file '//trim(annf)
  