# README for http://geo.mff.cuni.cz/~sramek/maruska/ #

## Cases from Choblet et al. 2007 GJI doi:10.1111/j.1365-246X.2007.03419.x ##


### Case 3 ###

nondim r_cmb=0.5
Rayleigh based on thickness of mantle = 1.4e6 <=>
<=> Rayleigh based on planet radius = 1.12e7
isoviscous
no internal heating
no-slip

Output data:

* choblet3.radial_r_T.70000
  choblet3.radial_r_Tmin.70000
  choblet3.radial_r_Tmax.70000
  choblet3.radial_r_vh.70000
  choblet3.radial_r_vr.70000
  ... average, minimum and maximum temperature and rms of horizontal and 
      vertical components of velocity as functions of radius at time step 70000

* choblet3.time_qtop ... surface heat flux as function of time
* choblet3.time_qbot ... bottom heat flux as function of time
* choblet3.time_vrms ... vrms as function of time

Few simple plots:

* choblet3.time_vrms.png
* choblet3.time_qs.png
* choblet3.radial_r_T.70000.png
* choblet3.radial_r_vrms.70000.png

(Case 3 on 2015-05-26)


### Case 4 ###

nondim r_cmb=0.5
surface Rayleigh based on thickness of mantle = 1.e3 <=>
<=> surface Rayleigh based on planet radius = 8.e3
Frank-Kamenetskii viscosity: eta = exp(-15.*T)
internal heating based on layer thickness = 10 <=>
<=> internal heating based on planet radius = 40
no-slip

Output quantities as in Case 3, prefix = choblet4, additionally:
* choblet3.radial_r_eta.60000
* choblet3.radial_r_eta.60000.png

(Case 3 on 2015-05-31)


### Case 3 Free Slip ###

As Case 3 but free-slip boundary conditions
File prefix = choblet3FS

(Case 3FS on 2015-05-31)


### Case 4 Free Slip ###

As Case 4 but free-slip boundary conditions
File prefix = choblet4FS

(Case 4FS on 2015-05-31)


### Notes ###

Mesh 12 x 64^3 elements, equidistant in r,theta,phi

Quantities are non-dimensional.

Bottom heat flux has to be multiplied by r_cmb**2 (0.5**2=1/4) to compare it 
to surface heat flux.