mars topography

cookbooks/mars_mola_topography/doc/mars_topo.md

@@ -0,0 +1,23 @@
+# Mars topography
+
+*This section was contributed by Cedric Thieulot and Bart Root.*
+
+From the original MOLA data (<https://attic.gsfc.nasa.gov/mola/>), we have created a topography file
+(present in the folder of this cookbook) with a 1 degree resolution in both longitude and latitude
+that is suitable for the Initial topography model.
+
+Note that for the time being the Initial topography model is not compatible with the sphere geometry
+model so the model geometry for this cookbook is a chunk with inner radius 2000 km and outer 
+radius 3390 km, the latter being the actual radius of the planet.
+
+The topography is then projected onto the surface of the planet, and as a first step both viscosity and
+density are kept constant in the domain. No-slip boundary conditions are prescribed on the bottom, and
+free slip is prescribed on all four sides while the top boundary is free.
+The model is isothermal, and there are no other buoyancy forces.
+
+```{figure-md} fig:mars-topo
+<img src="topo3.*" width="90%" />
+
+Mars topography.
+For visualization purposes the topography has been scaled up by a factor 10.
+```

cookbooks/mars_mola_topography/mars.prm

@@ -0,0 +1,116 @@
+# An input script to test the chunk geometry model
+# with additional initial topography.
+
+set Dimension = 3
+set CFL number                             = 1.0
+set End time                               = 0
+set Start time                             = 0
+set Adiabatic surface temperature          = 1613.0
+#set Surface pressure                       = 0
+set Use years in output instead of seconds = true
+set Nonlinear solver scheme                = no Advection, single Stokes
+set Output directory                       = output-mola
+set Pressure normalization                 = no
+
+subsection Gravity model
+  set Model name = radial constant
+  subsection Radial constant
+    set Magnitude = 3.71
+  end
+end
+
+subsection Solver parameters
+  subsection Stokes solver parameters
+    set Number of cheap Stokes solver steps = 0
+    set Maximum number of expensive Stokes solver steps = 3000
+  end
+end
+
+subsection Geometry model
+
+  set Model name = chunk
+  subsection Chunk
+    set Chunk minimum longitude = -150
+    set Chunk maximum longitude = -30
+    set Longitude repetitions = 40
+    set Chunk minimum latitude = -60
+    set Chunk maximum latitude = 60
+    set Latitude repetitions = 40
+    set Chunk inner radius = 2000e3
+    set Chunk outer radius = 3390e3
+    set Radius repetitions = 1
+  end
+
+  #set Model name = spherical shell
+  #subsection Spherical shell
+  #  set Inner radius = 2000e3
+  #  set Outer radius = 3390e3
+  #  set Cells along circumference = 12
+  #end
+
+  subsection Initial topography model
+    set Model name = ascii data
+    subsection Ascii data model
+      set Data directory       = /home/cedrict/ASPECT/pr6/aspect/build/
+      set Data file name       = MOLA_aspect.ascii
+      set Scale factor = 10
+    end
+  end
+end
+
+subsection Mesh refinement
+  set Initial adaptive refinement        = 0
+  set Initial global refinement          = 1
+end
+
+subsection Initial temperature model
+  set Model name = function
+  subsection Function
+    set Function expression = 1613.0
+  end
+end
+
+subsection Boundary temperature model
+  set List of model names = function
+  subsection Function
+    set Function expression = 1613.0
+  end
+end
+
+subsection Material model
+  set Model name = simple
+  subsection Simple model
+    set Reference density = 3340
+    set Reference specific heat = 1200
+    set Thermal expansion coefficient = 3e-5
+    set Viscosity = 1e22
+  end
+end
+
+subsection Boundary temperature model
+  set Fixed temperature boundary indicators = inner, outer
+end
+
+subsection Boundary velocity model
+  set Tangential velocity boundary indicators = east, west, north, south
+end
+
+subsection Boundary velocity model
+  set Zero velocity boundary indicators = inner
+end
+
+subsection Mesh deformation
+  set Mesh deformation boundary indicators = outer: free surface
+  subsection Free surface
+    set Surface velocity projection = normal
+  end
+end
+
+subsection Postprocess
+  set List of postprocessors = visualization, velocity statistics, basic statistics,  temperature statistics, topography
+  subsection Visualization
+    set Output format = vtu
+    set List of output variables = density, viscosity, strain rate, depth
+    set Time between graphical output = 0
+  end
+end

doc/sphinx/user/cookbooks/geophysical-setups.md

@@ -101,6 +101,7 @@ cookbooks/steinberger/doc/steinberger.md
 cookbooks/morency_doin_2004/doc/morency_doin_2004.md
 cookbooks/crustal_deformation/doc/crustal_deformation.md
 cookbooks/continental_extension/doc/continental_extension.md
+cookbooks/mars_mola_topography/doc/mars_topo.md
 cookbooks/inner_core_convection/doc/inner_core_convection.md
 cookbooks/global_melt/doc/global_melt.md
 cookbooks/mid_ocean_ridge/doc/mid_ocean_ridge.md