rowanc1

import numpy as np
import matplotlib.pyplot as plt



def syntheticSeismogram(v, rho, d):
    """
    function syntheticSeismogram()

    syntheicSeismogram generates a synthetic seismogram for a simple 1-D

rowanc1

In this tutorial we will create a simple linear problem from scratch using the SimPEG framework.

rowanc1

There are two different forms of Richards equation that differ on how they deal with the non-linearity in the time-stepping term. Here we reproduce results from a seminal paper, and show the ease with which you can do this in SimPEG and Richards python packages.

rowanc1

%% Print Table
%
% This function will print an ASCII table, for example:
%
%              Decrease in Error, t_fo = 100
% ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
%                mu=0.1        mu=0.5        mu=1.0      
%           .--------------------------------------------
%    h=0.01 | 8.692922e-01  8.530311e-01  4.018080e-01   
%    h=0.02 | 5.797046e-01  6.220551e-01  7.596669e-02