[CIG-SEISMO] Help with SW4 (Avinash Nayak)

Mon Apr 11 09:56:31 PDT 2016

Dear Anders,
Thanks for looking into this. I really appreciate your help. I have
attached my input file. I can try to send you my model (it is 1.2 GB)
and the videos too if you like. Please let me know if you see any
mistake.
Regards,
Avinash

On Mon, Apr 11, 2016 at 9:35 AM, Petersson, Anders <petersson1 at llnl.gov> wrote:
> Hi Avinash,
> Can you please attach your input file?
> Thanks,
> Anders
>
> On 4/10/16, 10:31 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. Help with SW4 (Avinash Nayak)
>>
>>
>>----------------------------------------------------------------------
>>
>>Message: 1
>>Date: Sun, 10 Apr 2016 22:31:02 -0700
>>From: Avinash Nayak <avinash07guddu at gmail.com>
>>To: cig-seismo at geodynamics.org
>>Subject: [CIG-SEISMO] Help with SW4
>>Message-ID:
>>       <CA+ebHScc3Sa3cYVrdjYE5=1_EiaxVf3qe-AWPQNf9_PG9QNjFQ at mail.gmail.com>
>>Content-Type: text/plain; charset="utf-8"
>>
>>Hi Everyone,
>>
>>I am having an issue with SW4 where I get non-zero values in velocity
>>at source position for a substantial amount of time (after source time
>>function becomes very small) and I am not sure what exactly I am doing
>>wrong.
>>
>>In the attached file (Summary.pdf), Figure 1 shows a screenshot of
>>velocity in YZ plane at x=1.393 km (passing through the source at
>>x=1.393, y=1.925 and z=0.469 km) at t=7.5 seconds. Color is supposed
>>to be scaled as red to blue from negative(peak absolute amplitude) to
>>positive(peak absolute amplitude) with white for zero. As you can see
>>the highest amplitude is a dot at the source position (y=1.925 km,
>>z=0.46 km). If I replace values at 7x7 grid points around the source
>>with zeros, then I can see the wavefield slightly better and the peak
>>amplitude decreases from 1.1e-5 m/s to 2.7e-6 m/s (Figure 2).
>>
>>In figure 3 and 4, I am plotting the velocity and displacement values
>>at the 7x7 grid points around the source, and you can see non-zero
>>very slowly decreasing values of velocity that persist even at 15-20
>>seconds. When I integrate them, I can see a slight increase of
>>displacements with time but I am not sure if they reach steady state
>>value (Figure 4).
>>
>>My grid spacing is 7 m for minimum shear wave velocity 150 m/s. I use
>>a Gaussian source time function with freq 3.1416 (so, f0=0.5 Hz,
>>fmax=f0 x 2.5 = 1.25 Hz, and 150/1.25/7 = 17 grid points per
>>wavelength for fmax). My t0 is 2.0 seconds( which is greater than
>>6/freq=1.91 s) so satisfies g(t=t0)=3.4e-9 which is less than 1e-8.
>>The source time function goes below 1e-15 after 5 seconds.
>>The source is a moment tensor source  very close to a strong material
>>discontinuity. I use supergrid size 73.
>>
>>I will be very grateful if someone can point out what I have done
>>wrong. Due to computation limitations, I cannot really go to lower
>>grid spacing. I am currently trying another run at higher value of t0
>>(2.5 s).
>>
>>Thanks.
>>Regards,
>>Avinash
>>
>>--
>>
>>Avinash Nayak,
>>PhD Candidate,
>>Berkeley Seismological Laboratory,
>>Department of Earth and Planetary Science,
>>University of California, Berkeley
>>-------------- next part --------------
>>A non-text attachment was scrubbed...
>>Name: Summary.pdf
>>Type: application/pdf
>>Size: 1461947 bytes
>>Desc: not available
>>URL:
>><http://lists.geodynamics.org/pipermail/cig-seismo/attachments/20160410/76
>>afc97a/attachment.pdf>
>>
>>------------------------------
>>
>>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 99, Issue 4
>>*****************************************
>
> _______________________________________________
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-seismo

-- 

Avinash Nayak,
PhD Candidate,
Berkeley Seismological Laboratory,
Department of Earth and Planetary Science,
University of California, Berkeley
-------------- next part --------------
fileio path=for1_te1_gau verbose=1 printcycle=500

grid nx=407 ny=483 nz=389 h=7.0 lat=29.9991677 lon=-91.1635866 az=0 proj=utm ellps=WGS84

time t=30.7

supergrid gp=73

source x=1393.0 y=1925.0 z=469.0 mxx=8.77489e10 myy=8.74982e10 mzz=8.50167e10 mxy=-5.45200e9 myz=5.82878e8 mxz=6.80339e9 type=Gaussian t0=2.0 freq=3.1416

#prefilter fc1=0.1 fc2=0.5 type=lowpass passes=1 order=4

pfile filename=../../../model_for2/model1.ppmod vsmin=151.0 vpmin=950.0 rhomin=1950.0 smoothingsize=2 style=cartesian

attenuation nmech=3 phasefreq=2.0 maxfreq=2.5

rec x=2045.2 y=2099.1 z=0.0 file=sta1 writeEvery=1000 
rec x=1974.4 y=803.5 z=0.0 file=sta2 writeEvery=1000 
rec x=900.9 y=2291.4 z=0.0 file=sta3 writeEvery=1000
rec x=1179.3 y=1684.7 z=0.0 file=sta8 writeEvery=1000 
rec x=1330.7 y=2473.5 z=0.0 file=sta9 writeEvery=1000 

image x=1393.0 mode=ux file=img_ux_x cycleInterval=75
image x=1393.0 mode=uy file=img_uy_x cycleInterval=75 
image x=1393.0 mode=uz file=img_uz_x cycleInterval=75 
image y=1925.0 mode=ux file=img_ux_y cycleInterval=75 
image y=1925.0 mode=uy file=img_uy_y cycleInterval=75 
image y=1925.0 mode=uz file=img_uz_y cycleInterval=75 
image z=469.0 mode=ux file=img_ux_z cycleInterval=75 
image z=469.0 mode=uy file=img_uy_z cycleInterval=75 
image z=469.0 mode=uz file=img_uz_z cycleInterval=75 

image z=0.0 mode=ux file=img_ux_z cycleInterval=75 
image z=0.0 mode=uy file=img_uy_z cycleInterval=75 
image z=0.0 mode=uz file=img_uz_z cycleInterval=75 

gmt file=final_config.gmt

#
#