anchor.diff:  C----
   C Read finite fault model
   C----
  -      call readffm(ifile,evlo,evla,evdp,str,dip,rak,dx,dy,slip,trup,
  -     1             hylo,hyla,seg,nseg,nflt)
  +      if (itype.eq.1) then
  +          call readffm(ifile,evlo,evla,evdp,str,dip,rak,dx,dy,slip,trup,
  +     1                 hylo,hyla,seg,nseg,nflt)
  +      elseif (itype.eq.2) then
  +          call readfsp(ifile,evlo,evla,evdp,str,dip,rak,dx,dy,slip,trup,
  +     1                 hylo,hyla,seg,nseg,nflt)
  +      endif
   
   C----
   C Write output files
  
positive.code:           write(*,*)
        endif
        write(*,*)
       1 'Usage: ff2gmt -f FFMFILE -slip SLIPFILE -time TIMEFILE',
       2                  ' -dep DEPFILE'
        write(*,*)
       1 '              -clip CLIPFILE|-clipseg CLIPFILE -epi EPIFILE ',
       2                 '-dc DCFILE'
  
negative.code:       PARAMETER (pi=4.0d0*datan(1.0d0),d2r=pi/180.0d0)
        CHARACTER*40 ifile,sfile,tfile,zfile,cfile,nfile,efile,dfile
        LOGICAL ex
        INTEGER i,nflt,FMAX,cf,itype
        PARAMETER (FMAX=1500)
        REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
       1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX),trup(FMAX),hylo,hyla
  