/[escript]/trunk/doc/examples/inversion/synthetic_VTI.py
ViewVC logotype

Diff of /trunk/doc/examples/inversion/synthetic_VTI.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 4621 by gross, Tue Dec 17 23:37:08 2013 UTC revision 4622 by sshaw, Fri Jan 17 04:55:41 2014 UTC
# Line 22  __url__="https://launchpad.net/escript-f Line 22  __url__="https://launchpad.net/escript-f
22  from esys.escript import *  from esys.escript import *
23  from esys.escript import unitsSI as U  from esys.escript import unitsSI as U
24  from esys.escript.pdetools import Locator  from esys.escript.pdetools import Locator
25  from esys.finley import Brick, Rectangle  from esys.ripley import Brick, Rectangle
26  from esys.weipa import saveSilo  from esys.weipa import saveSilo
27  from esys.downunder import Ricker, VTIWave, SimpleSEGYWriter  from esys.downunder import Ricker, VTIWave, SimpleSEGYWriter
28  from math import ceil  from math import ceil
# Line 117  delta=0.05*0 Line 117  delta=0.05*0
117  rho=2000*U.kg/U.m**3  rho=2000*U.kg/U.m**3
118  src_dir=[0,0,1]  src_dir=[0,0,1]
119    
120  sw=VTIWave(domain, v_p, v_s, wl, src_tags[0], source_vector = src_dir, eps=vareps, gamma=gamma, delta=delta, rho=rho,  \  sw=VTIWave(domain, v_p, v_s, wl, src_tags[0], source_vector = src_dir
121                       absorption_zone=300*U.m, absorption_cut=1e-2, lumping=True)                  eps=vareps, gamma=gamma, delta=delta, rho=rho,
122                    absorption_zone=300*U.m, absorption_cut=1e-2, lumping=True)
123    
124  locEW=Locator(domain,rcvEW_locations)  locEW=Locator(domain,rcvEW_locations)
125  tracerEW_x=SimpleSEGYWriter(receiver_group=rgEW, source=src_loc_2D, sampling_interval=sampling_interval, text='x-displacement - east-west line')  tracerEW_x=SimpleSEGYWriter(receiver_group=rgEW, source=src_loc_2D, sampling_interval=sampling_interval, text='x-displacement - east-west line')

Legend:
Removed from v.4621  
changed lines
  Added in v.4622

  ViewVC Help
Powered by ViewVC 1.1.26