# Event configuration
ec = sn.EventConfiguration(waveform_simulation_configuration=wsc,
wavelet=sn.simple_config.stf.Ricker(
center_frequency=0.5 * max_frequency,
time_shift_in_seconds=0.1))
# Background model
bm = sn.model.background.homogeneous.IsotropicAcoustic(vp=vp[1], rho=2000.0)
# Simulation Configuration
p += sn.SimulationConfiguration(name="initial_model",
elements_per_wavelength=2,
tensor_order=2,
max_frequency_in_hertz=max_frequency,
model_configuration=sn.ModelConfiguration(
background_model=None,
volume_models="layered_model"),
absorbing_boundaries=absorbing_par,
event_configuration=ec)
p += sn.SimulationConfiguration(
name="target_model",
elements_per_wavelength=2,
tensor_order=2,
max_frequency_in_hertz=max_frequency,
model_configuration=sn.ModelConfiguration(background_model=bm),
absorbing_boundaries=absorbing_par,
event_configuration=sn.EventConfiguration(
waveform_simulation_configuration=wsc,
wavelet=sn.simple_config.stf.Ricker(center_frequency=0.5 *
max_frequency,
# salvus simulation object
sim = config.simulation.Waveform(mesh=mesh, sources=srcs, receivers=recs)
wsc = sn.WaveformSimulationConfiguration(end_time_in_seconds=9.5e-6)
wsc.physics.wave_equation.time_step_in_seconds = 1.0e-10
ec = sn.EventConfiguration(
waveform_simulation_configuration=wsc,
wavelet=sn.simple_config.stf.Ricker(center_frequency=20.0e6),
)
p += sn.SimulationConfiguration(
name="true_model_new_EM",
elements_per_wavelength=1.5,
tensor_order=4,
max_frequency_in_hertz=f_max,
model_configuration=sn.ModelConfiguration(background_model=None,
volume_models="true_model_EM"),
# Potentially event dependent settings.
event_configuration=ec,
)
# %%
# Modeling
p.simulations.launch(
simulation_configuration="true_model_new_EM",
events=p.events.get_all(),
site_name=SALVUS_FLOW_SITE_NAME,
ranks_per_job=28,
wall_time_in_seconds_per_job=10,
)
# %%