[CIG-SEISMO] Anisotropy in SW4

Tue May 16 12:34:26 PDT 2017

Hi Lili,
Please provide your input file.
Thanks,
Anders

On 5/16/17, 12:00 PM, "CIG-SEISMO on behalf of cig-seismo-request at geodynamics.org" <cig-seismo-bounces at geodynamics.org on behalf of cig-seismo-request at geodynamics.org> wrote:

    Send CIG-SEISMO mailing list submissions to
    	cig-seismo at geodynamics.org

    To subscribe or unsubscribe via the World Wide Web, visit
    	http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo
    or, via email, send a message with subject or body 'help' to
    	cig-seismo-request at geodynamics.org

    You can reach the person managing the list at
    	cig-seismo-owner at geodynamics.org

    When replying, please edit your Subject line so it is more specific
    than "Re: Contents of CIG-SEISMO digest..."

    Today's Topics:

       1. Question about anisotropy for SW4 (Lili Feng)

    ----------------------------------------------------------------------

    Message: 1
    Date: Mon, 15 May 2017 14:54:45 -0600
    From: Lili Feng <Lili.Feng at colorado.edu>
    To: cig-seismo at geodynamics.org
    Subject: [CIG-SEISMO] Question about anisotropy for SW4
    Message-ID:
    	<CAD_j4bMkz_FAthbQV1etQhuri9U3jzqkSxT8Fuzvpgf2dr=zyA at mail.gmail.com>
    Content-Type: text/plain; charset="utf-8"

    Hi,
    I'm trying to run sw4 for anisotropic media, but the code just stuck and
    does not output anything.

    Here is the output information:
    "*----------------------------------------------------------------*
    *            sw4 version 1.1*

    * This program comes with ABSOLUTELY NO WARRANTY; released under GPL.*
    * This is free software, and you are welcome to redistribute     *
    * it under certain conditions, see LICENSE.txt for more details  *
    *----------------------------------------------------------------*
    *  Compiled on: Tue Apr  4 16:03:04 MDT 2017*
    *  By user:     life9360*
    *  Machine:     shas0137*
    *  Compiler:    /projects/life9360/anaconda2/bin/mpicxx*
    *  3rd party include dir: /include, and library dir: /lib*
    *----------------------------------------------------------------*

    *Input file: ak135_aniso_001.in <http://ak135_aniso_001.in>*"

    The code stuck here and did not output any more information...
    And below is my input file:

    *fileio path=ak135_aniso_001 pfs=1 verbose=1*
    *grid h=1000 nx=2001 ny=2001 nz=201 proj=utm ellps=WGS84*
    *time t=500*

    *source x=1000000 y=1000000 z=0 fz=1 t0=13.5 type=Dirac*
    *prefilter fc1=0.01 fc2=0.1*

    *rec x=1500000.000 y=1500000.000 depth=0.000 file=SW4.AAA sta=AAA
    writeEvery=1000 usgsformat=0 sacformat=1 nsew=1 variables=displacement*

    *anisotropy*
    *#------------------- ablock models -------------------*

    *ablock  z1=0.000000  z2=20000.000000  rho=2730.000000  c11=1.076668e+11
     c12=3.520350e+10  c13=3.458358e+10  c14=7.671698e+08  c22=1.033435e+11
     c23=3.484267e+10  c24=2.942917e+09  c33=9.909447e+10  c34=2.315423e+09
     c44=3.383794e+10  c55=3.397182e+10  c56=1.349769e+09  c66=3.506251e+10*

    *ablock  z1=20000.000000  z2=35000.000000  rho=2920.000000
     c11=1.233700e+11  c12=3.680660e+10  c13=3.680660e+10  c22=1.233700e+11
     c23=3.680660e+10  c33=1.233700e+11  c44=4.328170e+10  c55=4.328170e+10
     c66=4.328170e+10*

    *ablock  z1=35000.000000  z2=77000.000000  rho=3332.600000
     c11=2.155586e+11  c12=8.148659e+10  c13=8.148659e+10  c22=2.155586e+11
     c23=8.148659e+10  c33=2.155586e+11  c44=6.703600e+10  c55=6.703600e+10
     c66=6.703600e+10*

    *ablock  z1=77000.000000  z2=120000.000000  rho=3358.400000
     c11=2.174976e+11  c12=8.178445e+10  c13=8.178445e+10  c22=2.174976e+11
     c23=8.178445e+10  c33=2.174976e+11  c44=6.785656e+10  c55=6.785656e+10
     c66=6.785656e+10*

    *ablock  z1=120000.000000  z2=165000.000000  rho=3399.000000
     c11=2.236862e+11  c12=8.575123e+10  c13=8.575123e+10  c22=2.236862e+11
     c23=8.575123e+10  c33=2.236862e+11  c44=6.896748e+10  c55=6.896748e+10
     c66=6.896748e+10*

    *ablock  z1=165000.000000  z2=200000.000000  rho=3347.700000
     c11=2.271574e+11  c12=9.074873e+10  c13=9.074873e+10  c22=2.271574e+11
     c23=9.074873e+10  c33=2.271574e+11  c44=6.820432e+10  c55=6.820432e+10
     c66=6.820432e+10*

    I ran the code with 120 cores and waited for two days. The code is still
    running now without any warning or error messages, but no output at all...

    Lili

    -- 
    Leon(Lili) Feng
    Ph.D Candidate
    Center for Imaging the Earth's Interior
    Department of Physics
    University of Colorado Boulder
    +1-(720)-412-4379
    -------------- next part --------------
    An HTML attachment was scrubbed...
    URL: <http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20170515/238f4f09/attachment-0001.html>

    ------------------------------

    Subject: Digest Footer

    _______________________________________________
    CIG-SEISMO mailing list
    CIG-SEISMO at geodynamics.org
    http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo

    ------------------------------

    End of CIG-SEISMO Digest, Vol 112, Issue 4
    ******************************************