# specify input file, input.17.100.55  , for model cbl 
# Internal heat production: 100 ,  Rayleigh number: 10000000 ,Buoyance Ratio:0.55,  Resolution: 33 

# Parameters that Control Input Files

reference_state=1
refstate_file="refstate.dat"

mineral_physics_model=3

file_vbcs=off
vel_bound_file="bvel.dat"

coor=0
coor_file=""
coor_refine=0.1,0.15,0.1,0.2

lith_age=off
lith_age_file="age.dat"
lith_age_time=off

# Parameters that Control Output Files

output_format=ascii

output_optional="surf,botm,horiz_avg,comp_nd"

datadir="../output/cbl.17.100.55/"

datafile="cbl.17.100.55"

storage_spacing=200

checkpointFrequency=1000

output_ll_max=20

self_gravitation=off

use_cbf_topo=off

# Mesh and Processors Setup

nproc_surf=12

nprocx=1
nprocy=1
nprocz=4

nodex=33
nodey=33
nodez=33

mgunitx=32
mgunity=32
mgunitz=32
levels=1

#nodex=1+nprocx*mgunitx*2^(levels-1)

# Domain Size

theta_min=1.0708
theta_max=2.0708
fi_min=0
fi_max=1

radius_outer=1
radius_inner=0.54623

# Restarting the Code

restart=off
post_p=off

datadir_old=""
datafile_old=""
solution_cycles_init=0

# Run Length

minstep=1
maxstep=2001
maxtotstep=2001

cpu_limits_in_seconds=360000000

# Initial Conditions

tic_method=0
num_perturbations=1
perturbmag=0.05
perturbl=1
perturbm=1
perturblayer=5

# when tic_method=-1 ***
# datadir_old=""
# datafile_old=""
# solution_cycles_init=0
# zero_elapsed_time=on
# ***

# when tic_method=2  ***
# half_space_age=40
# mantle_temp=1.0
# blob_center=[-999,-999,-999]
# blob_radius=0.063
# blob_dT=0.18
# ***

# Boundary Conditions

topvbc=0
topvbxval=0
topvbyval=0

botvbc=0
botvbxval=0
botvbyval=0

side_sbcs=off
pseudo_free_surf=off

toptbc=1
toptbcval=0

bottbc=1
bottbcval=1

temperature_bound_adj=off
depth_bound_adj=0.157
width_bound_adj=0.08727
lith_age_depth=0.0314

# Non-dimensional Numbers

rayleigh=10000000

dissipation_number=0

gruneisen=0

surfaceT=0.1

Q0=100

# Depth Information

z_lith=0.014
z_410=0.06435
z_lmantle=0.105
z_cmb=0.45377

# Viscosity

Viscosity=system

visc_smooth_method=3

VISC_UPDATE=off

num_mat=4

visc0=1,1,1,1

TDEPV=off

rheol=3

viscE=1,1,1,1
viscT=1,1,1,1
viscZ=1,1,1,1

SDEPV=off

sdepv_expt=1,1,1,1
sdepv_misfit=0.001

PDEPV=off

pdepv_a=1e20,1e20,1e20,1d20
pdepv_b=0,0,0,0
pdepv_y=1e20,1e20,1e20,1d20

pdepv_eff=on

pdepv_offset=0

CDEPV=off

cdepv_ff=1,1

low_visc_channel=off
low_visc_wedge=off

lv_min_radius=0.9764
lv_max_radius=0.9921
lv_channel_thickness=0.0047
lv_reduction=0.5

VMIN=on
visc_min=0.001

VMAX=on
visc_max=1000

#Phase Change Information

Ra_410=0.0
clapeyron410=0.0235
transT410=0.78
width410=0.0058

Ra_670=0.0
clapeyron670=-0.0235
transT670=0.875
width670=0.0058

Ra_cmb=0.0
clapeyroncmb=-0.0235
transTcmb=0.875
widthcmb=0.0058

# Momentum Equation Solver Parameters

stokes_flow_only=off

remove_rigid_rotation=on

Solver=cgrad

node_assemble=on

mg_cycle=1
down_heavy=3
up_heavy=3
vlowstep=1000
vhighstep=3

piterations=1000

accuracy=1.0e-3
#tole_compressibility=1.0e-7

uzawa=cg
compress_iter_maxstep=100
#relative_err_accuracy=0.001

precond=on

aug_lagr=on
aug_number=2.0e3

# Energy Equation Solver Parameters

ADV=on

fixed_timestep=0.0

finetunedt=0.9

adv_sub_iterations=2
adv_gamma=0.5

filter_temp=off

monitor_max_T=on

inputdiffusivity=1

# Age Information

start_age=40.0
reset_startage=off

# Debugging Information

DESCRIBE=off
BEGINNER=off
VERBOSE=off

verbose=off

see_convergence=on

# HDF5 Output Parameters

#cb_block_size=1048576
#cb_buffer_size=4194304

#sieve_buf_size=1048576

#output_alignment=262144
#output_alignment_threshold=524288

#cache_mdc_nelmts=10330
#cache_rdcc_nelmts=521
#cache_rdcc_nbytes=1048576

# Tracer Parameters

tracer=on

tracer_ic_method = 0
tracers_per_element = 20
tracer_file = tracer.dat

tracer_flavors = 2

ic_method_for_flavors = 0
z_interface = 0.58547

chemical_buoyancy = on

buoy_type = 1

buoyancy_ratio =0.55

regular_grid_deltheta = 1.0
regular_grid_delphi = 1.0

itracer_warnings=on

tracer_enriched=off
Q0_enriched=0.0

# Dimensional Information

radius=6371e3
density=3340.0
thermdiff=1.0e-6
gravacc=9.81
thermexp=3.0e-5
refvisc=1.0e+21
cp=1200
density_above=1030.0
density_below=6600.0

# Required Information

Problem=convection
Geometry=sphere
Spacing=regular