anchor.diff:   1002 format(F12.6," =",F12.6," + (",F12.6,"*",F12.6,") + (",F12.6,"*",
        1       F12.6,")")
   
  -C Record solutions in solutions.log
  +! Record solutions in solutions.log
    1003 format('Iteration: ',I4,4X,'Temperature: ',F9.6,4X,'Objective=',
        1       F12.6)
    1004 format("  flt      slip      rake")
  
positive.code:  1002 format('    Misfit with damping = ',1PE10.3,' is ',1PE14.6)
        END
  
  C----------------------------------------------------------------------C
  
        SUBROUTINE anneal(soln,gf,obs,nobs,OBSMAX,nflt,FLTMAX,damp,smooth,
       1                  smoof,annf)
  C----
  C Perform a random search for the optimal solution, utilizing a
  C simulated annealing algorithm.
  C----
        IMPLICIT none
        REAL*8 pi,d2r
        PARAMETER (pi=4.0d0*atan(1.0d0),d2r=pi/1.8d2)
  C Inversion variables
        INTEGER OBSMAX,FLTMAX
        REAL*8 obs(OBSMAX,6),gf(OBSMAX,FLTMAX,6)
        REAL*8 soln(FLTMAX,2),best(FLTMAX,2)
  
negative.code:       RETURN
        END
  
  !----------------------------------------------------------------------C
  
        SUBROUTINE usage(str)
        IMPLICIT none
  