Archive of the iSALE simulation datasets for Jackson et al. 2022, "Impact generation of holes in the early lunar crust I: scaling relations"

Archive of the iSALE simulation datasets used in Jackson et al. 2022, "Impact generation of holes in the early lunar crust I: scaling relations".

This archive contains all necessary information to re-run any of the simulations in the dataset along with plots that can be used to check the outcomes of each simulation.  In addition, it also contains a full compilation of all of the derived data used in Jackson et al. 2022 as well as the analysis scripts used to extract that data from the simulations and the scripts that use the derived data to generate the figures in the manuscript.

The archive is provided in a standard compressed (.zip) format.

Detailed description of the contents

Note that this description can also be found in the ReadMe file provided with the archive

1) Simulation folders

The folder 'simsmain' contains files for the iSALE2D simulations that comprise the primary simulation database.

For each simulation we provide:

• asteroid.inp - the initialisation file for the simulation that would allow it to be re-run
• [].dat - a text file containing data derived from the final simulation output.  These values from all simulations are compiled in alldata.csv.  A description of the columns is provided below for alldata.csv
• []-Crthickprof.pdf - a figure showing the radial crust thickness profile at the end of the simulation. Light pink shows the cell-by-cell profile, solid red a moving average.
• []-Damprof.pdf - a figure showing the radial profile of the mean damage factor (averaged over the depth of the crust at each radial cell) at the end of the simulation.  Light grey shows the cell-by-cell profile, solid black a moving average.
• []-TmpDam.dpdf - a figure showing the temperature and damage structure of the entire high-resolution zone at the end of the simulation at 1:1 aspect ratio (i.e. square cells are square in the image).  The boundary of the crust is outlined in black.
• []-TmpDam-zoom.pdf - as above, but zoomed in on the central region of interest in the simulation.  This image is stretched to highlight the region of interest (i.e. square cells are not square in the image).
• []-temporal.dat - a text file containing data derived from the full simulation time series.  These values from all simulations are compiled in alldata_temporal.csv.  A description of the columns is provided below for alldata_temporal.csv
• []-crdepthtime.png - a figure showing the change in the maximum depression of the upper surface of the simulation over time.  This records the excavation and collapse of the initial transient cavity.
• []-holeradtime.png - a figure recording the radius of the hole (as measured by the locations at which the crust thickness passes 90 and 95 percent of the initial value) over time.  A final plateau is indicative of the simulation having settled.
• []-holevoltime.png - a figure recording the volume of the hole within the 95% crust thickness radius over time.
• []-trcavtime.png - a figure recording the volume of the surface depression over time.  This can be quite sensitive to small fluctuations at the edges of the surface depression and does not always plateau.

Within simsmain each sub-folder is named in the format dW_vX_cY_gZ, where d is the diameter of the impactor (in km), v is the impact velocity (in m/s), c is the crust thickness (in km) and g is the grid resolution of the high-resolution zone (in km).  W, X, Y and Z are all in exponential notation with two significant figures (i.e. N.NE+NN).

The folder 'simsres' contains files for the additional iSALE2D simulations run for resolution testing purposes.  These files have the same formats and naming conventions as simsmain.

The folder 'material' contains the file material.inp, the initialisation file for the material properties which is common to all of the simulations.

2) Analysis scripts

The folder 'analysisscripts' contains the Python scripts that were used to generate the files in simsmain and simsres.  There are two scripts, iSALEanalysis_batch.py which generates the end-of-simulation plots and data, and iSALEanalysis_batch_temporal.py which generates the time series plots and data.  The temporal script typically takes around 400 times longer to run than the end-of-simulation script since the former must examine all (usually 800) timesteps, whereas the end-of-simulation script examines only the first and last timestep.

3) Figures

The folder 'figures' contains CSV files summarising all of the data derived from the iSALE2D simulations along with plotting scripts that use this derived data to generate the plots in Jackson et al. 2022.

