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

<p>% ==================================================================== %<br> % Data directory for 3D spherical, thermal convection models presented %<br> % in Fuchs and Becker (2022) with a strain-dependent weakening and &nbsp; &nbsp; %&nbsp;<br> % hardening rheology following the formulation of &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % Fuchs and Becker (2019). &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % The directory contains an input_files, MATLAB_Scripts, and each &nbsp; &nbsp; &nbsp;%&nbsp;<br> % model directory with a certain name. &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % &nbsp;&nbsp; &nbsp;input_file &nbsp;&nbsp; &nbsp; &nbsp;&nbsp; &nbsp;Contains all the input files for CitcomS &nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %&nbsp;<br> % &nbsp;&nbsp; &nbsp;MATLAB_Scripts &nbsp;&nbsp; &nbsp;Contains all the MATLAB scripts required to &nbsp; &nbsp;%<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;reproduce the figures in the manuscript &nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % &nbsp;&nbsp; &nbsp;$ModelName &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;Contains a MATLAB directory for individual &nbsp; &nbsp; %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;data from each model, &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;a TPR directory for toroidal-poloidal data for %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;each degree, and, &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;a txt_data for data picked from CitcomS models %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;which is then visualized in MATLAB. &nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;For the models discussed in the paper, surface %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;maps plots and surface grd-files are available %<br> % &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;within those directories as well. &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; %&nbsp;<br> % ==================================================================== %</p> <p>MATLAB_Scripts directory:&nbsp;<br> -------------------------<br> To visualize certain data for each model you can run the script&nbsp;<br> Analyze_Citcom_Models in the MALTAB_Scripts directory, e.g.,&nbsp;<br> &nbsp;&nbsp; &nbsp;<br> &nbsp;&nbsp; &nbsp;Analyze_Citcom_Models(Name,PlotParam,PlotParam2,S)</p> <p>where,&nbsp;<br> &nbsp;&nbsp; &nbsp;Name is the model name as a string variable,<br> &nbsp;&nbsp; &nbsp;PlotParam a switch to save (1) or not save (0) the figures,<br> &nbsp;&nbsp; &nbsp;PlotParam2 a switch to activate (1) or deactivate (0) plotting,&nbsp;&nbsp; &nbsp;<br> &nbsp;&nbsp; &nbsp;S a switch to define the scaling of the model parameter,<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;no scaling (0), scaling with the diffusion time scale (1),&nbsp;<br> &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;or scaling with the overturn time OT (2).</p> <p>With this script one can visualize all the time-dependent data picked<br> from the CitcomS models. The CitcomS models can be reproduce with the&nbsp;<br> input-files given in the input_files directory.&nbsp;</p> <p>To reproduce the box whisker plots in figures 2, 3, S7, and S8 one&nbsp;<br> needs to run the script CompStat. This scripts reads in the data&nbsp;<br> from all models (from the txt_files directory in the $ModelName&nbsp;<br> directory) and creates a box whisker plot for each model period and&nbsp;<br> plots them again the average lithospheric damage (gamma_L).</p> <p>For more details to each MATLAB script see the help comments within the<br> script.&nbsp;</p> <p>In case of any questions, do not hesitate to contact me via email:&nbsp;</p> <p>lukas.fuchs84 at gmail dot com</p> <p>% ==================================================================== %<br> % =============================== END ================================ %<br> % ==================================================================== %</p>