[cig-commits] r18676 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: . EXAMPLES/splay_faults/DATA/FAULT Post-processing
percygalvez at geodynamics.org
percygalvez at geodynamics.org
Wed Jun 29 05:51:20 PDT 2011
Author: percygalvez
Date: 2011-06-29 05:51:20 -0700 (Wed, 29 Jun 2011)
New Revision: 18676
Added:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in
seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/
seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/FSEM3D_snapshot.m
seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m
seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/plotclr.m
Modified:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
Log:
New fault database done
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.barrier.in 2011-06-29 12:51:20 UTC (rev 18676)
@@ -0,0 +1,32 @@
+2
+0.00134 ! #tag #eta(damping) FAULT A
+0.00134 ! #tag #eta(damping) FAULT BC
+1 ! 1 = dyn 2=kin
+100 ! NTOUT : Number of time steps
+1000 ! NTSNAP: time interation of snapshots
+-1.0e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
+&BEGIN_FAULT dummy_idfault=1 /
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
+&DIST2D shape='rectangle', val = -2.16e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.88e6, valh =0.0e0, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, r=0e0, l=0.0e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
+&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
+&DIST2D shape='rectangle', val = 0.6e0, valh = 0.0e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, valh = 0.0e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+&BEGIN_FAULT dummy_idfault=2 /
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
+&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle', val = -3.54e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, valh=0.0e0, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, r=0e0, l=0e0, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
+&DIST2D shape='rectangle', val = 0.6e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&KINPAR kindt=10.0e-3
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.homogeneous.in 2011-06-29 12:51:20 UTC (rev 18676)
@@ -0,0 +1,32 @@
+2
+0.00134 ! #tag #eta(damping) FAULT A
+0.00134 ! #tag #eta(damping) FAULT BC
+1 ! 1 = dyn 2=kin
+100 ! NTOUT : Number of time steps
+1000 ! NTSNAP: time interation of snapshots
+-1.0e0 ! V_HEALING (-1 : Healing off)
+1e-3 ! V_RUPT
+&BEGIN_FAULT dummy_idfault=1 /
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
+&DIST2D shape='rectangle', val = -2.16e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.88e6, valh =0.0e0, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, r=0e0, l=0.0e0, lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
+&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, valh =0.0e0, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
+&DIST2D shape='rectangle', val = 0.6e0, valh = 0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, valh = 0.0e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, r=0e0, l=0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&BEGIN_FAULT dummy_idfault=2 /
+&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
+&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle', val = -3.54e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, valh=0.0e0, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, r=0e0, l=0e0, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3,r=0e0, l=0e0, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 55819.1065e0, yc =120.0e3, zc= -4.5e3, r=0e0, l=0e0, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, valh=0.0e0, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, r=0e0, l=0e0, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
+&DIST2D shape='rectangle', val = 0.6e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, valh=0.0e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, l=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&KINPAR kindt=10.0e-3
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/FSEM3D_snapshot.m
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/FSEM3D_snapshot.m (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/FSEM3D_snapshot.m 2011-06-29 12:51:20 UTC (rev 18676)
@@ -0,0 +1,66 @@
+%FSEM3D_SNAPSHOT reads and plots fault fields at a given time
+%
+% d = FSEM3D_snapshot(isnap, [fig, dir, fault])
+%
+% INPUTS isnap snapshot index, as in Snapshot*.bin file names
+% fig [1] plot horizontal slip rate or not
+% dir ["."] directory containing the SPECFEM3D output data Snapshot*.bin
+% fault [1] fault id
+%
+% OUTPUTS d structure containing fault fields:
+% X,Y,Z fault node coordinates (km)
+% Dx,Dz slip (m)
+% Vx,Vz slip rate (m/s)
+% Tx,Ty,Tz stress change (MPa)
+% S slip "state" variable in friction law (m)
+% Sg strength relative to initial stress (MPa)
+% Trup rupture time (s)
+% Tpz process zone time, when slip=Dc (s)
+%
+% Jean-Paul Ampuero ampuero at erdw.ethz.ch modified by
+% Percy Galvez percy.galvez at sed.ethz.ch 19/01/2011.
+%
+% WARNING : This script only do snapshots of single precision files so far.
+
+function d = FSEM3D_snapshot(isnap,fig,DATA_DIR,fault)
+
+NDAT = 14;
+VM = 3; % saturate the color scale of the slip rate figure at VM m/s
+
+if nargin<2, fig=1; end
+if nargin<3, DATA_DIR = '.'; end
+if nargin<4, fault = 1; end
+
+BinFile = sprintf('%s/Snapshot%u_F%u.bin',DATA_DIR,isnap,fault);
+
+if ~exist(BinFile,'file'), error(sprintf('File %s does not exist',BinFile)), end
+fid=fopen(BinFile);
+BinRead = fread(fid,[1,inf],'single')' ;
+fclose(fid);
+
+BinRead = reshape( BinRead(:),length(BinRead)/(NDAT),NDAT);
+BinRead = BinRead(2:end-1,:);
+
+% Reorder fault nodes (lexicographic z,x)
+%[LOC,IND] = sortrows( BinRead(:,[1 3]) );
+%BinRead = BinRead(IND,:);
+
+d.X = BinRead(:,1)/1e3; % in km
+d.Y = BinRead(:,2)/1e3; % in km
+d.Z = BinRead(:,3)/1e3; % in km
+d.Dx = BinRead(:,4);
+d.Dz = BinRead(:,5);
+d.Vx = BinRead(:,6);
+d.Vz = BinRead(:,7);
+d.Tx = BinRead(:,8); % in MPa
+d.Ty = BinRead(:,9);
+d.Tz = BinRead(:,10); % in MPa
+d.S = BinRead(:,11);
+d.Sg = BinRead(:,12); % in MPa
+d.Trup = BinRead(:,13);
+d.Tpz = BinRead(:,14);
+
+clear BinRead
+
+
+return
\ No newline at end of file
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m 2011-06-29 12:51:20 UTC (rev 18676)
@@ -0,0 +1,110 @@
+clear all
+close all
+ITSNAP = 33000;
+XLIM = [0 120]; % range along-strike
+SLIM = [-8 8]; % slip range
+
+dA = FSEM3D_snapshot(ITSNAP,1,'.',1);
+dBC = FSEM3D_snapshot(ITSNAP,1,'.',2);
+
+% along-dip slip is negative for thrust faulting
+dA.Dz = -dA.Dz;
+dBC.Dz = -dBC.Dz;
+
+% shift to Wendt's coordinate system
+dA.Y = 180-dA.Y;
+dA.X = 227-dA.X;
+dBC.Y = 180-dBC.Y;
+dBC.X = 90-dBC.X;
+
+figure(1)
+subplot(221)
+plot(dA.Y,dA.Dx,'.')
+set(gca,'XLim',XLIM)
+ylabel('Dx (m)')
+subplot(223)
+plot(dA.Y,dA.Dz,'.')
+set(gca,'XLim',XLIM)
+ylabel('Dz (m)')
+xlabel('Y along-strike (km)')
+subplot(222)
+plot(dBC.Y,dBC.Dx,'.')
+set(gca,'XLim',XLIM)
+subplot(224)
+plot(dBC.Y,dBC.Dz,'.')
+set(gca,'XLim',XLIM)
+
+figure(2)
+subplot(121)
+[hc,hf]=plotclr(dA.Y,dA.X,dA.Dz,'s',SLIM);
+set(hf,'YDir','reverse')
+xlabel('Along strike (km)')
+ylabel('Down dip (km)')
+axis equal
+axis tight
+set(hf,'XLim',XLIM)
+
+subplot(122)
+[hc,hf]=plotclr(dBC.Y,dBC.X,dBC.Dz,'s',SLIM);
+set(hf,'YDir','reverse')
+xlabel('Along strike (km)')
+ylabel('Down dip (km)')
+axis equal
+axis tight
+set(hf,'XLim',XLIM)
+
+%colormap('jet');
+%map1 = colormap;
+%l1 = length(map1);
+%map2 = map1(l1/2+1:end,:);
+%colormap(map2);
+
+
+%----------------------------------------------------------
+%----------------------------------------------------------
+%----------------------------------------------------------
+% old version, works for regular (cartesian) fault grids:
+
+% % Reorder fault nodes
+% [LOC,IND] = sortrows( BinRead(:,[3 1]) );
+% NX = length( unique(LOC(:,2)) );
+% NZ = length( unique(LOC(:,1)) );
+%
+% BinRead = reshape( BinRead(IND,:), NX,NZ,NDAT );
+%
+% d.X = BinRead(:,:,1)/1e3; % in km
+% d.Y = BinRead(:,:,2)/1e3; % in km
+% d.Z = BinRead(:,:,3)/1e3; % in km
+% d.Dx = BinRead(:,:,4);
+% d.Dz = BinRead(:,:,5);
+% d.Vx = BinRead(:,:,6);
+% d.Vz = BinRead(:,:,7);
+% d.Tx = BinRead(:,:,8); % in MPa
+% d.Ty = BinRead(:,:,9);
+% d.Tz = BinRead(:,:,10); % in MPa
+% d.S = BinRead(:,:,11);
+% d.Sg = BinRead(:,:,12); % in MPa
+% d.Trup = BinRead(:,:,13);
+% d.Tpz = BinRead(:,:,14);
+%
+% clear BinRead
+%
+% % plot
+% if fig
+% % clf
+% figure(1)
+% surf(d.X,d.Z,max(0,-d.Vx))
+% axis equal
+% view(2)
+% shading flat
+% zlim=get(gca,'ZLim');
+% set(gcf,'renderer','Zbuffer')
+% zlim=max(abs(zlim));
+% set(gca,'ZLim',[0 100])
+% caxis([0 VM])
+% title('Slip rate (m/s)')
+% xlabel('Distance along strike (km)')
+% ylabel('Depth (km)')
+% colorbar
+% end
+
Added: seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/plotclr.m
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/plotclr.m (rev 0)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/plotclr.m 2011-06-29 12:51:20 UTC (rev 18676)
@@ -0,0 +1,63 @@
+% function h = plotclr(x,y,v,marker,vlim)
+% plots the values of v colour coded
+% at the positions specified by x and y.
+% A colourbar is added on the right side of the figure.
+%
+% The colourbar strectches from the minimum value of v to its
+% maximum.
+%
+% 'marker' is optional to define the marker being used. The
+% default is a point. To use a different marker (such as circles, ...) send
+% its symbol to the function (which must be enclosed in '; see example).
+%
+% 'vlim' is optional, to define the limits of the colourbar.
+% v values outside vlim are not plotted
+%
+% modified by Stephanie Contardo, CSIRO, 2009
+% from 'plotc' by Uli Theune, University of Alberta, 2004
+%
+
+function [h,hf] = plotclr(x,y,v,marker,vlim)
+
+if nargin <4
+ marker='.';
+end
+
+map=colormap;
+if nargin >4
+ miv = vlim(1) ;
+ mav = vlim(2) ;
+else
+ miv=min(v);
+ mav=max(v);
+end
+clrstep = (mav-miv)/size(map,1) ;
+% Plot the points
+hold on
+for nc=1:size(map,1)
+ iv = find(v>miv+(nc-1)*clrstep & v<=miv+nc*clrstep) ;
+ plot(x(iv),y(iv),marker,'color',map(nc,:),'markerfacecolor',map(nc,:))
+end
+hf=gca;
+hold off
+
+% Re-format the colorbar
+h=colorbar;
+
+%set(h,'ylim',[1 length(map)]);
+yal=linspace(1,length(map),10);
+set(h,'ytick',yal);
+% Create the yticklabels
+ytl=linspace(miv,mav,10);
+s=char(10,4);
+for i=1:10
+ if min(abs(ytl)) >= 0.001
+ B=sprintf('%-4.3f',ytl(i));
+ else
+ B=sprintf('%-3.1E',ytl(i));
+ end
+ s(i,1:length(B))=B;
+end
+set(h,'yticklabel',s);
+grid on
+view(2)
\ No newline at end of file
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT 2011-06-29 10:38:57 UTC (rev 18675)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT 2011-06-29 12:51:20 UTC (rev 18676)
@@ -9,6 +9,7 @@
Details about the original package are described in its manual (manual_SPECFEM3D.pdf).
The current README file describes how to install and run our modified version the code.
+It is not a replacement for the original manual, which you should read first.
This document assumes that the root directory containing the code is ~/SPECFEM3D.
@@ -86,7 +87,7 @@
>> cd ~/SPECFEM3D/CUBIT
>> cubit
In CUBIT's menu "Tools", select "Play Journal File" and select a script file,
- for instance EXAMPLES/tpv5/fault_tpv5.py
+ for instance EXAMPLES/splay_faults/splay_faults.py
This creates several mesh files in directory CUBIT/MESH/:
absorbing_surface_file_* (5 files)
@@ -96,7 +97,7 @@
nodes_coords_file
nummaterial_velocity_file
- The CUBIT graphics window should show the mesh (e.g. EXAMPLES/tpv5/planar_fault.png)
+ The CUBIT graphics window should show the mesh (e.g. EXAMPLES/splay_faults/splay_faults.jpg)
2. Partition the mesh with the domain decomposition software SCOTCH.
@@ -112,13 +113,19 @@
and specify the path accordingly when executing xdecompose_mesh_SCOTCH.
-3. Edit the file DATA/Par_file:
+3. Copy the input files to the work directory:
+ >> cd ~/SPECFEM3D
+ >> cp -r EXAMPLES/splay_faults/DATA ./
+
+
+4. Edit the file DATA/Par_file:
+
LOCAL_PATH = should be the path to DATABASES_MPI
NPROC = number of processors. The same as the number of partitions in SCOTCH (step 2).
-4. Generate databases:
+5. Generate databases:
>> cd ~/SPECFEM3D
>> mpirun -np nproc ./xgenerate_databases
@@ -129,25 +136,34 @@
in directory DATABASES_MPI.
-5. Run the solver:
+6. If this is the first time you run a simulation with this mesh, edit DATA/Par_file to make sure that the variable DT (the time step)
+ is smaller than the "Maximum suggested time step" in OUTPUT_FILES/output_mesher.txt. ALso edit variables NSTEP and NTSTEP_* accordingly.
+
+7. Run the solver:
+
>> mpirun -np nproc ./xspecfem3D
or
>>./run/run.xspecfem3d
-
III. EXAMPLES
--------------
The package includes a few examples, the SCEC benchmark problems:
+ TPV5, a planar vertical strike-slip fault
- + TPV14 and TPV15, vertical strike-slip fault system with a fault branch
+ (+ TPV14 and TPV15, vertical strike-slip fault system with a fault branch)
+ + Splay fault models from Wendt et al. (2009)
-To run the examples first replace the contents of directory DATA/* by one of the EXAMPLES/tpv*/DATA,
-for instance:
+To run the examples:
+1. Read the documents in the directory EXAMPLES/*/description. They contain a description of the example
+ and additional instructions to run it.
+
+2. replace the contents of directory DATA/* by one of the EXAMPLES/tpv*/DATA,
+ for instance:
+
>> cp -r EXAMPLES/tpv5/DATA ./
Then follow all the steps in section II above.
More information about the CIG-COMMITS
mailing list