# 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, ) # %%