deletions | additions       diff --git a/untitled.tex b/untitled.tex  index 1c8e769..055d8a5 100644  --- a/untitled.tex  +++ b/untitled.tex 
...
%We numerically estimate pressure changes induced in the magma chamber by a nuclear explosion.  We use PyLith, ``a finite-element code for dynamic and quasistatic simulations of crustal deformation, primarily earthquakes and volcanoes'' \cite{Aagaard_2013} to compute pressure changes induced in the magma chamber by a nuclear explosion.  The region around Baekdu and the nuclear test site is represented by a 400 $\times$ 400 $\times$ 100 km domain. We adopt a velocity structure inferred for the region \cite{DUAN_2005} but for simplicity, assume that the entire domain is a single layer of crust with Vs = 3.25 km/s and Vp = 5.62 km/s. A magma chamber has been detected beneath Baekdu from low shear wave velocity in the depth range of 10-16 km \cite{ZHANG_2002}. Based on this observation, we assume that the magma chamber is a sphere of 3 km radius, of which center is at 13 km depth and 142 km away to the northeast of the center of the domain, where we put thenuclear  explosion site. source.  The S and P wave velocities assigned to the magma chamber are speeds of  3.0 km/s and 5.7 km/s. 5.19 km/s are assigned.  To represent anelastic attenuation in the crustal materials, materials and in the magma chamber,  the isotropic linear Maxwell viscoelastic rheology is assumed for the crust. assumed.  The constant viscosities of the crust and the magma chamber are 10$^{19}$ and 10$^{18}$ Pa$\cdot$s, respectively. The stresses are initially zero everywhere in the domain and the gravitational body force is not considered. These conditions are suitable to this study's purpose of estimating only dynamic stress changes. Absorbing boundary conditions are applied to all the boundaries except the top boundary, which is a free surface. 
...