[cig-commits] r4695 - in mc/3D/CitcomS/trunk/examples: Cookbook1
Cookbook2 Cookbook3 Cookbook4 Cookbook5 Cookbook6
tan2 at geodynamics.org
tan2 at geodynamics.org
Wed Oct 4 11:34:45 PDT 2006
Author: tan2
Date: 2006-10-04 11:34:45 -0700 (Wed, 04 Oct 2006)
New Revision: 4695
Modified:
mc/3D/CitcomS/trunk/examples/Cookbook1/cookbook1.cfg
mc/3D/CitcomS/trunk/examples/Cookbook2/cookbook2.cfg
mc/3D/CitcomS/trunk/examples/Cookbook3/cookbook3.cfg
mc/3D/CitcomS/trunk/examples/Cookbook4/cookbook4.cfg
mc/3D/CitcomS/trunk/examples/Cookbook5/cookbook5.cfg
mc/3D/CitcomS/trunk/examples/Cookbook6/cookbook6.cfg
Log:
* Removed some unnecessary parameters
* Changed the layout slightly
* Updated the cookbook scripts to be consistent with the manual
Modified: mc/3D/CitcomS/trunk/examples/Cookbook1/cookbook1.cfg
===================================================================
--- mc/3D/CitcomS/trunk/examples/Cookbook1/cookbook1.cfg 2006-10-04 18:33:24 UTC (rev 4694)
+++ mc/3D/CitcomS/trunk/examples/Cookbook1/cookbook1.cfg 2006-10-04 18:34:45 UTC (rev 4695)
@@ -14,3 +14,5 @@
num_perturbations = 1
perturbl = 3
perturbm = 2
+perturblayer = 5
+perturbmag = 0.05
Modified: mc/3D/CitcomS/trunk/examples/Cookbook2/cookbook2.cfg
===================================================================
--- mc/3D/CitcomS/trunk/examples/Cookbook2/cookbook2.cfg 2006-10-04 18:33:24 UTC (rev 4694)
+++ mc/3D/CitcomS/trunk/examples/Cookbook2/cookbook2.cfg 2006-10-04 18:34:45 UTC (rev 4695)
@@ -7,7 +7,7 @@
monitoringFrequency = 30 ; how often to monitor the computation
[CitcomS.solver]
-datafile = cookbook2 ; location/name of output files
+datafile = cookbook2 ; prefix of output filenames
# Modify the layout of the mesh.
[CitcomS.solver.mesher]
@@ -29,3 +29,4 @@
[CitcomS.solver.ic]
num_perturbations = 1
perturbmag = 0.0
+
Modified: mc/3D/CitcomS/trunk/examples/Cookbook3/cookbook3.cfg
===================================================================
--- mc/3D/CitcomS/trunk/examples/Cookbook3/cookbook3.cfg 2006-10-04 18:33:24 UTC (rev 4694)
+++ mc/3D/CitcomS/trunk/examples/Cookbook3/cookbook3.cfg 2006-10-04 18:34:45 UTC (rev 4695)
@@ -7,7 +7,7 @@
monitoringFrequency = 25 ; how often to monitor the computation
[CitcomS.solver]
-datafile = cookbook3 ; location/name of output files
+datafile = cookbook3 ; prefix of output filenames
rayleigh = 1e6 ; Rayleigh number
# Modify the layout of the mesh.
@@ -30,3 +30,4 @@
visc_min = 1.0
VMAX = on
visc_max = 100.0
+
Modified: mc/3D/CitcomS/trunk/examples/Cookbook4/cookbook4.cfg
===================================================================
--- mc/3D/CitcomS/trunk/examples/Cookbook4/cookbook4.cfg 2006-10-04 18:33:24 UTC (rev 4694)
+++ mc/3D/CitcomS/trunk/examples/Cookbook4/cookbook4.cfg 2006-10-04 18:34:45 UTC (rev 4695)
@@ -1,41 +1,30 @@
# Cookbook 4: Regionally Refined Meshes
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS]
-
steps = 250 ; number of time steps
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.controller]
-
monitoringFrequency = 10 ; how often to monitor the computation
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver]
-
rayleigh = 1e6 ; Rayleigh number
-datafile = cookbook4 ; location/name of output files
+datafile = cookbook4 ; prefix of output filenames
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.ic] ; Initial Conditions
-
+[CitcomS.solver.ic]
num_perturbations = 1
perturbmag = 0.05
+perturblayer = 10
# Perturb the initial temperature gradient in the longitudinal
# direction.
perturbl = 1
perturbm = 0
-perturblayer = 10
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.mesher] ; Mesh
-
+[CitcomS.solver.mesher]
# Read uneven mesh point coordinates from 'coor.dat'.
coor = on
coor_file = coor.dat
@@ -46,17 +35,15 @@
nodex = 33
nodey = 17
nodez = 17
-levels = 1
theta_min = 1
theta_max = 2
fi_min = 0
fi_max = 1
+radius_inner = 0.55
+radius_outer = 1
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.visc] ; Viscosity
-
-# Assign the viscosities.
+[CitcomS.solver.visc]
VISC_UPDATE = on
num_mat = 4
visc0 = 1,1,1,1
@@ -67,3 +54,4 @@
visc_min = 1.0
VMAX = on
visc_max = 100.0
+
Modified: mc/3D/CitcomS/trunk/examples/Cookbook5/cookbook5.cfg
===================================================================
--- mc/3D/CitcomS/trunk/examples/Cookbook5/cookbook5.cfg 2006-10-04 18:33:24 UTC (rev 4694)
+++ mc/3D/CitcomS/trunk/examples/Cookbook5/cookbook5.cfg 2006-10-04 18:34:45 UTC (rev 4695)
@@ -1,56 +1,32 @@
# Cookbook 5: Subduction Models with Trench Rollback
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS]
-
steps = 1000 ; number of time steps
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.controller]
-
monitoringFrequency = 10 ; how often to monitor the computation
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver]
-
-datafile = cookbook5_output
datafile = cookbook5
datafile_old = ./restart_files/cookbook5
-
-# Rayleigh number
rayleigh = 4.07e+08
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.bc] ; Boundary Conditions
-
+[CitcomS.solver.bc]
topvbc = 1
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver.param]
-
file_vbcs = on
-
start_age = 55
# Since the starting age is set to 55 Ma, there will be 56 velocity
# files, one for each Ma (bvel.dat0, bvel.dat1, ... bvel.dat56).
vel_bound_file = ./velocity/bvel.dat
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.ic] ; Initial Conditions
-
+[CitcomS.solver.ic]
restart = on
solution_cycles_init = 0
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver.mesher]
-
coor = on
coor_file = ./coor.dat
nprocx = 2
@@ -65,23 +41,14 @@
fi_max = 0.5
radius_inner = 0.7
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.const] ; Dimensional Information
-
-refvisc = 4e+21
-
-
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.visc] ; Viscosity
-
+[CitcomS.solver.visc]
num_mat = 4
visc0 = 100,0.003,1,2
TDEPV = on
viscE = 24,24,24,24
viscT = 0.182,0.182,0.182,0.182
-viscZ = 357.6,357.6,357.6,357.6
VMIN = on
visc_min = 0.01
VMAX = on
visc_max = 100.0
+
Modified: mc/3D/CitcomS/trunk/examples/Cookbook6/cookbook6.cfg
===================================================================
--- mc/3D/CitcomS/trunk/examples/Cookbook6/cookbook6.cfg 2006-10-04 18:33:24 UTC (rev 4694)
+++ mc/3D/CitcomS/trunk/examples/Cookbook6/cookbook6.cfg 2006-10-04 18:34:45 UTC (rev 4695)
@@ -1,27 +1,19 @@
-# Cookbook 6
+# Cookbook 6: Pseudo-Free-Surface Formulation
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS]
-
steps = 100
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.controller]
-
monitoringFrequency = 10
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver]
-
datafile = cookbook6
rayleigh = 4.312616e+08
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.mesher] ; Mesh
-
+[CitcomS.solver.mesher]
nprocx = 4
nprocy = 4
nprocz = 2
@@ -32,15 +24,11 @@
coor_file = ./coord.dat
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver.tsolver]
-
fixed_timestep = 7.77e-10
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver.vsolver]
-
precond = on
accuracy = 1e-10
tole_compressibility = 1e-06
@@ -52,12 +40,7 @@
piterations = 100000
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.ic] ; Initial Conditions
-
-num_perturbations = 0
-half_space_age = 100.0
-
+[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.
@@ -66,16 +49,14 @@
#
# 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] ; Boundary Conditions
-
+[CitcomS.solver.bc]
topvbc = 2
pseudo_free_surf = on
bottbcval = 0.82
@@ -87,28 +68,23 @@
#pseudo_free_surf = off
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[CitcomS.solver.param]
-
start_age = 60
mantle_temp = 0.82
-# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-[CitcomS.solver.visc] ; Viscosity
-
+[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] ; Dimensional Information
-
+[CitcomS.solver.const]
layerd = 6.371e+06
density = 3270.0
thermdiff = 1.0e-06
gravacc = 10.0
thermexp = 3.0e-05
refvisc = 1.0e+21
+
More information about the cig-commits
mailing list