anchor.diff:            if (kshr.eq.1.or.knor.eq.1.or.kcou.eq.1) then
                 call coulomb(coul,norml,shear,stress,trgstr,trgdip,trgrak,
        1                     frict)
  -              if (kshr.eq.1) then
  +              if (kshr.eq.1.and.noz.eq.0) then
                     write(25,9992) stlo,stla,stdp*1d-3,shear
  +              elseif (kshr.eq.1.and.noz.eq.1) then
  +                  write(25,*) stlo,stla,shear
                 endif
  -              if (knor.eq.1) then
  +              if (knor.eq.1.and.noz.eq.0) then
                     write(24,9992) stlo,stla,stdp*1d-3,norml
  +              elseif (knor.eq.1.and.noz.eq.1) then
  +                  write(24,*) stlo,stla,norml
                 endif
  -              if (kcou.eq.1) then
  +              if (kcou.eq.1.and.noz.eq.0) then
                     write(26,9992) stlo,stla,stdp*1d-3,coul
  +              elseif (kcou.eq.1.and.noz.eq.1) then
  +                  write(26,*) stlo,stla,coul
                 endif
             endif
   C         Update progress indicator
  
positive.code: C Parse command line
  C----
        call gcmdln(ffmf,fltf,magf,haff,staf,trgf,dspf,stnf,stsf,norf,
       1            shrf,coulf,gmtf,flttyp,auto,incr,adist,xy,prog,long)
  
  C----
  C Check for input fault source files, and whether output is specified
  
negative.code:       CHARACTER*40 staf,trgf,dspf,stnf,stsf,norf,shrf,coulf
        INTEGER i,nsta,ntrg,progout
        REAL*8 stlo,stla,stdp,trgstr,trgdip,trgrak,frict,vp,vs,dens
        INTEGER nflt,FMAX,flttyp,xy,prog,long,noz
        PARAMETER (FMAX=1500)
        REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
       1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX)
  