anchor.diff:                ostyle = 11
             elseif (tag(1:4).eq.'-sci') then
                 ostyle = 10
  +          elseif (tag(1:4).eq.'-sts') then
  +              i = i + 1
  +              call getarg(i,tag)
  +              read(tag,'(BN,F10.0)') stscon
  +              i = i + 1
  +              if (i.gt.narg) goto 102
  +              call getarg(i,tag)
  +              indx = index(tag,'-')
  +              if (indx.eq.0) then
  +                  read(tag,*) lkfile
  +              else
  +                  i = i - 1
  +              endif
             elseif (tag(1:7).eq.'-misfit') then
                 misfit = 1
             elseif (tag(1:9).eq.'-getdamp1') then
  
positive.code: C Run inversion
  C----
        if (invert.eq.0) then
            call lstsqr(soln,gf,obs,nobs,OBSMAX,nflt,FLTMAX,damp,smooth,
       1                smoof)
        elseif (invert.eq.1) then
            call anneal(soln,gf,obs,nobs,OBSMAX,nflt,FLTMAX,damp,smooth,
       1                smoof,annf)
        else
            call usage('!! Error: inversion option does not exist')
        endif
  
  C----
  C Misfit
  C----
        if (misfit.eq.1) then
            tot = 0.0d0
            do 905 i = 1,nobs
                pre(i,1) = 0.0d0
                pre(i,2) = 0.0d0
                pre(i,3) = 0.0d0
                do 906 j = 1,nflt
                    pre(i,1) = pre(i,1) + gf(i,j,1)*soln(j,1) 
       1                                             + gf(i,j,4)*soln(j,2)
  
negative.code:      1                                                      trim(smoof))
        endif
        call linect(nsmoo,smoof)
        if (nsmoo.gt.m) then
            write(*,*) '!! Error: more faults to smooth than m'
            call usage('!! Check smoothing linking file')
        endif
  C Initialize smoothing section of A and b arrays as zeros
        do 422 i = 1,nsmoo
            b(n+i,1) = 0.0d0
            do 421 j = 1,m
                A(n+i,j) = 0.0d0
    421     continue
    422 continue
  C Read smoothing element linking file and load into A
  