Every folder contains 5000 simulations. The parameters that define
our ensembles of simulations are specified in the folder name in the
following order:

 1) Plate geometry: Square / Circular.

 2) Plate-side length (S, in the case of a square plate) or plate-diameter length
    (D, in the case of a circular plate) in kilometers.
    e.g.) S550 (plate-side equal to 550 km); D770 (plate-diameter equal to 770 km).

 3) Viscosity of the asthenosphere (muA) in Pa*s 
    e.g.) muA_5e19: Viscosity of the asthenosphere is equal to 5*10^19 Pa*s.

 4) Viscosity of the lower part of the upper mantle (muM) in Pa*s.
    e.g.) muM_1p5e19: Viscosity of the lower part of the upper mantle is equal
    to 1.5*10^21 Pa*s.

 5) b-value of the Gutenberg-Richter scaling law.
    e.g.) b1p2: b-value is equal to 1.2.

 6) Amount of stress that slowly builds-up (charge phase) and suddenly drops 
    (discharge phase) that is transferred into a plate-motion (Euler-vector) change.
    Two possibilities: fs (full stress, all the stress is transferred into a plate-motion
    change); hs (half stress, only half of the stress is transferred into a plate-motion
    change).

Every simulation .mat file contains the following variables:

 EP: Initial microplate's Euler-pole.

 AV: Initial microplate's Angular Velocity.

 EV: Evolution of the microplate's Euler-vector through model time.
 We present five workspaces, depending on the earthquakes we account for:
   EV_allMw: we account for the Euler-vector changes caused by the charge
   and discharge of all the earthquakes in the set.
   EV_just_maxMw: we account for only the largest earthquake in the set.
   EV_up2Mw4p5: we account for the earthquakes whose magnitudes are smaller
   than 4.5 (Mw < 4.5).
   EV_up2Mw5p5: we account for the earthquakes whose magnitudes are smaller
   than 5.5 (Mw < 5.5).
   EV_up2Mw6p5: we account for the earthquakes whose magnitudes are smaller
   than 6.5 (Mw < 6.5), i.e. all the earthquakes but the largest one in the set.

Mw_DTB: Matrix with 9 columns that contains the magnitude, location and stress
associated to every earthquake in the simulation.

 1st col: Earthquake's moment magnitude (Mw).
 2nd col: Plate's side (only valid for a square plate, NaN for a circular plate).
 3rd col: Earthquake's rupture surface, in km^2.
 4th col: Earthquake's rupture length, in km.

The four following parameters define the location of the rupture area by using
spherical coordinates:
 5th col: Rupture's starting angle, in radians.
 6th col: Rupture's ending angle, in radians.
 7th col: Outer radius, in m
 8th col: Inner radius, in m.

 9th col: Stress drop, in Pa.