-----------
alldata.csv is a compilation of the data in the [].dat files for all 252 simulations in the primary simulation database.  There are 14 columns which in order are:

1. dimp        - impactor diameter in km
2. vimp        - impact velocity in km/s
3. crd         - initial crust thickness in km
4. tdat        - simulation time (in seconds) at which the data was taken
5. gres        - grid resolution in the high-resolution zone in km
6. rhires      - radius of the high-resolution zone in km
7. r_thincr0_95- radial distance (in km) at which the crust thickness exceeds 95% of the initial value
8. r_thincr0_9 - radial distance (in km) at which the crust thickness exceeds 90% of the initial value
9. r_magcr     - radius (in km) within which there is direct exposure of magma to the surface (i.e. crust thickness is zero)
10. r_frac09    - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.9
11. r_frac08    - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.8
12. r_frac05    - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.5
13. holevolume  - volume (in km^3) of crust excavated from within r_thincr0_95
14. comment2    - impact regime classification, 'partial'= partial penetration, 'complete' = complete penetration, 'cfd' = cratering with full-depth fracturing, 'classical' = classical crater.

-----------
alldata_temporal.csv is a compilation of the data in the []-temporal.dat files for all 252 simulations in the primary simulation database.  There are 20 columns which in order are:

1. dimp        - impactor diameter in km
2. vimp        - impact velocity in m/s (Note the difference from alldata.csv)
3. crd         - initial crust thickness in km
4. tdat        - simulation time (in seconds) at which the data was taken
5. gres        - grid resolution in the high-resolution zone in km
6. rhires      - radius of the high-resolution zone in km
7. r_thincr0_95- radial distance (in km) at which the crust thickness exceeds 95% of the initial value
8. r_thincr0_9 - radial distance (in km) at which the crust thickness exceeds 90% of the initial value
9. r_magcr     - radius (in km) within which there is direct exposure of magma to the surface (i.e. crust thickness is zero)
10. r_frac09    - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.9
11. r_frac08    - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.8
12. r_frac05    - radial distance (in km) within which the vertically averaged damage factor in the crust exceeds 0.5
13. holevolume  - volume (in km^3) of crust excavated from within r_thincr0_95
14. mcrd        - maximum depth of the transient cavity in km
15. t_mcrd      - simulation time (in seconds) at which transient cavity reaches maximum depth
16. trcav_mcrd  - volume (in km^3) of the transient cavity at the time of maximum depth
17. mtrcav      - maximum volume of the transient cavity
18. t_mtrcav    - simulation time (in seconds) at which transient cavity reaches maximum volume
19. crd_mtrcav  - depth of the transient cavity (in km) at the time of maximum volume
20. comment2    - impact regime classification, 'partial'= partial penetration, 'complete' = complete penetration, 'cfd' = cratering with full-depth fracturing, 'classical' = classical crater.

-----------
restest-rthin.csv and restest-frac.csv are in the same format as alldata.csv and contain the data for the sets of simulations used to test the resolution dependence of r_thin and r_frac respectively.

-----------

There are 7 Python plotting scripts:

1. plotter_main.py - the main plotting script, produces plots that don't fall under another heading
2. plotter_KE.py   - produces plots as a function of impact kinetic energy
3. plotter_temporal.py - produces plots using time series data from alldata_temporal
4. plotter_cox.py  - produces plots comparing our results to those of Bray et al. 2014, Cox et al. 2008, Cox & Bauer 2015 and Miljkovic et al. 2015.  Note that we do not include the CSV files containing the data from those papers since it is not ours to distribute, however they can be easily extracted from the relevant publications.
5. plotter_res.py  - produces resolution comparison plots
6. plotter_pi.py   - produces Pi-scaling plots
7. kplot-sims.py   - produces a plot of the 'k' factor (hole size/impactor mass) for comparison with Perera et al. 2018

All of the plots and fitting summaries produced by the plotting scripts are also included