anchor.diff:            write (13,9999) stlo,stla,norml
             write (14,*) stlo,stla,enet(1,1),enet(2,2),enet(3,3),
        1                           enet(1,2),enet(1,3),enet(2,3)
  +
  +          prog = prog + 1
  +          call progbar(prog,prog100,progtag)
             goto 101
     107 continue
   
  
positive.code: C----
  C Calculate net strain, stress, coulomb stress change at each station
  C----
    101 read (22,*,end=107) stlo,stla,stdp
            stdp = 1.0d3*stdp
            do 103 i = 1,3
                do 102 j = 1,3
  
negative.code:       open (unit=23,file=haffile,status='old')
        open (unit=24,file=trgfile,status='old')
        read (23,*) vp,vs,dens
  
  C----
  C Compare number of lines in stations.txt and targetparams.txt
  C and use line count of stations.txt to initialize progress indicator
  C----
        stact = 0
    112 read (22,*,end=113)
            stact = stact + 1
            goto 112
    113 continue
        trgct = 0
    114 read (24,*,end=115)
             trgct = trgct + 1
             goto 114
    115 continue
        rewind 22
        rewind 24
  
        prog100 = stact
        if (trgct.eq.1) then
            if (verbose) write (*,8997)
        elseif (trgct.eq.stact) then
            if (verbose) write (*,8998)
            trgmatchsta = .true.
        else
            write (*,8999)
            call usage()
        endif
   8997 format ('NL in targetparams.txt = 1')
   8998 format ('NL in targetparams.txt = NL in stations.txt')
   8999 format ('Error: NL in targetparams.txt must be 1 or',
       1        ' equal to NL in stations.txt')
  
  C----
  C Output files
  C----
  