anchor.diff:       2                                    'compression.'
             write(*,*)
        1 '    This differs from many geological applications, where ',
  -     2          'greater normal'
  +     2          'positive normal'
             write(*,*)
        1 '    stress is compressive.'
             write(*,*)
  
positive.code:           elseif (tag(1:4).eq.'-gmt') then
                kgmt = 1
                i = i + 1
                call getarg(i,gmtfile)
            elseif (tag(1:2).eq.'-h'.or.tag(1:5).eq.'-help') then
                call usage(' ')
            elseif (tag(1:2).eq.'-d'.or.tag(1:8).eq.'-details') then
  
negative.code:       COMMON /ICHECK/ kffm,kflt,kmag,ksta,khaf,ktrg
        COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
  C     Open files specified for output
        if (kdsp.eq.1) open(unit=21,file=dspf,status='unknown')
        if (kstn.eq.1) open(unit=22,file=stnf,status='unknown')
        if (ksts.eq.1) open(unit=23,file=stsf,status='unknown')
        if (knor.eq.1) open(unit=24,file=norf,status='unknown')
        if (kshr.eq.1) open(unit=25,file=shrf,status='unknown')
        if (kcou.eq.1) open(unit=26,file=coulf,status='unknown')
  C     Initialize progress indicator
        if (prog.eq.1) then
            i = 0
            progout = 0
            call progbar(i,nsta,progout)
        endif
  C     If NTRG=1, read target fault parameters here
        if (ktrg.eq.2) read(trgf,*) trgstr,trgdip,trgrak,frict
        if (ktrg.eq.1) open(unit=13,file=trgf,status='old')
        if (ntrg.eq.1.and.ktrg.eq.1) read(13,*) trgstr,trgdip,trgrak,frict
  C     Compute deformation for each receiver location
        open(unit=12,file=staf,status='old')
        do 201 i = 1,nsta
            read(12,*) stlo,stla,stdp
            stdp = stdp*1.0d3
  C         If NTRG=1, read target fault parameters here
            if (ntrg.eq.nsta.and.ktrg.eq.1.and.nsta.ne.1) then
                read(13,*) trgstr,trgdip,trgrak,frict
            endif
  C         Compute displacements
            if (kdsp.eq.1) then
                call calcdisp(uN,uE,uZ,stlo,stla,stdp,nflt,evlo,evla,evdp,
       1                      str,dip,rak,dx,dy,slip,vp,vs,dens,flttyp,xy)
                if (long.eq.0) then
                    write(21,9990) stlo,stla,uE,uN,uZ
                else
                    write(21,9993) stlo,stla,uE,uN,uZ
                endif
            endif
  C         Compute (3x3) strain matrix
            if (kstn.ge.1) then
  