[cig-commits] commit: fixes geometry of rfaults.xy for shifted origin coordinates, and color of cover picture.

[Mercurial]

changeset:   67:db5d0c6c7ab3
tag:         tip
user:        Sylvain Barbot <sbarbot at caltech.edu>
date:        Thu Jan 05 00:27:28 2012 -0800
files:       examples/run2.input export.f90 input.f90 latex/graphics/cover.eps latex/relax.tex
description:
fixes geometry of rfaults.xy for shifted origin coordinates, and color of cover picture.


diff -r 561a8af1a923 -r db5d0c6c7ab3 examples/run2.input
--- a/examples/run2.input	Fri Dec 30 17:19:40 2011 -0800
+++ b/examples/run2.input	Thu Jan 05 00:27:28 2012 -0800
@@ -1,3 +1,9 @@
+# nonlinear viscoelastic relaxation (power exponent n=3)
+# following slip on a strike-slip fault
+#
+# output every two computational steps (Dt = -2). 
+# The time step is computed automatically and the value unchanged (scale = 1)
+#
 # grid size (sx1,sx2,sx3)
 256 256 256
 # sampling size defines the grid spacing in units of distance
@@ -18,7 +24,7 @@ 1 1 8.33e-4
 1 1 8.33e-4
 # integration time refers to the length in units of time of the simulation.
 # integration time (t1)
-100 -1 1
+100 -1 0.5
 # observation planes are arbitrary planes where the inelastic strain rate
 # is sampled for output.
 # number of observation planes
@@ -66,4 +72,3 @@ 0
 # surface loads can be used to model loading of lakes or dams
 # number of surface loads
 0
-
diff -r 561a8af1a923 -r db5d0c6c7ab3 export.f90
--- a/export.f90	Fri Dec 30 17:19:40 2011 -0800
+++ b/export.f90	Thu Jan 05 00:27:28 2012 -0800
@@ -1077,8 +1077,9 @@ END SUBROUTINE exportcreep
   !!
   !! \author sylvain barbot 03/05/11 - original form
   !------------------------------------------------------------------
-  SUBROUTINE exportxy_rfaults(e,rffilename)
+  SUBROUTINE exportxy_rfaults(e,x0,y0,rffilename)
     TYPE(EVENT_STRUC), INTENT(IN) :: e
+    REAL*8, INTENT(IN) :: x0, y0
     CHARACTER(80), INTENT(IN) :: rffilename
 
     INTEGER :: iostatus,k
@@ -1109,8 +1110,8 @@ END SUBROUTINE exportcreep
        dip=e%s(k)%dip
 
        ! fault center position
-       x1=e%s(k)%x
-       x2=e%s(k)%y
+       x1=e%s(k)%x+x0
+       x2=e%s(k)%y+y0
        x3=e%s(k)%z
 
        ! fault dimension
diff -r 561a8af1a923 -r db5d0c6c7ab3 input.f90
--- a/input.f90	Fri Dec 30 17:19:40 2011 -0800
+++ b/input.f90	Thu Jan 05 00:27:28 2012 -0800
@@ -217,7 +217,7 @@ CONTAINS
        PRINT '("     * export to VTK format cancelled                     (--",a,")")', trim(opts(4)%name)
     END IF
     IF (in%isoutputvtkrelax) THEN
-       PRINT '("     * export relaxation component to VTK format          (--",a,")")', trim(opts(10)%name)
+       PRINT '("     * export relaxation component to VTK format   (--",a,")")', trim(opts(10)%name)
     END IF
 #endif
     PRINT 2000
@@ -1079,7 +1079,7 @@ CONTAINS
           CALL exportvtk_rfaults(in%events(e),rffilename)
 #endif
           rffilename=trim(in%wdir)//"/rfaults-"//digit//".xy"
-          CALL exportxy_rfaults(in%events(e),rffilename)
+          CALL exportxy_rfaults(in%events(e),in%x0,in%y0,rffilename)
 
           PRINT 2000
        END IF
diff -r 561a8af1a923 -r db5d0c6c7ab3 latex/graphics/cover.eps
Binary file latex/graphics/cover.eps has changed
diff -r 561a8af1a923 -r db5d0c6c7ab3 latex/relax.tex
--- a/latex/relax.tex	Fri Dec 30 17:19:40 2011 -0800
+++ b/latex/relax.tex	Thu Jan 05 00:27:28 2012 -0800
@@ -56,7 +56,7 @@
 
 \thispagestyle{empty}
 \begin{figure}[!h]
-\centering\includegraphics[width=15.1cm]{cover.eps}
+\centering\includegraphics[width=15.092cm]{cover.eps}
 \end{figure}
 \pagebreak
 
@@ -76,7 +76,7 @@ The possible applications for the earthq
 The possible applications for the earthquake-cycle modeling include i) co-seismic displacement and Coulomb stress calculation, ii) quasi-static stress transfer between earthquakes due to a postseismic transient, iii) modeling of postseismic transients including nonlinear rheologies and multiple mechanisms, iv) cycle of multiple earthquakes and spin-up models, v) loading cycle of lakes or the monsoon.
 
 \section{Acknowledgment}
-We are greatly thankful for the help of Yuri Fialko and Walter Landry, who contributed to the coming about of the software. We appreciate the efforts of Lucile Bruhat, Yaru Hsu, Mikhail Kogan and Baptiste Rousset, for testing an earlier version of the code. The support of CIG is greatly appreciated.
+We are greatly thankful for the help of Yuri Fialko and Walter Landry, who contributed to the coming about of the software. We appreciate the efforts of Lucile Bruhat, Yaru Hsu, Mikhail Kogan, Zhen Liu and Baptiste Rousset for testing an earlier version of the code. The support of CIG is greatly appreciated.
 
 
 \vspace{2.5cm}
@@ -226,7 +226,7 @@ In RELAX, the fault thickness is chosen 
 \section{Setting up the program}
 
 \subsection{Introduction}
-The RELAX code is written in Fortran90 with a few I/O functions written in C. The performance of the code depends greatly on the efficiency of the discrete Fourier transform being used. The program can work with the Cooley-Tukey FFT algorithm, for which the source code is provided. For better performance, it is recommended to use the FFT native to the computer environment. The program can readily use the SGI, the FFTW and the Intel MKL FFTs. While we have found that the Intel MLK FFT provides the most efficient calculation, the package provided by the CIG web site implements FFTW.
+The RELAX code is written in Fortran90 with a few I/O functions written in C. The performance of the code depends greatly on the efficiency of the discrete Fourier transform being used. The program can work with the Cooley-Tukey FFT algorithm, for which the source code is provided. For better performance, it is recommended to use the FFT native to the computer environment. The program can readily use the SGI, the FFTW and the Intel MKL FFTs. While we have found that the Intel MKL FFT provides the most efficient calculation, the package provided by the CIG web site implements FFTW.
 
 Both the post-processing and the storage of the simulation are greatly facilitated by writing output files in the cross-platform NetCDF binary format used by the Generic Mapping Tools (GMT). GMT is convenient to rapidly display the simulation results as it is computed, transform the output into other formats or projections (for example, to project the displacement into the Radar line of sight of a satellite to compare with synthetic aperture radar data), make animations and communicate results. Although RELAX can output in ASCII format, it is recommended to link the code to the GMT 4.5 libraries. A suite of GMT-based post-processing scripts are available in the \verb`util` directory and required GMT to be installed in your system.