anchor.diff:        PARAMETER (FMAX=1500)
         REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
        1       rak(FMAX),slip(FMAX),area(FMAX),dx(FMAX),dy(FMAX)
  -      INTEGER flttyp(FMAX)
  +      INTEGER flttyp
         REAL*8 vp,vs,dens
   
         uNnet = 0.0d0
  
positive.code: C----
  C Parse the command line, check for FFM file, and read FFM
  C----
        call gcmdln(fltfile,stafile,haffile,dspfile,auto)
        inquire(file=fltfile,exist=ex)
        if (.not.ex) call usage('!! Error: no file named '//fltfile)
        call readfaults(fltfile,nflt,evlo,evla,evdp,str,dip,rak,slip,dy,
       1                dx,area,flttyp)
  
  C----
  C If auto=0: use existing input files
  
negative.code:           else
                open (unit=23,file=haffile,status='old')
                read (23,*) vp,vs,dens
                close(23)
            endif
            inquire(file=stafile,exist=ex)
            if (.not.ex) call usage('!! Error: no file named '//stafile)
            open(unit=22,file=stafile,status='old')
        elseif (auto.eq.1) then
            call sethafdef(vp,vs,dens)
            call autogrid(stafile,nflt,evlo,evla,evdp,str,
       1                  dip,rak,dx,dy,area,slip,vp,vs,dens,flttyp)
            open (unit=22,file=stafile,status='old')
        endif
  
  