On the role of rheological memory for convection-driven plate reorganizations

% ==================================================================== %
% Data directory for 3D spherical, thermal convection models presented %
% in Fuchs and Becker (2022) with a strain-dependent weakening and     % 
% hardening rheology following the formulation of                   %
% Fuchs and Becker (2019).                                           %
%                                                  %
% The directory contains an input_files, MATLAB_Scripts, and each      % 
% model directory with a certain name.                           %
%     input_file          Contains all the input files for CitcomS          % 
%     MATLAB_Scripts     Contains all the MATLAB scripts required to    %
%                 reproduce the figures in the manuscript          %
%     $ModelName         Contains a MATLAB directory for individual     %
%                 data from each model,                      %
%                 a TPR directory for toroidal-poloidal data for %
%                 each degree, and,                      %
%                 a txt_data for data picked from CitcomS models %
%                 which is then visualized in MATLAB.          %
%                 For the models discussed in the paper, surface %
%                 maps plots and surface grd-files are available %
%                 within those directories as well.               % 
% ==================================================================== %

MATLAB_Scripts directory: 
-------------------------
To visualize certain data for each model you can run the script 
Analyze_Citcom_Models in the MALTAB_Scripts directory, e.g., 

    Analyze_Citcom_Models(Name,PlotParam,PlotParam2,S)

where, 
    Name is the model name as a string variable,
    PlotParam a switch to save (1) or not save (0) the figures,
    PlotParam2 a switch to activate (1) or deactivate (0) plotting,    
    S a switch to define the scaling of the model parameter,
        no scaling (0), scaling with the diffusion time scale (1), 
        or scaling with the overturn time OT (2).

With this script one can visualize all the time-dependent data picked
from the CitcomS models. The CitcomS models can be reproduce with the 
input-files given in the input_files directory. 

To reproduce the box whisker plots in figures 2, 3, S7, and S8 one 
needs to run the script CompStat. This scripts reads in the data 
from all models (from the txt_files directory in the $ModelName 
directory) and creates a box whisker plot for each model period and 
plots them again the average lithospheric damage (gamma_L).

For more details to each MATLAB script see the help comments within the
script. 

In case of any questions, do not hesitate to contact me via email: 

lukas.fuchs84 at gmail dot com

% ==================================================================== %
% =============================== END ================================ %
% ==================================================================== %