anchor.diff:  C      do 101 i=1,nflt
   C          print *,evlo(i),evla(i),evdp(i),slip(i)
   C  101 continue
  -      close(34) 
  +      close(34)
         RETURN
         END
   
  
positive.code:               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.and.kdsp.eq.1) then
                    write(21,9990) stlo,stla,stdp*1d-3,uE,uN,uZ
                elseif (long.eq.1.and.kdsp.eq.1) then
                    write(21,9993) stlo,stla,stdp*1d-3,uE,uN,uZ
                elseif (kdsp.eq.2) then
                    az = atan2(uE,uN)*r2d
                    hdsp = dsqrt(uE*uE+uN*uN)
  
negative.code:      1                                         'SUBFAULTS'
            write(0,*)
        endif
  
  C----
  C Check for auto flag to define stations
  C----
  