anchor.diff:  C     Open files specified for output
         if (kdsp.ne.0) open(unit=21,file=dspf,status='unknown')
         if (kstn.eq.1) open(unit=22,file=stnf,status='unknown')
  -      if (ksts.eq.1.or.ksts.eq.3) 
  +      if (ksts.eq.1.or.ksts.eq.3)
        1    open(unit=23,file=stsf,status='unknown')
         if (knor.eq.1) open(unit=24,file=norf,status='unknown')
         if (kshr.ge.1) open(unit=25,file=shrf,status='unknown')
  
positive.code:       endif
        write(*,*)
       1 '-coul COULFILE          Coulomb stress (',
       2                            'requires target fault kinematics)' 
        if (str.eq.'long') then
            write(*,*)
            write(*,*)
  
negative.code:      1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX)
        REAL*8 uE,uN,uZ,stdp,vp,vs,dens,adist
        INTEGER xy,dir,flttyp
        if (evdp(1).gt.1.0d5) then
            xmin = x0 - 3.0d0/dcos(evla(1)*0.01745d0)
            xmax = x0 + 3.0d0/dcos(evla(1)*0.01745d0)
            ymin = y0 - 3.0d0
            ymax = y0 + 3.0d0
            return
        endif
        if (xy.eq.0) then
            dl(1) = 1.0d0 ! Amounts to move test location (deg)
            dl(2) = 0.1d0
  