anchor.diff:        if (kmag.eq.1) then
             call readmag(magf,evlo,evla,evdp,str,dip,rak,dx,dy,slip,
        1                 nflt)
  +      endif
  +      if (nflt.lt.1) then
  +          write(*,*) '!! Error: no input faults read'
  +          call usage('!! Check input files for correct format')
         endif
   C     If FFM not used, use unweighted mean coordinates as hylo and hyla
         if (kffm.ne.1) then
  
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: C----
        IMPLICIT NONE
        CHARACTER*40 fltf
        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
  