{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example 2: Lower crustal anisotropy" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In this example we generate P receiver functions for a model that includes a lower-crustal anisotropic layer. This example follows that of Figure 2 in [Porter et al. (2011)](#references), which uses the Raysum software developed by [Frederiksen and Bostock (2000)](#references])\n", "\n", "Start by importing the necessary modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from obspy.core import Stream\n", "from telewavesim import utils as ut\n", "from telewavesim import wiggle as wg" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Select the model file:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "modfile = '../models/model_Porter2011.txt'" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Select the type of incident wave - options are `'P'`, `'SV'`, `'SH'`, or `'Si'`, which is an isotropic S-wave source\n", "\n", "
\n", " Danger! Using 'SH' will not work properly for modeling receiver functions as the code will think you want plane-wave displacements (see below). Do not use 'SH' if you want S-wave receiver functions.\n", "
" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "wvtype = 'P'" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next we use variables to define the desired time series. \n", "\n", "