anchor.diff:            dy(i) = 1.0d3*dy(i)
             goto 321
     322 continue
  +      if (f-1.ne.ln) then
  +          call usage('!! Error: NFLT in -flt file not equal to NLINE')
  +      endif
         nflt = i - 1
         close(32)
         RETURN
  
positive.code:       else
            mu = 4.3d10
        endif
  C      print *,'SHEAR MOD:',mu
        open (unit=33,file=magf,status='old')
        f = 0
    331 f = f + 1
  
negative.code: C----
        IMPLICIT NONE
        CHARACTER*40 magf
        INTEGER i,f,nflt,FMAX,ln
        PARAMETER (FMAX=1500)
        REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
       1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX),mag
  