anchor.diff:            dx(i) = len*1.0d3
             dy(i) = wid*1.0d3
             evdp(i) = 1.0d3*evdp(i)
  -          goto 201
  -  202 continue
  +          goto 331
  +  332 continue
         nflt = i - 1
         close(33)
         RETURN
  
positive.code:      1       rak(FMAX),slip(FMAX),dx(FMAX),dy(FMAX),area
        INTEGER flttyp
        REAL*8 vp,vs,dens
  C----
  C Compute NEZ displacements at (stlo,stla,stdp)
  C----
        uNnet = 0.0d0
        uEnet = 0.0d0
        uZnet = 0.0d0
        do 102 f = 1,nflt
            call ddistaz(dist,az,evlo(f),evla(f),stlo,stla)
            dist = dist*6.371d6
            x = dist*( dcos(az-d2r*str(f)))
  
negative.code:           read (dyc,*) dyr
            dxr = 1.0d3*dxr
            dyr = 1.0d3*dyr
            read (31,*) du,du,du,du,du,du,du,du,du,du,hylo,du,hyla
            do 311 i=1,7
                read (31,*) du
    311     continue
            do 312 i = 1,nx*ny
                read (31,*) evla(ct+i),evlo(ct+i),evdp(ct+i),slip(ct+i),
       1                    rak(ct+i),str(ct+i),dip(ct+i)
                slip(ct+i) = 1.0d-2*slip(ct+i)
                dx(ct+i) = dxr
                dy(ct+i) = dyr
                evdp(ct+i) = 1.0d3*evdp(ct+i)
    312     continue
            ct = ct + nx*ny
    313 continue
        nflt = nflt + ct
        close(31)
        RETURN
  