Code in R for 'Several kilometers of global contraction on Mercury: A sample-size independent assessment of fault strain.'
Description
"Code for Mercury Contraction in R.txt" contains the code used to calculate the amount Mercury contracted using multiple thrust fault data with the R coding language. This is the methodology used in the paper titled 'Several kilometers of global contraction on Mercury: A sample-size independent assessment of fault strain.' which follows the works of Twiss and Marrett (2010a&b). To use, copy and paste the text in this file into your R command line.
Files
Steps to reproduce
This code serves as a framework to calculate the strain of a fault population across an entire planet's lithosphere. To use it copy and paste the code in an R command window. Set the directory at the beginning of the code. Note that values are populated with parameters used for the paper titled 'Several kilometers of global contraction on Mercury: A sample-size independent assessment of fault strain.' and should be changed based on the planetary body one is assessing. This code needs displacement-length and cumulative fault length frequency distribution samples in order to produce results. This code used data from three different data sets which are referenced as Byrne, Loveless, Watters throughout the code. These should be changed based on the fault data one wishes to use. The dimensions and displacement of the largest fault in the lithosphere also need to be constrained. These values have been populated with the dimensions of Enterprise Rupes on Mercury and will also need to be changed based on the planetary body one is assessing.