Skip to content

Benchmarks for the effects of ocean. Compared between global solver SPECFEM3D Globe, regional solver RegSEM, and normal-mode calculations

Notifications You must be signed in to change notification settings

earthinversion/Compare-Synthetics-for-Earth-Models-Ocean-Effects

Repository files navigation

How to run?

  1. Edit the input_params.yml
  2. Run the script check_traces_single_station.py
    python check_traces_single_station.py
    

Dependencies

  1. Numpy
  2. Matplotlib
  3. Obspy
  4. Pyyaml
  5. Pandas

Tested Cases

  • Compared the Specfem Layered with Polynomial: matches fairly well
  • Compared the Specfem Layered with ocean vs Specfem layered without ocean
  • Compared RegSEM with Normal-mode synthetics
  • Compared RegSEM with SPECFEM for no ocean, no topography: matches very well
  • Compared the SPECFEM runs for R_EARTH = 6371 and 6368 km
  • Compared SPECFEM with/without elevated receivers
  • Compared SPECFEM with SPECFEM version with etopo values 0

Input parameters

Model: Preliminary Earth Reference Model, with the ocean layer of thickness 3 km

  • Spherical Earth: ON
  • Anisotropy: ON
  • Oceans: ON
  • Topography: ON
  • Gravity: OFF
  • Attenuation: OFF
  • Ellipticity: OFF

Input Event and Stations Locations


Event-Stations Map

Event-Stations Map


Results

  • Comparison between SPECFEM, RegSEM and NMS (Normal-mode synthetics) with OCEAN, anisotropy and Topography. No mirror for the cases of SPECFEM and RegSEM.

BBB Time

BBB Time


BBB Frequency

BBB Frequency


About

Benchmarks for the effects of ocean. Compared between global solver SPECFEM3D Globe, regional solver RegSEM, and normal-mode calculations

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published