Skip to content

Commit

Permalink
address comments and update test
Browse files Browse the repository at this point in the history
  • Loading branch information
@jdannberg
jdannberg committed Jun 12, 2024
1 parent 9814c9f commit d691b1f
Show file tree
Hide file tree
Showing 4 changed files with 36 additions and 119 deletions.
2 changes: 1 addition & 1 deletion ...antle_convection_with_continents_in_annulus/doc/mantle_convection_in_annulus.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -83,7 +83,7 @@ Results of the reference model (R) after 500 Ma of mantle convection simulation.
temperature (b) within the modelled domain, showing five slabs actively subducting below the continents (pink).
(c) Solid lines indicate average surface velocity (blue)
and mobility M (red) of model R as well as their dashed time-averages Vsurf=2.1 cm/a and M = 2.218.
Mobility is defined as the ratio of rms surface velocity to rms velocity averaged over the entire 3D domain{cite}`tackley:2000`.
Mobility is defined as the ratio of rms surface velocity to rms velocity averaged over the entire 3D domain {cite}`tackley:2000`.
(d) Radial velocity of all tracers defined as
slabs plotted at their position in the model where blue indicates sinking slabs and red slowly rising. Shown tracers are 300 K colder than the radial averaged
temperature at similar depth and automatically obtained (see methods section of the publication).
Expand Down
3 changes: 2 additions & 1 deletion cookbooks/mantle_convection_with_continents_in_annulus/modelR.prm
View file
Original file line number Diff line number Diff line change
Expand Up @@ -15,6 +15,7 @@ subsection Solver parameters
subsection Stokes solver parameters
set Stokes solver type = block GMG
set Number of cheap Stokes solver steps = 5000
set Maximum number of expensive Stokes solver steps = 0
set Linear solver tolerance = 1e-6
end
end
Expand Down Expand Up @@ -204,7 +205,7 @@ subsection Postprocess

subsection Generator
subsection Reference cell
set Number of particles per cell per direction = 3
set Number of particles per cell per direction = 4
end
end

Expand Down
103 changes: 13 additions & 90 deletions tests/cookbook_mantle_convection_annulus/screen-output
View file
Original file line number Diff line number Diff line change
@@ -1,68 +1,24 @@
-----------------------------------------------------------------------------
-- This is ASPECT --
-- The Advanced Solver for Planetary Evolution, Convection, and Tectonics. --
-----------------------------------------------------------------------------
-- . version 2.6.0-pre (add_mantle_conv_annulus_cookbook, 31f9f36ad)
-- . using deal.II 9.5.1
-- . with 32 bit indices and vectorization level 1 (128 bits)
-- . using Trilinos 13.2.0
-- . using p4est 2.3.2
-- . using Geodynamic World Builder 0.6.0
-- . running in DEBUG mode
-- . running with 1 MPI process
-----------------------------------------------------------------------------

-----------------------------------------------------------------------------
-- For information on how to cite ASPECT, see:
-- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&GWB=1&particles=1&sha=31f9f36ad&src=code
-----------------------------------------------------------------------------
Number of active cells: 48 (on 2 levels)
Number of degrees of freedom: 1,272 (480+72+240+240+240)

*** Timestep 0: t=0 years, dt=0 years
Solving temperature system... 0 iterations.
Advecting particles... done.
Solving 1_mantle system ... 0 iterations.
Solving 3_continent system ... 0 iterations.
Rebuilding Stokes preconditioner...
Solving Stokes system... 0+17 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 1
Solving Stokes system... 14+0 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 1: 0.0570421

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+17 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.21808
Solving Stokes system... 14+0 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 2: 0.0124028

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+16 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 0.024248
Solving Stokes system... 12+0 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 3: 0.00146411

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+15 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 0.00816953
Solving Stokes system... 11+0 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 4: 0.00017388

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+14 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 5: 0.0059109

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+14 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 6: 0.00430698

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+14 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 7: 0.00320187

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+14 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 8: 0.002389

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+13 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 9: 0.00179171

Rebuilding Stokes preconditioner...
Solving Stokes system... 0+13 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 10: 0.00135323
Solving Stokes system... 8+0 iterations.
Relative nonlinear residual (Stokes system) after nonlinear iteration 5: 2.75785e-05


Postprocessing:
Expand All @@ -71,51 +27,18 @@ Number of degrees of freedom: 1,272 (480+72+240+240+240)
that does not have any quadrature points in it.
Consider reducing the number of depth layers for averaging.

RMS, max velocity: 0.0131 m/year, 0.0283 m/year
RMS, max velocity: 0.0139 m/year, 0.023 m/year
Temperature min/avg/max: 300 K, 2067 K, 3700 K
Writing graphical output: output-cookbook_mantle_convection_annulus/solution/solution-00000
Heat fluxes through boundary parts: 9.327e+06 W, 1.851e+07 W
Heat fluxes through boundary parts: 1.168e+07 W, 2.059e+07 W
Writing depth average: output-cookbook_mantle_convection_annulus/depth_average
Writing particle output: output-cookbook_mantle_convection_annulus/particles/particles-00000
Writing stokes residuals output-cookbook_mantle_convection_annulus/stokes_residuals.txt
Mobility: 0.399
Mobility: 3.04e-06

