[cig-commits] r18679 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: . EXAMPLES/splay_faults/DATA/FAULT Post-processing src
ampuero at geodynamics.org
ampuero at geodynamics.org
Thu Jun 30 00:50:49 PDT 2011
Author: ampuero
Date: 2011-06-30 00:50:49 -0700 (Thu, 30 Jun 2011)
New Revision: 18679
Modified:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in
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/README_SPECFEM3D_FAULT
seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90
seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90
Log:
more cleanup
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in 2011-06-29 14:03:50 UTC (rev 18678)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/DATA/FAULT/Par_file_faults.in 2011-06-30 07:50:49 UTC (rev 18679)
@@ -6,27 +6,26 @@
1000 ! NTSNAP: time interation of snapshots
-1.0e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
-&BEGIN_FAULT dummy_idfault=1 /
+&BEGIN_FAULT /
&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 /
+&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.88e6, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0, 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, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.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, 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, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =140.0e3, zc= -15.0e3, lx=227.0e3, ly=40.0e3, lz=30000.0e0 /
+&BEGIN_FAULT /
&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 /
+&DIST2D shape='rectangle', val = -5.69e6, valh=0.0e0, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, 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, 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, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, 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, lx=45574.5247e3, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, 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
+&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3,r=0e0, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
Modified: 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 2011-06-29 14:03:50 UTC (rev 18678)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/FSEM3D_snapshot.m 2011-06-30 07:50:49 UTC (rev 18679)
@@ -4,6 +4,7 @@
%
% INPUTS isnap snapshot index, as in Snapshot*.bin file names
% fig [1] plot horizontal slip rate or not
+% *** this option has been temporarily disabled, the value is ignored ***
% dir ["."] directory containing the SPECFEM3D output data Snapshot*.bin
% fault [1] fault id
%
@@ -63,4 +64,4 @@
clear BinRead
-return
\ No newline at end of file
+return
Modified: 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 2011-06-29 14:03:50 UTC (rev 18678)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/Post-processing/Plot_snapshot_splay_faults.m 2011-06-30 07:50:49 UTC (rev 18679)
@@ -53,58 +53,10 @@
axis tight
set(hf,'XLim',XLIM)
+% This is the colormap of Wendt et al, but it does not work well yet with plotclr:
%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
-
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 14:03:50 UTC (rev 18678)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT 2011-06-30 07:50:49 UTC (rev 18679)
@@ -170,9 +170,11 @@
3. follow all the steps in section II above.
+4. Visualize the results with the matlab scripts in the directory EXAMPLES/*/post
+
IV. INPUT FILES
----------------
@@ -181,85 +183,142 @@
DATA/STATIONS List of stations outside the fault (see manual page 23).
-DATA/FAULT/Par_file_faults.in contains parameters of the fault:
+DATA/FAULT/Par_file_faults contains parameters of the fault. The first part of this file
+ has a strict format:
+
Line 1: Number of faults (NF)
- Lines 2 to NF+1: Two columns: domain tag on side #1 of the fault
- and domain tag on side #2 of the fault. These domain tags are assigned
- to blocks of elements on each side of the fault during mesh generation in CUBIT.
- #3 eta (Kelvin Voigt factor).
- Line NF+2: Number of time steps between updates of the time series outputs at selected
- fault points (see DATA/FAULT/FAULT_STATIONS.in), usually a large number (100s or 1000s).
+ Lines 2 to NF+1: Kelvin Voigt damping (in seconds) for each fault.
+ Line NF+2: Type of simulation (1=dynamic , 2 = kinematic)
+ Line NF+3: Number of time steps between updates of the time series outputs at selected
+ fault points (see DATA/FAULT/FAULT_STATIONS), usually a large number (100s or 1000s).
Note that the sampling rate of the time series is usually much higher.
- Line NF+3: Number of time steps between fault snapshot outputs (quantities at every fault
- point exported at selected times), usually a large number (100s or 1000s).
- Line NF+4: Slip velocity threshold below which frictional healing is set (friction coefficient
+ Line NF+4: Number of time steps between fault snapshot outputs (quantities at every fault
+ point exported at regular times), usually a large number (100s or 1000s).
+ Line NF+5: Slip velocity threshold below which frictional healing is set (friction coefficient
is reset to its static value). If this value is negative healing is disabled.
- Line NF+5: Slip velocity threshold to define the rupture front. Only used for outputs.
- Line NF+6: Fault parameters id.
- Line NF+7: Initial fault stress for fault #1:
- S1 = along-strike shear
- S2 = along-dip shear
- S3 = normal stress (negative in compresion)
- See figure 1 below.
- Line NF+8: Static and dynamic friction coefficients.
- Repeat NF+6 to NF+8 for each fault.
+ Line NF+6: Slip velocity threshold to define the rupture front. Only used for outputs.
-DATA/FAULT/FAULT_STATIONS.in Stations in the fault plane.
+ The rest of this file is made of namelist input blocks (see "namelist" in a Fortran 9x manual).
+ The input for each fault has the following sequence (arguments in [brackets] are optional):
+
+ &BEGIN_FAULT /
+ &INIT_STRESS S1, S2, S3 [,n1, n2, n3] /
+ followed by (n1+n2+n3) &DIST2D blocks
+ &SWF mus, mud, dc [, nmus, nmud, ndc] /
+ followed by (nmus+nmud+ndc) &DIST2D blocks
+
+ We adopt the usual coordinate system in seismology: (strike,dip,normal) where
+ the dip direction points downwards. This implies that the normal direction points out of
+ the footwall. See the figure below:
+
+ ^(normal)
+ .
+ .
+ ___________ . . > (strike)
+ | |
+ | Fault |
+ |_________|
+ .
+ .
+ v (dip)
+
+ The &INIT_STRESS input block sets the initial fault stresses relative to the foot-wall side of
+ the fault. Initial stresses are composed of a constant background value possibly overwritten
+ in prescribed regions by heterogeneous distributions (see &DIST2D blocks below):
+ S1 = constant background value of along-strike shear stress
+ (positive in the usual strike direction)
+ S2 = constant background value of along-dip shear
+ (positive is down-dip, normal faulting)
+ S3 = constant background value of normal stress (negative in compresion)
+ n1 = number of heterogeneous items for along-strike shear stress [default is 0]
+ n2 = number of heterogeneous items for along-dip shear stress [default is 0]
+ n3 = number of heterogeneous items for normal stress [default is 0]
+
+ The &SWF input block sets the slip-weakening friction parameters of the fault:
+ mus = constant background value of static friction coefficient
+ mud = constant background value of dynamic friction coefficient
+ dc = constant background value of critical slip-weakening distance
+ nmus = number of heterogeneous items for static friction coefficient [default is 0]
+ nmud = number of heterogeneous items for dynamic friction coefficient [default is 0]
+ ndc = number of heterogeneous items for critical slip-weakening distance [default is 0]
+
+ The &DIST2D input blocks modify (overwrite) the value of a fault parameter by a heterogeneous
+ spatial distribution:
+
+ &DIST2D shape='square', val, xc, yc, zc, l /
+ sets a constant value (val) within a cube with center (xc,yc,zc) and edge size l.
+
+ &DIST2D shape='rectangle', val, xc, yc, zc, lx, ly, lz /
+ sets a constant value (val) within a rectangular prism with center (xc,yc,zc)
+ and edge sizes (lx,ly,lz).
+
+ &DIST2D shape='rectangle-taper', val, valh, xc, yc, zc, lx, ly, lz /
+ sets a vertical linear gradient
+ within a rectangular prism with center (xc,yc,zc) and edge sizes (lx,ly,lz).
+ Values vary linearly as a function of vertical position z
+ between value val at z = zc-lz/2 and value valh at z = zc+lz/2 .
+
+ &DIST2D shape='circular', val, xc, yc, zc, r /
+ sets a constant value (val) within a sphere with center (xc,yc,zc) and radius r.
+
+
+
+DATA/FAULT/FAULT_STATIONS Stations in the fault plane.
Line 1: number of stations.
Line 2 to end: 5 columns: X, Y, Z (-depth), station name, fault-id
- The fault-id identifies the fault that contains the station.
- It is the index of appearance in the faults list after line 2 of Par_file_faults.in
+ The fault-id identifies the fault that contains the station.
+ It is the index of appearance in the faults list after line 2 of Par_file_faults
Several files are generated automatically by xgenerate_databases in directory DATABASES_MPI
and do not need to be modified by the user.
-
IV. OUTPUT FILES
-----------------
Several output files are saved in ~/SPECFEM3D/OUTPUT_FILES:
-1. Stations on the fault plane. Their locations and names are given in
- DATA/FAULT/FAULT_STATIONS.in. Their output files are named after the station.
+1. Seismograms for each station on the fault plane given in DATA/FAULT/FAULT_STATIONS.
+ One output file is generated for each station, named after the station. The files
+ are ascii and start with a header (22 lines long) followed by a data block with the
+ following format, one line per time sample:
-# Column #1 = Time (s)
-# Column #2 = horizontal right-lateral slip (m)
-# Column #3 = horizontal right-lateral slip rate (m/s)
-# Column #4 = horizontal right-lateral shear stress (MPa)
-# Column #5 = vertical up-dip slip (m)
-# Column #6 = vertical up-dip slip rate (m/s)
-# Column #7 = vertical up-dip shear stress (MPa)
-# Column #8 = normal stress (MPa)
+ # Column #1 = Time (s)
+ # Column #2 = horizontal right-lateral slip (m)
+ # Column #3 = horizontal right-lateral slip rate (m/s)
+ # Column #4 = horizontal right-lateral shear stress (MPa)
+ # Column #5 = vertical up-dip slip (m)
+ # Column #6 = vertical up-dip slip rate (m/s)
+ # Column #7 = vertical up-dip shear stress (MPa)
+ # Column #8 = normal stress (MPa)
-Fault coordinates system : (strike,dip,normal)
+ The stresses are relative to the footwall side of the fault (this convention controls
+ their sign, but not their amplitude). Slip is defined as displacement of the hanging
+ wall relative to the footwall.
- .(normal)
- .
- .
- -------------> (strike)
- | |
- | Fault |
- |---------
- .
- . (dip)
- Figure 1.
+2. Seismograms at stations in the bulk (out of the fault plane) given in DATA/STATIONS.
+ The name and format of these output files is described in page 51 of the manual.
-2. Stations in the bulk, outside the fault plane.
- Their locations and names are given in DATA/STATIONS.
- Output format: see manual page 51.
+3. Rupture time files are named Rupture_time_FAULT-id. One file is generated for each fault.
+ The files are ascii and start with a header (12 lines long) followed by a data block
+ with the following format, one line per fault node:
+ # Column #1 = X position (m)
+ # Column #2 = Y position (m)
+ # Column #3 = Z position (m)
+ # Column #3 = rupture time (s)
-3. Rupture time files are named Rupture_time_FAULT-id. Their format is 4 columns:
- X, Y, Z and rupture time.
-4. The face FAULT plane reference used in the computatinos takes fault_down as reference.
-
-
V. POST-PROCESSING AND VISUALIZATION
-------------------------------------
Some Matlab functions for post-processing and visualization are included in directory
Post-processing. The functions are internally documented (see their matlab help).
+
+ FSEM3D_snapshot reads a fault data snapshot
+
+The directories EXAMPLES/*/post contain additional Matlab scripts to generate figures
+specific to each example.
+
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90 2011-06-29 14:03:50 UTC (rev 18678)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90 2011-06-30 07:50:49 UTC (rev 18679)
@@ -89,12 +89,12 @@
! read fault input file
nb = 0
- open(unit=IIN,file='DATA/FAULT/Par_file_faults.in',status='old',action='read',iostat=ier)
+ open(unit=IIN,file='DATA/FAULT/Par_file_faults',status='old',action='read',iostat=ier)
if (ier==0) then
read(IIN,*) nb
else
write(6,*) 'No faults in the domain'
- write(6,*) 'Par_file_faults.in does not exist '
+ write(6,*) 'Par_file_faults does not exist '
close(IIN)
end if
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90 2011-06-29 14:03:50 UTC (rev 18678)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90 2011-06-30 07:50:49 UTC (rev 18679)
@@ -127,9 +127,9 @@
dummy_idfault = 0
- open(unit=IIN_PAR,file='DATA/FAULT/Par_file_faults.in',status='old',iostat=ier)
+ open(unit=IIN_PAR,file='DATA/FAULT/Par_file_faults',status='old',iostat=ier)
if( ier /= 0 ) then
- write(6,*) 'File Par_file_faults.in not found: assume no faults'
+ write(6,*) 'File Par_file_faults not found: assume no faults'
close(IIN_PAR)
return
endif
@@ -687,7 +687,7 @@
! requested coordinate
IIN = 251
- open(IIN,file='DATA/FAULT/FAULT_STATIONS.in',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
read(IIN,*) np
DataT%npoin =0
do i=1,np
@@ -701,7 +701,7 @@
allocate(DataT%iglob(DataT%npoin))
allocate(DataT%name(DataT%npoin))
- open(IIN,file='DATA/FAULT/FAULT_STATIONS.in',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
if( ier /= 0 ) stop 'error opening FAULT_STATIONS file'
read(IIN,*) np
k = 0
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90 2011-06-29 14:03:50 UTC (rev 18678)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90 2011-06-30 07:50:49 UTC (rev 18679)
@@ -107,9 +107,9 @@
dummy_idfault = 0
- open(unit=IIN_PAR,file='DATA/FAULT/Par_file_faults.in',status='old',iostat=ier)
+ open(unit=IIN_PAR,file='DATA/FAULT/Par_file_faults',status='old',iostat=ier)
if( ier /= 0 ) then
- write(6,*) 'Have not found Par_file_faults.in: assume no faults'
+ write(6,*) 'Have not found Par_file_faults: assume no faults'
return
endif
@@ -491,7 +491,7 @@
! requested coordinate
IIN = 251
- open(IIN,file='DATA/FAULT/FAULT_STATIONS.in',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
read(IIN,*) np
DataT%npoin =0
do i=1,np
@@ -505,7 +505,7 @@
allocate(DataT%iglob(DataT%npoin))
allocate(DataT%name(DataT%npoin))
- open(IIN,file='DATA/FAULT/FAULT_STATIONS.in',status='old',action='read',iostat=ier)
+ open(IIN,file='DATA/FAULT/FAULT_STATIONS',status='old',action='read',iostat=ier)
if( ier /= 0 ) stop 'error opening FAULT_STATIONS file'
read(IIN,*) np
k = 0
More information about the CIG-COMMITS
mailing list