[cig-commits] commit: cancel reading of geographic location altogether when --no-proj-output is used

[Mercurial]

changeset:   13:6355ba48e6e7
user:        Sylvain Barbot <sylbar.vainbot at gmail.com>
date:        Thu May 12 14:28:43 2011 -0700
files:       examples/run1.sh
description:
cancel reading of geographic location altogether when --no-proj-output is used


diff -r 54d6b6185e08 -r 6355ba48e6e7 examples/run1.sh
--- a/examples/run1.sh	Tue Apr 26 17:20:05 2011 -0700
+++ b/examples/run1.sh	Thu May 12 14:28:43 2011 -0700
@@ -1,41 +1,60 @@
 #!/bin/sh
 
+#
+# comments:
+#
 # linear viscoelastic relaxation following a strike-slip fault
 # followed by a second strike-slip event and a relaxation of
 # both stress perturbations.
-
+#
 # viscous substrate starts at 3 times the seismogenic depth.
 # output at each computation step (Dt = -1)
-
+#
 # to visualize coseismic deformation (requires GRD output, or manual conversion to GRD format):
 # map.sh -b -3/3/-3/3 output1/000
-
+#
 # to convert .xyz output in .grd output in post processing, type (requires GMT installed on your machine)
 # xyz2grd.sh output1/000
-
+#
 # to visualize a time series of postseismic deformation (requires GRD output): 
 # map.sh -b -5/5/-5/5 -p -0.002/0.002/0.0001 -v 0.005 output1/0{01,02,03,04,05,06,07,08,09,10}-relax
 # or simply
 # map.sh -b -5/5/-5/5 -p -0.002/0.002/0.0001 -v 0.005 output1/0??-relax
-
+#
 # type map.sh for a description of command-line options.
+#
+# to visualize in 3-D with Paraview (www.paraview.org),
+# load the following files:
+# 
+#   output1/cgrid.vtp           for the computational grid spatial extent
+#   output1/rfaults-001.vtp     for the slip distribution of the first event
+#   output1/rfaults-002.vtp     for the slip distribution of the second event
+#   output1/disp-000.vtr        for the coseismic displacement field
+#   output1/linearlayer-001.vtp for the depth of brittle-ductile transition
+#   output1/vel-001.vtr         for the instantaneous velocity field
+#
+# note that we save the entire output to output1/in.param. this is convenient
+# to document how the simulation was computed.
+#
+# type relax --help to display options.
+#
 
 time ../relax <<EOF | tee output1/in.param
 # grid size (sx1,sx2,sx3)
 256 256 256
-# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
+# sampling size (in unit of length), smoothing (0-0.5) & nyquist (dx1,dx2,dx3,beta,nq)
 0.05 0.05 0.05 0.2 2
 # origin position & rotation
 0 0 0
-# geographic origin, UTM zone and length unit
+# geographic origin (longitude and latitude), UTM zone and real length unit (usually m or km)
 120 22 51 1e3
-# observation depth (displacement and stress) (stress in only exported in GRD)
+# observation depth for displacement and for stress (stress in only exported in GRD)
 0 0.5
-# output directory
+# output directory (all output written here)
 ./output1
 # elastic parameters and gamma = (1-nu) rho g / mu = 8.33e-7 /m = 8.33e-4 /km
 1 1 8.33e-4
-# integration time (t1), step and scaling
+# integration time (in unit of time), step (negative for automatic) and scaling of computed value
 20 -1 1
 # number of observation planes
 0
@@ -43,9 +62,9 @@ 0
 0
 # number of prestress interfaces
 0
-# number of linear viscous interfaces
+# number of linear viscous interfaces (where viscosity changes)
 2
-# no depth gammadot0 cohesion
+# no depth gammadot0 cohesion (gammadot0 is shear modulus divided by viscosity)
    1   3.0       1.0      0.0
    2   9.0       1.0      0.0
 # number of linear ductile zones
@@ -57,9 +76,9 @@ 0
 # number of interseismic loading strike-slip and opening
 0
 0
-# number of coseismic events
+# number of coseismic events (when slip distribution is prescribed)
 2
-# number of shear dislocations
+# number of shear dislocations (strike-slip and dip-slip faulting)
 1
 # no slip x1 x2 x3 length width strike dip rake
    1    1 -1  0  0      1     1      0  90    0