# Cookbook 5: Subduction Models with Trench Rollback
[CitcomS]
steps = 1000 ; number of time steps
[CitcomS.controller]
monitoringFrequency = 10 ; how often outputs are created
[CitcomS.solver]
datafile = cookbook5
datadir_old = ./restart_file
datafile_old = cookbook5
rayleigh = 4.07e+08
[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]
restart = on
solution_cycles_init = 0
[CitcomS.solver.mesher]
coor = on
coor_file = ./coor.dat
nprocx = 2
nprocy = 8
nprocz = 4
nodex = 17
nodey = 65
nodez = 33
theta_min = 1.47
theta_max = 1.67
fi_min = 0
fi_max = 0.5
radius_inner = 0.7
[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
VMIN = on
visc_min = 0.01
VMAX = on
visc_max = 100.0