Ground Motion Simulation of Moon Basin

As the initial trial of the project, Earthquake Ground Motion Modeling in Large Basins, the half moon basin was used to test our theories, methodologies, and techniques.
There are four cases done for a given input Ricker's pulse
Here are the figures we have been preparing for a paper to be published to BSSA later
Figure 1: geometry of the sythetic basin;
Figure 2:geometry and the input waves;
Figure 3: 3D mesh;
Figure 4: damping considerations;
Figure 5: computational results in time domain: 3 D and 1 D;
Figure 6: computational results in transfer functions: 3 D and 1 D;
Figure 7: computational results in response spectrums: 3 D and 1 D;
Figure 8: The time history of some nodes along axes from both BEM and FEM, case 1-a;
Figure 9: The time history of some nodes along axes for the other cases than case 1-a;
Figure 10: horizontal velocity distribution and history on plan surface;
Figure 11: transfer functions of displacement x-component (U) along y axis for the four cases;
Figure 12: transfer functions of U along x axis for the four cases;
Figure 13: transfer functions of U on the plane surface, f = 0.28Hz for the four cases;
Figure 14: transfer functions of U on the plane surface, f = 0.40Hz for the four cases;
Figure 15: transfer functions of U on the plane surface, f = 0.48Hz for the four cases;
Figure 16: response spectrums of U along y-axis for the four cases;
Figure 17: response spectrums of U along x-axis for the four cases;
Figure 18: response spectrum of U on the plane surface, fn = 0.25Hz for the four cases;
Figure 19: response spectrum of U on the plane surface, fn = 0.50Hz for the four cases;
Figure 20: response spectrum of U on the plane surface, fn = 0.75Hz for the four cases;
Figure 21: response spectrum of V on the plane surface, fn = 0.50Hz for the four cases;
Figure 22: maximum horizontal displacement on the plan surface;
Figure 23: maximum horizontal velocity on the plan surface;
Figure 24: maximum horizontal acceleration on the plan surface;
Figure 25: Fourier Transformation of U along y-axis;
Figure 26: Response spectrum for 2 locations along y for both 3D and 1D case.
Figure 27: Arias intensity of horizontal motion from the four cases;
Figure 28: Arias intensity of vertical motion from the four cases;
Figure 29: Local directivity of horizontal motion.