
To use these files, you will need:
- authorized access to the iSALE shock physics code (iSALE-Dellen version) re-compiled with our modifications. Add reference to the manual in your work
- access to the pySALEPlot tool for iSALE users made by T. Davison. Add acknowledgement in your work
- running the iSALE models to generate the different jdata.dat files (average size of a jdata.file is 7 Go)
- python
- Ubuntu or macOS

++++++++++++++++++++++
 CONTENT OF DATA
++++++++++++++++++++++

Ready-to-use Aneos table (dunite_.aneos) and aneos input file (dunite_.input) to create aneos table for dunit

iSALE_mod.txt -  The iSALE-Dellen version used in this paper was slightly modified. Modifications are listed in the file "iSALE_mod.txt"
 		  Compile iSALE with these modifications before running the models.
 

Rubble_pile_model - contains configuration files for iSALE (asteroid.inp, additional.inp, material.inp, the BMP file of the sample)
                    note: the asteroid.inp represents an example for a specific setup (impactor size, impact velocity, boulder and matrix porosity)
                          depending on matrix porosity, you have to choose the bmp file   
		    
 			  asteroid.inp      : model general configuration
 			  material.inp      : model material properties
 			  additional.inp    : particule classes for specific host objects

Monolithic_model - contains configuration files for iSALE (asteroid.inp, material.inp)
                   note: the asteroid.inp represents an example for a specific setup (impactor size, impact velocity, bulk porosity of the target)

pysaleplot_scripts - contains following post-processing scripts:
 
PrePor.py		 - Pressure/Porosity contour plots
Compile.py		 - to create shock stage plots 
material_plot.py 	 - material plot with time for projectile, matrix, boulders
ShckMetMatFinalPos.py	 - pressure contour plot, final position of tracers
ShckPorMatFinalPos.py 	 - porosity contour plot, final position of tracers
ShckPorPreFinalPos.py 	 - porosity and pressure panel contour plot, final position of tracers
ShckPorPreInitialPos.py  - porosity and pressure panel contour plot, tracer position set back to initial position
ShockMetamorphism_mod.py - to create table with fraction and volumes of boulders, matrix, projectile for certain shock stages                          
EjektaInterpol_Ppeak_MY.py - Ejecta analysis , creates tables with detailed tracer analysis

note: change the path to pySALEPlot in the beginning of script, regarding your own directories
      the jdata.dat file from iSALE must be two folder down as the script reads it, so [script_folder]/xx/xx/jdata.dat