Skipping mesh refinement, because the mesh did not change.

Termination requested by criterion: end time


+----------------------------------------------+------------+------------+
| Total wallclock time elapsed since start | 3.35s | |
| | | |
| Section | no. calls | wall time | % of total |
+----------------------------------+-----------+------------+------------+
| Assemble Stokes system | 10 | 0.869s | 26% |
| Assemble composition system | 2 | 0.153s | 4.6% |
| Assemble temperature system | 1 | 0.0464s | 1.4% |
| Build Stokes preconditioner | 10 | 0.529s | 16% |
| Build composition preconditioner | 2 | 0.000633s | 0% |
| Build temperature preconditioner | 1 | 0.000487s | 0% |
| Initialization | 1 | 0.628s | 19% |
| Particles: Advect | 2 | 0.00416s | 0.12% |
| Particles: Copy | 1 | 0.000155s | 0% |
| Particles: Generate | 1 | 0.0334s | 1% |
| Particles: Initialization | 1 | 0.000207s | 0% |
| Particles: Initialize properties | 1 | 0.00762s | 0.23% |
| Particles: Sort | 2 | 0.0321s | 0.96% |
| Particles: Update properties | 1 | 0.00139s | 0% |
| Postprocessing | 1 | 0.355s | 11% |
| Refine mesh structure, part 1 | 1 | 0.11s | 3.3% |
| Setup dof systems | 1 | 0.0237s | 0.71% |
| Setup initial conditions | 1 | 0.135s | 4% |
| Setup matrices | 1 | 0.0523s | 1.6% |
| Solve Stokes system | 10 | 0.276s | 8.2% |
| Solve composition system | 2 | 0.00116s | 0% |
| Solve temperature system | 1 | 0.000882s | 0% |
+----------------------------------+-----------+------------+------------+

-- Total wallclock time elapsed including restarts: 3s
-----------------------------------------------------------------------------
-- For information on how to cite ASPECT, see:
-- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&GWB=1&particles=1&sha=31f9f36ad&src=code
-----------------------------------------------------------------------------
47 changes: 20 additions & 27 deletions tests/cookbook_mantle_convection_annulus/statistics
View file
Original file line number Diff line number Diff line change
Expand Up @@ -7,30 +7,23 @@
# 7: Number of degrees of freedom for all compositions
# 8: Nonlinear iteration number
# 9: Iterations for temperature solver
# 10: Iterations for composition solver 1
# 11: Iterations for composition solver 2
# 12: Iterations for Stokes solver
# 13: Velocity iterations in Stokes preconditioner
# 14: Schur complement iterations in Stokes preconditioner
# 15: RMS velocity (m/year)
# 16: Max. velocity (m/year)
# 17: Minimal temperature (K)
# 18: Average temperature (K)
# 19: Maximal temperature (K)
# 20: Average nondimensional temperature (K)
# 21: Visualization file name
# 22: Outward heat flux through boundary with indicator 0 ("bottom") (W)
# 23: Outward heat flux through boundary with indicator 1 ("top") (W)
# 24: Number of advected particles
# 25: Particle file name
# 26: Mobility
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 0 0 0 0 16 154 84 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 1 0 0 0 16 156 90 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 2 0 0 0 15 158 85 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 3 0 0 0 14 127 80 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 4 0 0 0 13 153 75 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 5 0 0 0 13 140 75 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 6 0 0 0 13 167 75 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 7 0 0 0 13 142 75 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 8 0 0 0 12 134 70 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 9 0 0 0 12 123 70 1.30678462e-02 2.82680195e-02 3.00000000e+02 2.06745340e+03 3.70000000e+03 5.19839236e-01 output-cookbook_mantle_convection_annulus/solution/solution-00000 9.32654848e+06 1.85052632e+07 48 output-cookbook_mantle_convection_annulus/particles/particles-00000 3.98946605e-01
# 10: Iterations for Stokes solver
# 11: Velocity iterations in Stokes preconditioner
# 12: Schur complement iterations in Stokes preconditioner
# 13: RMS velocity (m/year)
# 14: Max. velocity (m/year)
# 15: Minimal temperature (K)
# 16: Average temperature (K)
# 17: Maximal temperature (K)
# 18: Average nondimensional temperature (K)
# 19: Visualization file name
# 20: Outward heat flux through boundary with indicator 0 ("bottom") (W)
# 21: Outward heat flux through boundary with indicator 1 ("top") (W)
# 22: Number of advected particles
# 23: Particle file name
# 24: Mobility
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 0 0 13 15 15 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 1 0 13 15 15 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 2 0 11 13 13 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 3 0 10 12 12 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00 "" 0.00000000e+00 0.00000000e+00 0 "" 0.00000000e+00
0 0.000000000000e+00 0.000000000000e+00 48 552 240 480 4 0 7 9 9 1.39049456e-02 2.29754397e-02 3.00000000e+02 2.06745340e+03 3.70000000e+03 5.19839236e-01 output-cookbook_mantle_convection_annulus/solution/solution-00000 1.16815394e+07 2.05949390e+07 3840 output-cookbook_mantle_convection_annulus/particles/particles-00000 3.03648058e-06

0 comments on commit d691b1f

Please sign in to comment.