anchor.diff:  
   C----------------------------------------------------------------------C
   
  -      SUBROUTINE smoothing(A,b,nobs,nflt,smooth,smoof)
  +      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
  
positive.code: C----
  C Add smoothing
  C----
        call smoothing(A,b,nobs,nflt,smooth,smoof)
        
  C----
  C Solve generalized least squares problem
  C----
        call lsq(A,b,x,OBSMAX,PARMAX,nobs,nflt)
        open(unit=23,file='fltinv.out',status='unknown')
        do 103 i = 1,nflt
            write(23,*) x(i,1)
    103 continue
  
negative.code: 
  C----------------------------------------------------------------------C
  
        SUBROUTINE gcmdln(fltf,obsf,smooth,smoof,damp,dimen,ofile,geo)
        IMPLICIT NONE
        CHARACTER*80 fltf,obsf,tag,smoof,ofile
        INTEGER i,narg
        REAL*8 smooth,damp
        INTEGER dimen,geo
        fltf = 'none'
        obsf = 'none'
        smoof = 'none'
        ofile = 'none'
        dimen = 0
        smooth = -1.0d0
        damp = -1.0d0
        geo = 0
        narg = iargc()
        if (narg.eq.0) call usage('!! Error: no command line arguments')
        i = 0
  