Revision 4cab4ca49c9e5ed08de64533ca79429e2f41491a authored by Eh Tan on 12 November 2007, 22:49:46 UTC, committed by Eh Tan on 12 November 2007, 22:49:46 UTC
........ r8275 | becker | 2007-11-11 12:37:22 -0800 (Sun, 11 Nov 2007) | 3 lines Synching with the new release, none of those changes should be relevant. ........
1 parent 0461d2e
cookbook6.cfg
# Cookbook 6: Pseudo-Free-Surface Formulation
[CitcomS]
steps = 100
[CitcomS.controller]
monitoringFrequency = 10
[CitcomS.solver]
datafile = cookbook6
rayleigh = 4.312616e+08
[CitcomS.solver.mesher]
nprocx = 2
nprocy = 2
nprocz = 1
nodex = 61
nodey = 61
nodez = 25
coor = 1
coor_file = ./coord.dat
[CitcomS.solver.tsolver]
fixed_timestep = 7.77e-10
[CitcomS.solver.vsolver]
precond = on
accuracy = 1e-10
tole_compressibility = 1e-06
vlowstep = 100000
piterations = 100000
[CitcomS.solver.ic]
# To try a different hot blob, one can specify three parameters:
# 1. blob_center: Coordinates(theta,fi,r) of the center.
# 2. blob_radius: Radius as non-dimensional distance.
# 3. blob_dT: The amount of temperature increase/decrease
# within the blob as non-dimensional temperature.
#
# Note: It is required to keep "tic_method = 2" to have the initial
# hot blob.
tic_method = 2
blob_center = 1.570800e+00,1.570800e+00,9.246600e-01
blob_radius = 6.278334e-02
blob_dT = 0.18
half_space_age = 100.0
[CitcomS.solver.bc]
topvbc = 2
pseudo_free_surf = on
bottbcval = 0.82
# To convert pseudo-free-surface B.C. to free-slip B.C., do the
# following:
#topvbc = 0
#pseudo_free_surf = off
[CitcomS.solver.param]
start_age = 60
mantle_temp = 0.82
[CitcomS.solver.visc]
TDEPV = on
visc0 = 1,1,1,1
viscE = 9.50614,9.50614,9.50614,9.50614
viscT = 1.02126,1.02126,1.02126,1.02126
[CitcomS.solver.const]
radius = 6.371e+06
density = 3270.0
thermdiff = 1.0e-06
gravacc = 10.0
thermexp = 3.0e-05
refvisc = 1.0e+21
Computing file changes ...