anchor.diff:                call calcstn(strain,stlo,stla,stdp,nflt,
        1                     evlo,evla,evdp,str,dip,rak,dx,dy,slip,
        2                     vp,vs,dens,flttyp,xy)
  -              if (kstn.eq.1.and.noz.eq.0) then
  +              if (kstn.eq.1) then
                     write(22,9991) stlo,stla,stdp*1d-3,
        1                           strain(1,1),strain(2,2),strain(3,3),
        2                           strain(1,2),strain(1,3),strain(2,3)
  -              elseif(kstn.eq.1.and.noz.eq.1) then
  -                  write(22,*) stlo,stla,
  -     1                           strain(1,1),strain(2,2),strain(3,3),
  -     2                           strain(1,2),strain(1,3),strain(2,3)
                 endif
             endif
   C         Compute (3x3) stress matrix
             if (ksts.ge.1) then
                 call stn2sts(stress,strain,vp,vs,dens)
  -              if (ksts.eq.1.and.noz.eq.0) then
  +              if (ksts.eq.1) then
                     write(23,9991) stlo,stla,stdp*1d-3,
        1                           stress(1,1),stress(2,2),stress(3,3),
        2                           stress(1,2),stress(1,3),stress(2,3)
  -              elseif (ksts.eq.1.and.noz.eq.1) then
  -                  write(23,*) stlo,stla,
  -     1                           stress(1,1),stress(2,2),stress(3,3),
  -     2                           stress(1,2),stress(1,3),stress(2,3)
  -              elseif (ksts.eq.3.and.noz.eq.0) then
  +              elseif (ksts.eq.3) then
                     write(23,9991) stlo,stla,stdp*1d-3,estrs(stress)
  -              elseif (ksts.eq.3.and.noz.eq.1) then
  -                  write(23,*) stlo,stla,estrs(stress)
                 endif
             endif
   C         Resolve stresses onto target fault planes
             if (kshr.ge.1.or.knor.eq.1.or.kcou.eq.1) then
                 call coulomb(coul,norml,shear,stress,trgstr,trgdip,trgrak,
        1                     frict)
  -              if (kshr.ge.1.and.noz.eq.0) then
  +              if (kshr.ge.1) then
                     write(25,9992) stlo,stla,stdp*1d-3,shear
  -              elseif (kshr.ge.1.and.noz.eq.1) then
  -                  write(25,*) stlo,stla,shear
                 endif
  -              if (knor.eq.1.and.noz.eq.0) then
  +              if (knor.eq.1) 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.and.noz.eq.0) then
  +              if (kcou.eq.1) 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:       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
  
negative.code:      2                               '[-haf HAFFILE] [-xy]'
        write(*,*)
       1 '               -disp DSPFILE [-az] -strain STNFILE -stress ',
       2                             'STSFILE -estress ESTSFILE'
        write(*,*)
       1 '               -normal NORFILE -shear SHRFILE -coul COULFILE'
        write(*,*)
       1 '               [-auto h|d|s DIST INCR] [-long] [-prog] ',
       2                   '[-gmt GMTFILE]'
        write(*,*)
       1 '               [-h|-help] [-d|-details]'
  