deletions | additions       diff --git a/untitled.tex b/untitled.tex  index f3b947d..ce9278b 100644  --- a/untitled.tex  +++ b/untitled.tex 
...
A nuclear explosion is approximated by time-dependent normal tractions applied on the surface of a 1 km-radius sphere carved out of the crust(Fig. \ref{fig:model_setup}c). The hollow sphere is at the center of the domain and buried at the depth of 2 km. The actual explosion must have occurred at a shallower depth but we used this geometry to prevent the mesh around the explosion site from being refined too much. The magnitude of normal traction is spatially uniform over the spherical boundary but varies in time. The magnitude is 0 Pa from 0 to 0.1 sec, linearly increases to 75 MPa from 0.1 to 2.1 sec and linearly decreases to 0 Pa from 2.1 to 4.1 sec. It remains at 0 Pa afterwards. The maximum magnitude and time variation of the tractions are adjusted such that model outputs can reasonably fit waveforms recorded at the two nearby seismic stations, DNH and YNB, in terms of the magnitude and duration. The monitoring points are set at the locations corresponding to the relative coordinates of these two stations with respect to the explosion site (Fig. \ref{fig:model_setup}b). DNH is at about 210 km NNW and YNB at about 160 km NNE from the test site. Displacement solutions interpolated at the monitoring points are compared with the records of the nuclear test-generated ground motions (Fig. \ref{fig:displ_mag}).  \subsection{Results}  The process of wave generation and propagation in the model with Vp/Vs equal to 1.7  is shown in Fig. \ref{fig:walkthrough}. Iso-magnitude surfaces for the displacement field reflects the main features included in the models such as the isotropic source, the layered structure composed of crust and mantle and the absorbing boundary conditions. As in the displacement magnitude plot (Fig. \ref{fig:displ_mag}), the dynamic pressure change is seen to attenuate slowly  Dynamic pressure change in the magma chamber, of which Vp/Vs is 1.7, is traced in terms of the mean and the maximum pressure in the chamber (Fig. \ref{fig:P_change_mean_max}). The peak values of the mean and the maximum dynamic pressure are 1.48 and 1.70 kPa and are reached at about 32 seconds since the initiation of the simulation, which corresponds to the arrival of surface waves(?). Post-peak fluctuation of the dynamic pressure change is persistent as expected from the lack of attenuation mechanism in the model.  Peak values of the dynamic pressure change are positively correlated with Vp/Vs ratios (Fig. \ref{fig:pressure_vs_VpVs}). The maximum increases by about 0.25 kPa from 1.65 kPa at Vp/Vs = 1.6 to 1.90 KPa at Vp/Vs = 2.0. Increase in the peak mean value is about 0.15 KPa: 1.45 kPa at Vp/Vs = 1.6 to 1.6 kPa at Vp/Vs = 2.0.   
...