anchor.diff:  C----
         IMPLICIT NONE
         CHARACTER*40 fltf
  -      INTEGER i,f,nflt,FMAX
  +      INTEGER i,f,nflt,FMAX,ln
         PARAMETER (FMAX=1500)
         REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
        1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX)
         open (unit=32,file=fltf,status='old')
  +      call linecount(fltf,ln)
         f = 0
     321 f = f + 1
         i = nflt + f
  
positive.code:           elseif (-135.0d0.lt.rak(i).and.rak(i).lt.-45.0d0) then
                typ = 3 ! Normal fault (-135<rake<-45)
            else
                typ = 1 ! Strike-slip faul (otherwise)
            endif
            call wellscoppersmith(wid,len,mag,typ)
            slip(i) = ((1.0d1**(1.5d0*(mag+10.7d0)))*1.0d-7)/
  
negative.code:       else
            mu = 4.3d10
        endif
        call linecount(magf,ln)
        open (unit=33,file=magf,status='old')
        f = 0
    331 f = f + 1
  