anchor.diff:       2                               '[-haf HAFFILE] [-xy]'
         write(*,*)
        1 '               -disp DSPFILE [-az] -strain STNFILE -stress ',
  -     2                             'STSFILE'
  +     2                             'STSFILE -estress ESTSFILE'
         write(*,*)
        1 '               -normal NORFILE -shear SHRFILE -coul COULFILE'
         write(*,*)
  -     1 '               [-auto AB VALUE INCR] [-long] [-prog] ',
  +     1 '               [-auto h|d|s DIST INCR] [-long] [-prog] ',
        2                   '[-gmt GMTFILE]'
         write(*,*)
        1 '               [-h|-help] [-d|-details]'
  
positive.code:               typ = 1 ! Strike-slip fault (otherwise)
            endif
            call wellscoppersmith(wid,len,mag,typ)
            slip(i) = ((1.0d1**(1.5d0*(mag+10.7d0)))*1.0d-7)/
       1                                            (mu*wid*len*1.0d6)
            dx(i) = len*1.0d3
  
negative.code: C----
        call calcdefm(staf,trgf,nsta,ntrg,nflt,evlo,evla,evdp,str,dip,rak,
       1              dx,dy,slip,flttyp,vp,vs,dens,dspf,stnf,stsf,norf,
       2              shrf,coulf,xy,prog,long)
  
  C----
  C Write fault input to GMT-compatible file
  