============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 1: 7:58 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174373E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984348 2 1.416236E-09 3.1749E-09 0.397013 0.29 7.0003E+13 0.907625 3 1.439403E-09 1.8516E-09 0.091336 0.73 4.3943E+13 0.696277 4 1.450008E-09 1.3739E-09 0.036102 0.90 5.0291E+13 0.419761 5 1.455032E-09 8.6402E-10 0.017929 0.96 2.2904E+13 0.257770 6 1.457713E-09 7.8453E-10 0.009790 0.98 1.5983E+13 0.176929 7 1.459204E-09 4.9423E-10 0.006334 0.98 1.5046E+13 0.130361 8 1.460097E-09 4.1571E-10 0.004674 0.99 3.2232E+13 0.100547 9 1.460696E-09 3.9804E-10 0.003638 0.99 1.2557E+13 0.079348 10 1.461066E-09 2.8661E-10 0.002978 0.99 4.1745E+13 0.066351 11 1.461335E-09 2.1282E-10 0.002466 0.98 7.8352E+12 0.053579 12 1.461555E-09 1.6792E-10 0.002075 1.00 7.2000E+12 0.045053 13 1.461725E-09 1.4184E-10 0.001785 1.00 7.1228E+12 0.038364 14 1.461869E-09 1.2246E-10 0.001554 1.00 6.2440E+12 0.033071 15 1.461990E-09 1.0590E-10 0.001357 1.00 6.1823E+12 0.028691 16 1.462100E-09 9.1816E-11 0.001193 1.00 6.1307E+12 0.025000 17 1.462199E-09 7.8635E-11 0.001052 1.00 5.6715E+12 0.021915 18 1.462292E-09 6.7667E-11 0.000927 1.00 5.0936E+12 0.019367 19 1.462377E-09 5.8515E-11 0.000820 1.00 3.8106E+12 0.017022 20 1.462453E-09 5.0551E-11 0.000727 1.00 3.5117E+12 0.015064 21 1.462523E-09 4.3654E-11 0.000647 1.00 3.2034E+12 0.013390 22 1.462587E-09 3.7710E-11 0.000576 1.00 2.9060E+12 0.011927 23 1.462645E-09 3.1988E-11 0.000513 1.00 2.6054E+12 0.010639 24 1.462699E-09 2.7722E-11 0.000458 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 1:20:53 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 1:20:53 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174395E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984357 2 1.420107E-09 3.1756E-09 0.396820 0.30 6.9387E+13 0.907646 3 1.443556E-09 2.0153E-09 0.089909 0.72 4.4100E+13 0.695969 4 1.453379E-09 1.3684E-09 0.034766 0.90 4.8830E+13 0.420646 5 1.457498E-09 8.5281E-10 0.016981 0.96 2.0152E+13 0.257381 6 1.459636E-09 7.9021E-10 0.009459 0.97 1.5526E+13 0.176956 7 1.460800E-09 4.8971E-10 0.006175 0.99 1.3686E+13 0.129402 8 1.461508E-09 4.0512E-10 0.004507 0.99 1.5478E+13 0.099584 9 1.461980E-09 3.8401E-10 0.003517 1.00 1.4131E+13 0.078877 10 1.462255E-09 2.8596E-10 0.002865 0.99 4.9883E+13 0.067581 11 1.462399E-09 2.1038E-10 0.002435 0.95 4.6640E+13 0.057035 12 1.462542E-09 1.6391E-10 0.001977 0.99 8.2748E+12 0.045379 13 1.462642E-09 1.3471E-10 0.001680 1.00 1.1835E+13 0.038396 14 1.462726E-09 1.1403E-10 0.001450 1.00 1.1079E+13 0.033144 15 1.462798E-09 9.7671E-11 0.001260 1.00 9.0371E+12 0.028746 16 1.462863E-09 8.2871E-11 0.001105 1.00 7.7888E+12 0.024973 17 1.462920E-09 7.1100E-11 0.000976 1.00 4.7103E+12 0.021613 18 1.462976E-09 6.1085E-11 0.000866 1.00 4.6394E+12 0.018998 19 1.463027E-09 5.2431E-11 0.000767 1.00 4.0068E+12 0.016837 20 1.463073E-09 4.4823E-11 0.000683 1.00 3.3044E+12 0.014979 21 1.463115E-09 3.8659E-11 0.000609 1.00 2.8650E+12 0.013332 22 1.463154E-09 3.2806E-11 0.000546 1.00 2.4919E+12 0.011952 23 1.463189E-09 2.8541E-11 0.000491 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 1:34:30 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 1:34:30 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174402E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984358 2 1.420478E-09 3.1743E-09 0.396702 0.30 6.9268E+13 0.907508 3 1.443984E-09 2.0059E-09 0.089659 0.72 4.5775E+13 0.694675 4 1.453748E-09 1.3690E-09 0.034427 0.90 4.8778E+13 0.419575 5 1.457876E-09 8.7028E-10 0.016847 0.96 2.0204E+13 0.256626 6 1.459997E-09 7.7788E-10 0.009408 0.97 1.5621E+13 0.176173 7 1.461159E-09 4.9483E-10 0.006144 0.99 1.3385E+13 0.128960 8 1.461860E-09 4.1615E-10 0.004485 0.99 1.3800E+13 0.099006 9 1.462336E-09 3.9547E-10 0.003489 1.00 1.3104E+13 0.078619 10 1.462595E-09 2.8504E-10 0.002823 1.00 3.8429E+13 0.065443 11 1.462790E-09 2.1148E-10 0.002293 0.99 1.1355E+13 0.053350 12 1.462902E-09 1.6739E-10 0.001946 1.00 8.8662E+13 0.047847 13 1.463023E-09 1.4075E-10 0.001622 0.98 1.0743E+13 0.038285 14 1.463070E-09 1.2059E-10 0.001470 0.99 4.1467E+13 0.034690 15 1.463142E-09 1.0398E-10 0.001234 1.00 8.5775E+12 0.028888 16 1.463207E-09 9.0040E-11 0.001076 1.00 9.5205E+12 0.024800 17 1.463272E-09 7.7665E-11 0.000946 1.00 4.5529E+12 0.021682 18 1.463327E-09 6.6813E-11 0.000838 1.00 4.2068E+12 0.019093 19 1.463377E-09 5.7290E-11 0.000741 1.00 3.5807E+12 0.016850 20 1.463421E-09 4.9448E-11 0.000659 1.00 2.8217E+12 0.014906 21 1.463461E-09 4.2603E-11 0.000589 1.00 2.6716E+12 0.013229 22 1.463497E-09 3.6802E-11 0.000528 1.00 2.5188E+12 0.011761 23 1.463529E-09 3.1851E-11 0.000475 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 1:48: 8 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 1:48: 8 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174403E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984359 2 1.420496E-09 3.1743E-09 0.396720 0.30 6.9234E+13 0.907440 3 1.444056E-09 2.0051E-09 0.089581 0.72 4.5904E+13 0.694202 4 1.453861E-09 1.3690E-09 0.034353 0.89 4.8731E+13 0.419348 5 1.458027E-09 8.7186E-10 0.016820 0.96 2.0797E+13 0.256462 6 1.460175E-09 7.7674E-10 0.009404 0.97 1.5671E+13 0.175998 7 1.461366E-09 4.9537E-10 0.006128 0.99 1.3367E+13 0.128823 8 1.462083E-09 4.1794E-10 0.004471 0.99 1.3437E+13 0.098935 9 1.462572E-09 3.9854E-10 0.003471 0.99 1.3022E+13 0.078621 10 1.462841E-09 2.8555E-10 0.002802 1.00 3.8428E+13 0.065407 11 1.463042E-09 2.1257E-10 0.002274 0.98 1.1345E+13 0.053358 12 1.463194E-09 1.6904E-10 0.001877 1.00 1.0748E+13 0.045010 13 1.463304E-09 1.4241E-10 0.001602 1.00 1.0324E+13 0.037939 14 1.463371E-09 1.2304E-10 0.001415 1.00 3.8041E+13 0.034004 15 1.463439E-09 1.0694E-10 0.001216 1.00 8.5392E+12 0.028705 16 1.463498E-09 9.2536E-11 0.001059 1.00 9.4443E+12 0.024735 17 1.463551E-09 8.0046E-11 0.000930 1.00 4.4514E+12 0.021568 18 1.463599E-09 6.9001E-11 0.000821 1.00 3.8703E+12 0.018996 19 1.463643E-09 5.9026E-11 0.000726 1.00 3.4936E+12 0.016799 20 1.463682E-09 5.0494E-11 0.000645 1.00 2.8011E+12 0.014865 21 1.463717E-09 4.3376E-11 0.000575 1.00 2.6720E+12 0.013241 22 1.463749E-09 3.7032E-11 0.000515 1.00 2.5196E+12 0.011774 23 1.463777E-09 3.2038E-11 0.000462 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 2: 1:47 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 2: 1:47 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174404E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984359 2 1.420486E-09 3.1743E-09 0.396772 0.30 6.9224E+13 0.907420 3 1.444079E-09 2.0051E-09 0.089544 0.72 4.5818E+13 0.694100 4 1.453919E-09 1.3690E-09 0.034320 0.89 4.8711E+13 0.419370 5 1.458116E-09 8.7039E-10 0.016827 0.96 2.0789E+13 0.256268 6 1.460283E-09 7.7782E-10 0.009400 0.97 1.5687E+13 0.176035 7 1.461493E-09 4.9534E-10 0.006123 0.99 1.3389E+13 0.128945 8 1.462221E-09 4.1765E-10 0.004461 0.99 1.3250E+13 0.098976 9 1.462717E-09 3.9886E-10 0.003454 0.99 1.3151E+13 0.078660 10 1.462990E-09 2.8563E-10 0.002784 1.00 3.8302E+13 0.065459 11 1.463190E-09 2.1270E-10 0.002254 0.98 1.1346E+13 0.053356 12 1.463340E-09 1.6875E-10 0.001858 1.00 1.0703E+13 0.045012 13 1.463445E-09 1.4269E-10 0.001590 1.00 1.0243E+13 0.038117 14 1.463507E-09 1.2332E-10 0.001415 1.00 3.5566E+13 0.034049 15 1.463577E-09 1.0713E-10 0.001196 1.00 8.5232E+12 0.028669 16 1.463636E-09 9.2605E-11 0.001040 1.00 9.4516E+12 0.024686 17 1.463689E-09 8.0108E-11 0.000915 1.00 4.4005E+12 0.021542 18 1.463737E-09 6.8436E-11 0.000807 1.00 3.8753E+12 0.018963 19 1.463782E-09 5.8480E-11 0.000713 1.00 3.4163E+12 0.016768 20 1.463821E-09 5.0497E-11 0.000634 1.00 2.7936E+12 0.014860 21 1.463856E-09 4.3586E-11 0.000565 1.00 2.6662E+12 0.013228 22 1.463888E-09 3.7641E-11 0.000505 1.00 2.5151E+12 0.011786 23 1.463916E-09 3.2542E-11 0.000453 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 2:15:22 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 2:15:22 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174404E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984359 2 1.420495E-09 3.1743E-09 0.396767 0.30 6.9221E+13 0.907412 3 1.444105E-09 2.0051E-09 0.089509 0.72 4.5884E+13 0.694071 4 1.453962E-09 1.3690E-09 0.034301 0.89 4.8705E+13 0.419362 5 1.458169E-09 8.6895E-10 0.016820 0.96 2.0786E+13 0.256225 6 1.460344E-09 7.7962E-10 0.009399 0.97 1.5691E+13 0.176057 7 1.461559E-09 4.9480E-10 0.006120 0.99 1.3364E+13 0.128925 8 1.462291E-09 4.1676E-10 0.004455 0.99 1.3227E+13 0.098945 9 1.462789E-09 3.9823E-10 0.003445 0.99 1.3145E+13 0.078597 10 1.463062E-09 2.8569E-10 0.002773 1.00 3.8273E+13 0.065438 11 1.463263E-09 2.1263E-10 0.002246 0.98 1.1339E+13 0.053334 12 1.463414E-09 1.6843E-10 0.001849 1.00 1.0692E+13 0.044990 13 1.463523E-09 1.4226E-10 0.001577 1.00 1.0166E+13 0.037971 14 1.463582E-09 1.2287E-10 0.001411 1.00 4.7402E+13 0.034351 15 1.463653E-09 1.0672E-10 0.001187 1.00 8.3993E+12 0.028669 16 1.463712E-09 9.2210E-11 0.001034 1.00 9.4117E+12 0.024661 17 1.463765E-09 7.9943E-11 0.000908 1.00 4.3942E+12 0.021524 18 1.463814E-09 6.8659E-11 0.000801 1.00 3.8753E+12 0.018932 19 1.463859E-09 5.8564E-11 0.000708 1.00 3.3974E+12 0.016763 20 1.463897E-09 5.0093E-11 0.000628 1.00 2.7942E+12 0.014853 21 1.463932E-09 4.3252E-11 0.000560 1.00 2.6666E+12 0.013214 22 1.463964E-09 3.7370E-11 0.000500 1.00 2.5153E+12 0.011769 23 1.463992E-09 3.2322E-11 0.000449 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 2:28:55 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 2:28:55 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174404E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984359 2 1.420498E-09 3.1743E-09 0.396762 0.30 6.9221E+13 0.907395 3 1.444124E-09 2.0051E-09 0.089481 0.72 4.5917E+13 0.694066 4 1.453993E-09 1.3690E-09 0.034274 0.89 4.8683E+13 0.419382 5 1.458206E-09 8.6775E-10 0.016824 0.96 2.0790E+13 0.256287 6 1.460386E-09 7.8080E-10 0.009402 0.97 1.5694E+13 0.176117 7 1.461600E-09 4.9439E-10 0.006134 0.99 1.3369E+13 0.128954 8 1.462337E-09 4.1605E-10 0.004444 0.99 1.3286E+13 0.098930 9 1.462836E-09 3.9795E-10 0.003439 0.99 1.3133E+13 0.078697 10 1.463110E-09 2.8600E-10 0.002767 1.00 3.8175E+13 0.065473 11 1.463313E-09 2.1238E-10 0.002242 0.98 1.1322E+13 0.053378 12 1.463466E-09 1.6810E-10 0.001845 1.00 1.0723E+13 0.045052 13 1.463576E-09 1.4184E-10 0.001574 1.00 1.0161E+13 0.037950 14 1.463641E-09 1.2243E-10 0.001397 1.00 4.5509E+13 0.034183 15 1.463711E-09 1.0636E-10 0.001187 1.00 8.4277E+12 0.028635 16 1.463770E-09 9.1945E-11 0.001031 1.00 9.3952E+12 0.024643 17 1.463823E-09 7.9623E-11 0.000905 1.00 4.4055E+12 0.021521 18 1.463871E-09 6.8732E-11 0.000798 1.00 3.8752E+12 0.018945 19 1.463916E-09 5.8609E-11 0.000705 1.00 3.3742E+12 0.016752 20 1.463955E-09 4.9982E-11 0.000625 1.00 2.7971E+12 0.014840 21 1.463989E-09 4.2880E-11 0.000557 1.00 2.6688E+12 0.013225 22 1.464021E-09 3.7056E-11 0.000498 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 2:41:58 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 2:41:58 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174405E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984359 2 1.420508E-09 3.1743E-09 0.396768 0.30 6.9221E+13 0.907385 3 1.444152E-09 2.0052E-09 0.089456 0.72 4.5942E+13 0.694057 4 1.454026E-09 1.3690E-09 0.034268 0.89 4.8678E+13 0.419445 5 1.458236E-09 8.6734E-10 0.016830 0.96 2.0787E+13 0.256313 6 1.460420E-09 7.8104E-10 0.009396 0.97 1.5691E+13 0.176090 7 1.461637E-09 4.9424E-10 0.006133 0.99 1.3358E+13 0.128996 8 1.462375E-09 4.1577E-10 0.004443 0.99 1.3293E+13 0.098942 9 1.462873E-09 3.9773E-10 0.003438 0.99 1.3126E+13 0.078711 10 1.463148E-09 2.8602E-10 0.002765 1.00 3.8165E+13 0.065484 11 1.463352E-09 2.1260E-10 0.002241 0.98 1.1322E+13 0.053368 12 1.463507E-09 1.6794E-10 0.001847 1.00 1.0705E+13 0.045057 13 1.463619E-09 1.4166E-10 0.001573 1.00 1.0160E+13 0.037946 14 1.463686E-09 1.2225E-10 0.001389 1.00 4.1722E+13 0.034071 15 1.463757E-09 1.0618E-10 0.001188 1.00 8.4520E+12 0.028622 16 1.463816E-09 9.1820E-11 0.001032 1.00 9.3848E+12 0.024642 17 1.463870E-09 7.9474E-11 0.000905 1.00 4.4096E+12 0.021513 18 1.463919E-09 6.8734E-11 0.000797 1.00 3.8777E+12 0.018925 19 1.463964E-09 5.8475E-11 0.000704 1.00 3.3765E+12 0.016755 20 1.464003E-09 5.0156E-11 0.000625 1.00 2.7987E+12 0.014849 21 1.464038E-09 4.2813E-11 0.000556 1.00 2.6701E+12 0.013230 22 1.464070E-09 3.6924E-11 0.000497 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 2:54:58 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 2:54:58 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174405E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984359 2 1.420504E-09 3.1743E-09 0.396764 0.30 6.9221E+13 0.907373 3 1.444158E-09 2.0052E-09 0.089440 0.72 4.5956E+13 0.694052 4 1.454041E-09 1.3690E-09 0.034270 0.89 4.8674E+13 0.419435 5 1.458255E-09 8.6694E-10 0.016826 0.96 2.0792E+13 0.256300 6 1.460441E-09 7.8133E-10 0.009402 0.97 1.5683E+13 0.176110 7 1.461664E-09 4.9410E-10 0.006139 0.99 1.3354E+13 0.129012 8 1.462403E-09 4.1545E-10 0.004442 0.99 1.3309E+13 0.098966 9 1.462902E-09 3.9759E-10 0.003436 0.99 1.3088E+13 0.078719 10 1.463179E-09 2.8602E-10 0.002765 1.00 3.8119E+13 0.065484 11 1.463384E-09 2.1282E-10 0.002242 0.98 1.1321E+13 0.053389 12 1.463541E-09 1.6792E-10 0.001848 1.00 1.0706E+13 0.045072 13 1.463656E-09 1.4174E-10 0.001574 1.00 1.0159E+13 0.037949 14 1.463726E-09 1.2204E-10 0.001387 1.00 3.8970E+13 0.034007 15 1.463797E-09 1.0593E-10 0.001191 1.00 8.4646E+12 0.028639 16 1.463858E-09 9.1945E-11 0.001032 1.00 9.3712E+12 0.024639 17 1.463912E-09 7.9208E-11 0.000905 1.00 4.4138E+12 0.021516 18 1.463962E-09 6.8581E-11 0.000797 1.00 3.8758E+12 0.018926 19 1.464008E-09 5.8751E-11 0.000704 1.00 3.3709E+12 0.016755 20 1.464047E-09 4.9898E-11 0.000624 1.00 2.7993E+12 0.014849 21 1.464083E-09 4.2806E-11 0.000556 1.00 2.6705E+12 0.013222 22 1.464115E-09 3.6782E-11 0.000497 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 3: 8: 4 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 3: 8: 4 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.174405E-09 4.2959E-09 1.000000 ---- 2.3250E+16 0.984359 2 1.420505E-09 3.1743E-09 0.396767 0.30 6.9221E+13 0.907372 3 1.444168E-09 2.0052E-09 0.089430 0.72 4.5968E+13 0.694055 4 1.454057E-09 1.3690E-09 0.034270 0.89 4.8674E+13 0.419448 5 1.458274E-09 8.6717E-10 0.016826 0.96 2.0788E+13 0.256267 6 1.460463E-09 7.8120E-10 0.009405 0.97 1.5684E+13 0.176140 7 1.461688E-09 4.9396E-10 0.006141 0.99 1.3352E+13 0.128998 8 1.462428E-09 4.1523E-10 0.004443 0.99 1.3214E+13 0.098975 9 1.462928E-09 3.9740E-10 0.003435 0.99 1.3117E+13 0.078717 10 1.463206E-09 2.8605E-10 0.002766 1.00 3.8102E+13 0.065479 11 1.463413E-09 2.1273E-10 0.002244 0.98 1.1321E+13 0.053379 12 1.463572E-09 1.6787E-10 0.001850 1.00 1.0662E+13 0.045062 13 1.463688E-09 1.4160E-10 0.001575 1.00 1.0133E+13 0.037932 14 1.463761E-09 1.2188E-10 0.001387 1.00 3.7226E+13 0.033964 15 1.463833E-09 1.0578E-10 0.001193 1.00 8.4729E+12 0.028630 16 1.463894E-09 9.1805E-11 0.001033 1.00 9.3683E+12 0.024639 17 1.463950E-09 7.9081E-11 0.000906 1.00 4.3915E+12 0.021519 18 1.464001E-09 6.8469E-11 0.000797 1.00 3.8771E+12 0.018917 19 1.464047E-09 5.8729E-11 0.000704 1.00 3.3651E+12 0.016751 20 1.464087E-09 5.0033E-11 0.000625 1.00 2.7997E+12 0.014839 21 1.464123E-09 4.2681E-11 0.000556 1.00 2.6708E+12 0.013219 22 1.464155E-09 3.6680E-11 0.000497 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 3:21:10 ----------------------------------------------------------- ============================================================== I Output from program Shells, I I a spherical-Earth, thin-shell program for computing time- I I averaged (non-elastic) neotectonic deformation of plates I I with realistic frictional/dislocation-creep rheology. I I Distinct thicknesses and thermal and mechanical I I properties are read for the crust and mantle layers I I of the lithosphere. I I Faults may be included, with specified dip and friction. I I Also, different elements *may* have different rheologies. I I (*This is the primary new feature in version 5.0+.) I I The velocity below the base of the model may be fixed, I I (to represent subduction and other convection), I I or shear traction on the base of the lithosphere may I I be set to zero, or basal shear traction may be adjusted I I (in a series of runs) to obtain desired plate velocities.I I by I I Peter Bird & Xianghong Kong I I Department of Earth, Planetary, and Space Sciences I I University of California I I Los Angeles, CA 90095-1567 I I Peter Bird's version 5.0* of 29 January 2018 I ================================================================ ----------------------------------------------------------- Run began on 2023. 2.27 at 3:21:10 ----------------------------------------------------------- =========================================================================================================================== Table of alternative Lithospheric Rheologies defined and used: LR# fFric cFric Biot Byerly aCreep(1) aCreep(2) bCreep(1) bCreep(2) cCreep(1) cCreep(2) dCreep(1) dCreep(2) eCreep 0 0.085 0.205 1.000 0.000 1.32E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 1 0.085 0.850 1.000 0.000 1.38E+06 2.23E+04 10524. 16494. 0.0079 0.0104 5.00E+08 5.00E+08 0.26667 =========================================================================================================================== Iteration history: Relative Corre- Maximum mean Relative lation vertically vertically Maximum mean with integrated integrated Iter- RMS velocity velocity last stress stress ation velocity change change change error error 1 1.323284E-09 3.3225E-09 1.000000 ---- 8.7655E+15 0.977484 2 1.395134E-09 2.2125E-09 0.524465 -0.20 9.2568E+13 0.852628 3 1.417396E-09 1.0018E-09 0.198095 0.86 4.6808E+13 0.711788 4 1.436134E-09 8.6057E-10 0.090850 0.95 3.3080E+13 0.459668 5 1.440736E-09 5.2968E-10 0.035295 0.97 2.4204E+13 0.274983 6 1.442338E-09 3.6520E-10 0.019889 0.97 2.3874E+13 0.180146 7 1.444109E-09 3.4285E-10 0.014597 0.99 1.9363E+13 0.130653 8 1.445822E-09 3.2511E-10 0.011103 1.00 1.5530E+13 0.100619 9 1.447359E-09 4.0619E-10 0.008609 1.00 1.1233E+13 0.080255 10 1.448748E-09 4.7747E-10 0.006873 1.00 1.0046E+13 0.065365 11 1.450052E-09 5.2967E-10 0.005655 1.00 8.6806E+12 0.053947 12 1.451256E-09 5.4796E-10 0.004716 1.00 8.3072E+12 0.044998 13 1.452348E-09 5.3027E-10 0.003945 1.00 7.2773E+12 0.037764 14 1.453329E-09 4.8179E-10 0.003340 1.00 5.6089E+12 0.032025 15 1.454265E-09 4.2065E-10 0.002876 1.00 5.5623E+12 0.027428 16 1.455151E-09 3.5734E-10 0.002514 1.00 6.5779E+12 0.023937 17 1.455933E-09 2.9911E-10 0.002204 1.00 6.3350E+12 0.020702 18 1.456641E-09 2.4727E-10 0.001941 1.00 8.4120E+12 0.018278 19 1.457298E-09 2.0292E-10 0.001757 1.00 3.1099E+13 0.016797 20 1.457917E-09 1.6644E-10 0.001527 1.00 4.6319E+12 0.014142 21 1.458491E-09 1.3481E-10 0.001363 1.00 2.5132E+12 0.012457 22 1.459011E-09 1.0937E-10 0.001219 1.00 2.3311E+12 0.011120 23 1.459474E-09 8.7269E-11 0.001095 1.00 2.1478E+12 0.009980 24 1.459891E-09 7.0249E-11 0.000980 1.00 1.9641E+12 0.009009 25 1.460267E-09 6.3083E-11 0.000880 1.00 1.7873E+12 0.008144 26 1.460606E-09 5.6572E-11 0.000791 1.00 1.7414E+12 0.007380 27 1.460914E-09 5.0620E-11 0.000713 1.00 1.6231E+12 0.006657 28 1.461194E-09 4.5223E-11 0.000643 1.00 1.4874E+12 0.006047 29 1.461450E-09 4.0171E-11 0.000582 1.00 1.4023E+12 0.005507 30 1.461682E-09 3.5761E-11 0.000527 1.00 1.4757E+12 0.005007 31 1.461892E-09 3.1721E-11 0.000477 1.00 CONVERGED !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ----------------------------------------------------------- Run ended on 2023. 2.27 at 3:39: 8 -----------------------------------------------------------