figure
hold on	
        
        %plot projections of fault planes
        [nf, ten] = size(dis_geom);
        for i=1:nf
         displot(dis_geom(i,:))
        end
        title('Geometry in Local Cartesian Coordiantes'), drawnow
       plotfault
       
       
       epxy=[18.9473  -41.5610];
         plot(epxy(1),epxy(2),'*')
         text(epxy(1),epxy(2),'epicenter')
         
         
         scale=2; 
         quiver(xy(:,1),xy(:,2),scale*d(1:3:length(d)),scale*d(2:3:length(d)),0)
         
         quiver(xy(:,1),xy(:,2),scale*dhat(1:3:length(d)),scale*dhat(2:3:length(d)),0,'r')
         
         axis tight
         axis equal
 