anchor.diff:        REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
        1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX),hylo,hyla
         REAL*8 vp,vs,dens
  -      INTEGER nsta,ntrg
  -      INTEGER kffm,kflt,kmag,ksta,khaf,ktrg        ! kvar definitions
  +      INTEGER nsta,ntrg,typ(FMAX)
  +      INTEGER kffm,kflt,kmag,ksta,khaf,ktrg,kvol        ! kvar definitions
         INTEGER kdsp,kstn,ksts,knor,kshr,kcou,kgmt   ! in gcmdln comments
  -      COMMON /ICHECK/ kffm,kflt,kmag,ksta,khaf,ktrg
  +      COMMON /ICHECK/ kffm,kflt,kmag,ksta,khaf,ktrg,kvol
         COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
   
   C----
  
positive.code:       endif
        if (kmag.eq.1) then
            call readmag(magf,evlo,evla,evdp,str,dip,rak,dx,dy,slip,
       1                 nflt,emprel)
        endif
        if (nflt.lt.1) then
            write(*,*) '!! Error: no input faults read'
  
negative.code: C----
        call gcmdln(ffmf,fltf,magf,haff,staf,trgf,dspf,stnf,stsf,norf,
       1            shrf,coulf,gmtf,flttyp,auto,incr,adist,xy,prog,long,
       2            emprel,volf)
  
  C----
  C Check for input fault source files, and whether output is specified
  C----
        call chksrc(ffmf,fltf,magf,volf)
        call chkout()
        if (kvol.eq.1) goto 101
  C----
  C Read input fault files
  C----
        call readsrc(ffmf,fltf,magf,nflt,evlo,evla,evdp,str,dip,rak,dx,dy,
       1             slip,hylo,hyla,emprel,haff)
  
  C----
  C Check for auto flag to define stations
  