Opening versus lengthening effects in moment rate as observed in far field body wave amplitudes

opening-lengthening-tradeoff.py

import numpy as np
import matplotlib.pyplot as plt

Lmax = 10.5e3

mustar = 3e9
tmax = 300

D_list = np.logspace(-3,1,100)
H_list = np.logspace(-1,np.log10(400),50)
umax = np.zeros((100,50))

for i,Dmax in enumerate(D_list):
    for j,H in enumerate(H_list):

        t = np.linspace(0,tmax)

        dLdt = Lmax/tmax
        dDdt = Dmax/tmax

        D = dDdt*t
        L = dLdt*t
        Mdot = mustar * H * ( dLdt * D + dDdt*L )

        r = 10e3
        alpha = 4000
        rho = 900

        u = Mdot / (4*np.pi*rho*alpha**3 * r)
        umax[i,j] = max(u)

c=plt.contourf(H_list,D_list,umax*1e6,cmap='Blues')
plt.contour(H_list,D_list,np.log10(umax*1e6),[-2,-1,0],colors='k')
plt.colorbar(c)
plt.yscale('log')
plt.xscale('log')
plt.ylabel('Final rift width')
plt.xlabel('Ice thickness')
plt.grid()
plt.show()