[cig-commits] commit: adding examples
Mercurial
hg at geodynamics.org
Tue Sep 20 12:12:58 PDT 2011
changeset: 3:6f8a209b1e84
user: Sylvain Barbot <sylbar.vainbot at gmail.com>
date: Thu Jan 06 17:46:27 2011 -0800
files: examples/mogi.sh examples/relax.sh examples/run1.sh examples/run1_proj.sh examples/run2.sh examples/run3.sh examples/run4.sh examples/run5.sh
description:
adding examples
diff -r 5eba818786ac -r 6f8a209b1e84 examples/mogi.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/mogi.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,53 @@
+#!/bin/sh
+
+# viscoelastic relaxation following a dilatation source
+# viscous flow governed by a power-law viscosity with
+# stress exponent of n = 3.
+
+time ../relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size & smoothing (dx1,dx2,dx3,beta) & nyquist
+0.05 0.05 0.05 0.2 1
+# origin position and rotation
+0 0 0
+# observation depth (displacement and stress)
+0 0.5
+# output directory
+./mogi
+# elastic parameters (lambda,mu) and gravity effect
+1 1 0
+# integration time (t1)
+10 0.257
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+0
+# number of powerlaw viscous interfaces
+1
+# no depth gammadot0 power cohesion
+ 1 2 20 3 0
+# number of power-law ductile zones
+0
+# number of friction faults
+0
+# number of interseismic loading stuff
+0
+0
+# number of coseismic events
+1
+# number of shear dislocations
+0
+# index slip x1 x2 x3 length width strike dip rake
+# 1 0 -1.0 0 0 1 0.8 0 90 0
+# number of tensile cracks
+0
+# number of dilatation point sources
+1
+# no dilatation xs ys zs
+ 1 1 0 0 1
+EOF
diff -r 5eba818786ac -r 6f8a209b1e84 examples/relax.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/relax.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,53 @@
+#!/bin/sh
+
+time ./relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size & smoothing (dx1,dx2,dx3,beta)
+0.05 0.05 0.05 0.2
+# origin position
+0 0
+# observation depth
+0
+# output directory
+./output
+# elastic parameters (lambda,mu)
+1 1
+# integration time (t1)
+20 0.5
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+2
+1 1.0 0 0.0
+2 9.0 0 0.0
+# number of powerlaw viscous interfaces
+2
+1 1.0 1e1 3.0 0.0
+2 9.0 1e1 3.0 0.0
+# number of friction faults
+0
+# number of interseismic loading stuff
+0
+0
+# number of coseismic events
+2
+# number of shear dislocations
+1
+# index slip x1 x2 x3 length width strike dip rake
+ 1 1 -1.0 0 0 1 0.8 0 90 0
+# number of tensile cracks
+0
+# time of second event
+10
+# number of shear dislocations
+1
+# index slip x1 x2 x3 length width strike dip rake
+ 1 0.02 0.0 0 0 1 0.8 0 90 0
+# number of tensile cracks
+0
+EOF
diff -r 5eba818786ac -r 6f8a209b1e84 examples/run1.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/run1.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,78 @@
+#!/bin/sh
+
+# 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.
+
+time ../relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
+0.05 0.05 0.05 0.2 2
+# origin position & rotation
+0 0 0
+# observation depth (displacement and stress) (stress in only exported in GRD)
+0 0.5
+# output directory
+./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)
+20 -1
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+2
+# no depth gammadot0 cohesion
+ 1 3.0 1.0 0.0
+ 2 9.0 1.0 0.0
+# number of linear ductile zones
+0
+# number of powerlaw viscous interfaces
+0
+# number of friction faults
+0
+# number of interseismic loading strike-slip and opening
+0
+0
+# number of coseismic events
+2
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 1 -1 0 0 1 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+# time of second event
+1
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 1 0 0 0 1 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+EOF
diff -r 5eba818786ac -r 6f8a209b1e84 examples/run1_proj.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/run1_proj.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,86 @@
+#!/bin/sh
+
+# linear viscoelastic relaxation following a strike-slip fault
+# followed by a second strike-slip event and a relaxation of
+# both stress perturbations.
+
+# output is exported in geographic coordintes in .xyz format
+
+# 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.
+
+time ../relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
+0.05 0.05 0.05 0.2 2
+# origin position & rotation
+0 0 0
+# geographic origin (longitude, latitude, UTM zone, unit)
+# displacements and stress are converted to lon/lat geographic coordinates
+# unit corresponds to a scaling from dx1,dx2,dx3 to real unit
+# use unit = 1 if dimensions are described in units of m
+# use unit = 1e3 if dimensions are described in units of km
+-115 37.5 11 1e3
+# observation depth (displacement and stress) (stress in only exported in GRD)
+0 0.5
+# output directory
+./output1_proj
+# 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)
+20 -1
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+2
+# no depth gammadot0 cohesion
+ 1 3.0 1.0 0.0
+ 2 9.0 1.0 0.0
+# number of linear ductile zones
+0
+# number of powerlaw viscous interfaces
+0
+# number of friction faults
+0
+# number of interseismic loading strike-slip and opening
+0
+0
+# number of coseismic events
+2
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 1 -1 0 0 1 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+# time of second event
+1
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 1 0 0 0 1 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+EOF
diff -r 5eba818786ac -r 6f8a209b1e84 examples/run2.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/run2.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,61 @@
+#!/bin/sh
+
+# nonlinear viscoelastic relaxation (power exponent n=3)
+# following a strike-slip fault
+
+# output every computational step (Dt = -2)
+
+# to visualize coseismic deformation (requires GRD output, or manual conversion to GRD format):
+# map.sh -b -3/3/-3/3 output2/000
+
+# to visualize postseismic deformation (requires GRD output):
+# map.sh -b -5/5/-5/5 output2/0{02,04,06,08}-relax
+
+# type map.sh for a description of command-line options.
+
+time ../relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
+0.05 0.05 0.05 0.2 2
+# origin position & rotation
+0 0 0
+# observation depth (displacements and stress) (stress is only exported in GRD)
+0 0.5
+# output directory
+./output2
+# 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)
+50 -2
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+0
+# number of powerlaw viscous interfaces
+2
+# no depth gammadot0 power cohesion
+ 1 3.0 5e3 3.0 0.0
+ 2 8.0 5e3 3.0 0.0
+# number of power-law ductile zones
+0
+# number of friction faults
+0
+# number of interseismic loading strike-slip and opening
+0
+0
+# number of coseismic events
+1
+# number of shear dislocations
+1
+# index slip x1 x2 x3 length width strike dip rake
+ 1 1 1 0 0 2 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+EOF
diff -r 5eba818786ac -r 6f8a209b1e84 examples/run3.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/run3.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,78 @@
+#!/bin/sh
+
+# linear viscoelastic relaxation following a strike-slip fault
+# followed by a second strike-slip event and a relaxation of
+# both stress perturbations.
+
+# output is exported in geographic coordintes in .xyz format
+
+# 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.
+
+time ../relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
+0.05 0.05 0.05 0.2 1
+# origin position & rotation
+0 0 0
+# geographic origin (longitude, latitude, UTM zone, unit)
+# displacements and stress are converted to lon/lat geographic coordinates
+# unit corresponds to a scaling from dx1,dx2,dx3 to real unit
+# use unit = 1 if dimensions are described in units of m
+# use unit = 1e3 if dimensions are described in units of km
+#-115 37.5 11 1e3
+# observation depth (displacement and stress) (stress in only exported in GRD)
+0 0
+# output directory
+./output3
+# 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)
+2 -1
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+0
+# number of powerlaw viscous interfaces
+0
+# number of friction faults
+1
+# no depth gamma0 (a-b)sig friction cohesion
+ 1 1 1e3 1e3 0.6 0
+# number of creeping faults
+1
+# no x1 x2 x3 length width strike dip
+ 1 -1 0 1 2 1 0 90
+# number of interseismic loading strike-slip
+0
+# number of interseismic loading opening cracks
+0
+# number of coseismic events
+1
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 -1 -1 0 0 2 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+EOF
diff -r 5eba818786ac -r 6f8a209b1e84 examples/run4.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/run4.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,76 @@
+#!/bin/sh
+
+# linear viscoelastic relaxation following a strike-slip fault
+# followed by a second strike-slip event and a relaxation of
+# both stress perturbations.
+
+# output is exported in geographic coordintes in .xyz format
+
+# 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.
+
+time ../relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
+0.05 0.05 0.05 0.2 1
+# origin position & rotation
+0 0 0
+# geographic origin (longitude, latitude, UTM zone, unit)
+-120 34 11 1e3
+# observation depth (displacement and stress) (stress in only exported in GRD)
+0 0
+# output directory
+./output4
+# 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)
+2 -1
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+3
+# no depth gammadot0 cohesion
+ 1 2 1 0
+ 2 3 1 0
+ 3 3 0 0
+# number of ductile zones
+1
+# no gammadot0 x1 x2 x3 length width thickness strike dip
+ 1 1 -1 0 1.1 2 0.9 1 0 90
+# number of powerlaw viscous interfaces
+0
+# number of friction faults
+0
+# number of interseismic loading strike-slip
+0
+# number of interseismic loading opening cracks
+0
+# number of coseismic events
+1
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 -1 -1 0 0 2 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+EOF
diff -r 5eba818786ac -r 6f8a209b1e84 examples/run5.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/examples/run5.sh Thu Jan 06 17:46:27 2011 -0800
@@ -0,0 +1,80 @@
+#!/bin/sh
+
+# 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.
+
+time ../relax <<EOF
+# grid size (sx1,sx2,sx3)
+256 256 256
+# sampling size, smoothing & nyquist (dx1,dx2,dx3,beta,nq)
+0.05 0.05 0.05 0.2 2
+# origin position & rotation
+0 0 0
+# geographic origin (longitude, latitude, UTM zone, unit)
+131.9 43.0 52 1e3
+# observation depth (displacement and stress) (stress in only exported in GRD)
+0 0.5
+# output directory
+./output5
+# 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)
+20 -1
+# number of observation planes
+0
+# number of observation points
+0
+# number of prestress interfaces
+0
+# number of linear viscous interfaces
+2
+# no depth gammadot0 cohesion
+ 1 3.0 1.0 0.0
+ 2 9.0 1.0 0.0
+# number of linear ductile zones
+0
+# number of powerlaw viscous interfaces
+0
+# number of friction faults
+0
+# number of interseismic loading strike-slip and opening
+0
+0
+# number of coseismic events
+2
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 1 -1 0 0 1 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+# time of second event
+1
+# number of shear dislocations
+1
+# no slip x1 x2 x3 length width strike dip rake
+ 1 1 0 0 0 1 1 0 90 0
+# number of tensile cracks
+0
+# number of dilatation sources
+0
+EOF
More information about the CIG-COMMITS
mailing list