anchor.diff:   9991 format(3F12.4,X,6E16.6)
    9992 format(3F12.4,X,1E16.6)
    9993 format(3F12.4,X,3F16.8)
  + 9994 format(3F12.4,X,1P3E16.6)
         RETURN
         END
   
  
positive.code:       SUBROUTINE chksta(sta,staf,nsta,xy)
        IMPLICIT none
        REAL*8 sta(3)
        CHARACTER*40 staf 
        INTEGER i,nsta,chk,xy,ios
        chk = 0
        open(unit=18,file=staf,status='old')
  
negative.code: C      do 101 i=1,nflt
  C          print *,evlo(i),evla(i),evdp(i),slip(i)
  C  101 continue
        close(34)
        RETURN
        END
  
  