anchor.diff:            write(*,*)
         endif
         write(*,*)
  -     1 '-geo                Use geographic coordinates in 3-D problem'
  +     1 '-geo                Use geographic coordinates'
         if (d.eq.1) then
             write(*,*)
  -     1 '                       X=E, Y=N, Z pos down, uz pos up'
  -          write(*,*)
  -      endif
  -      write(*,*)
  -     1 '-haf HAFSPC         Half-space parameters'
  +     1 '                       X->E, Y->N'
  +          write(*,*)
  +      endif
  +      write(*,*)
  +     1 '-haf HAFSPC         Half-space parameter file'
         if (d.eq.1) then
             write(*,*)
        1 '                       vp(m/s) vs(m/s) dens(kg/m^3)'
             write(*,*)
  -     1 '                       Lame lamda mu'
  +     1 '                       Lame lamda(Pa) mu(Pa)'
             write(*,*)
         endif
         write(*,*)
        1 '-damp CONST         Minimize length of model, with weight CONST'
  +      if (d.eq.1) then
  +          write(*,*)
  +      endif
         write(*,*)
        1 '-smooth CONST FILE  Smooth model, with weight CONST, linking ',
  -     2                         'file FILE'
  +     2                         'FILE'
         if (d.eq.1) then
             write(*,*)
  -     1 '                      Linking file format: flt #_neighbors ',
  -     2                                  'neighbor_1 neighbor_2...'
  +     1 '                        Linking file format: flt nghbrs ',
  +     2                                  'nghbr_1-N'
             write(*,*)
        1 '                        1    2     2  4     '
             write(*,*)
  
positive.code:       SUBROUTINE damping(A,b,nobs,nflt,damp)
  C----
  C Minimize length of solution to Ax = b with weighting factor damp.
  C m=nobs n=nflt d=damp
  C         [ A11 A12 ... A1n ]
  C     A = [ A21 A22 ... A2n ]
  C         [  :           :  ]
  C         [ Am1 Am2 ... Amn ]
  C
  C         [ d*d  0  ...  0  ]
  C  Append [  0  d*d ...  0  ]
  C         [                 ]
  C         [  0   0  ... d*d ]
  C----
        IMPLICIT none
        REAL*8 damp
        INTEGER OBSMAX,FLTMAX
        PARAMETER (OBSMAX=2500,FLTMAX=2000)
        REAL*8 A(OBSMAX,FLTMAX),b(OBSMAX,1)
        INTEGER i,j,nobs,nflt
        if (damp.le.0.0d0) return
        do 101 i = 1,nflt
            b(nobs+i,1) = 0.0d0
            do 102 j = 1,nflt
                if (i.eq.j) then
                    A(nobs+i,j) = damp*damp*1.0d0
                else
                    A(nobs+i,j) = 0.0d0
                endif
    102     continue
    101 continue
        nobs = nobs + nflt
        RETURN
        END
  
  C----------------------------------------------------------------------C
  
        SUBROUTINE smoothing(A,b,nobs,nflt,smooth,smoof,dimen)
  C----
  C Minimize roughness of solution to Ax = b with weighting factor smooth.
  C m=nobs n=nflt d=damp
  C         [ A11 A12 ... A1n ]
  C     A = [ A21 A22 ... A2n ]
  C         [  :           :  ]
  C         [ Am1 Am2 ... Amn ]
  C
  C  (e.g.) [ d*d -d*d    0     0  ...   0   0  ]
  C  Append [-d*d 2*d*d -d*d    0  ...   0   0  ]
  C         [  0  -d*d  2*d*d -d*d ...   0   0  ]
  C         [  0    0     0     0  ... -d*d d*d ]
  C----
        IMPLICIT none
        CHARACTER*80 smoof
        LOGICAL ex
        REAL*8 smooth
        INTEGER OBSMAX,FLTMAX
        PARAMETER (OBSMAX=2500,FLTMAX=2000)
        REAL*8 A(OBSMAX,FLTMAX),b(OBSMAX,1)
        INTEGER nsmoo,el,nnei,nei(FLTMAX)
        INTEGER i,j,nobs,nflt,dimen
        if (smooth.le.0.0d0) return
  C Check for input smoothing file
        if (smoof.eq.'none') then
  
negative.code:      1 '-o OFILE            Output slip values'
        if (d.eq.1) then
            write(*,*)
       1 '                       str_slip(m) dip_slip(m)'
            write(*,*)
        endif
        write(*,*)
  