Repository navigation
Expand file tree
/
Copy pathdata-example
More file actions
82 lines (73 loc) · 3.47 KB
/
Copy pathdata-example
File metadata and controls
82 lines (73 loc) · 3.47 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
# -*- getpot -*- (GetPot mode activation for emacs)
#-------------------------------------------------
# Data file for the Bidomain solver
#-------------------------------------------------
# FasterCard V.10 Data file
[schaf]
[./physics]
tau_in = 0.3
tau_out = 6.0
tau_open = 120.0
tau_close = 150.0
tau_close_hetero = 0 # tells if heterog. close is considered (on cases where exists)
tau_close_endo = 140.0
tau_close_mcell = 105.0
tau_close_epi = 105.0
tau_close_RV = 120.0
tau_in_hithm = 0.3
tau_out_hithm = 6.0
tau_open_hithm = 120.0
tau_close_hithm = 150.0
v_min = 0.0
v_max = 1.0
vcrit = 0.13
gateOutput = 0 #tells if using mMS ionic model instead of MS
solver = 3 #1: explicit; 2: semi-implicit on linear part; 3: Newton
ionic_params_from_file = 0 #tells if reading ionic parameters from a file
ionic_params_file_name = mMS.ipar
[bidomain]
[./physics]
isotropic = 0 #tells to use isotropic models; transverse conductivities set euql to longitudinal; fake fibres are generated
conductivity_i_l = 3.0e-3
conductivity_i_t = 3.0e-4
conductivity_i_l_hitm = 3.0e-3 # conductivity on hithsmus region, longitudinal
conductivity_i_t_hitm = 3.0e-4 # conductivity on hithsmus region, transversal
lambda = 1.0 # ratio between intra and extra conductivity tensors
read_cond_values = 1 #tells to read conductivities from a file
cond_file_name = conductivities.cpar #conductivity file
delta = 1.0 # correction term for conduction velocity
Cm = 1.0 # membrane capacitance mF/cm^2= 1 micro F/cm^2
test_case = 1 # test case to run
CVadaption = 1 # adapts the Conduction velocity to the underpinning electro physiology
DiffCurr = 0 # switches diffusive current on (experimental)
solveEiko = 1 #toggles eikonal solver on/off. To use with DiffCurr on
tactonly = 0 # 0: run simulation; 1: eval activation times only for the first stimulus
Time_period = 800.0
t_stim = 0.6
t_s1s2 = 256.0
end_time = 700.0
nbprepacing = 1
left_radius = 1.5
right_radius = 1.5
IappMod = 1.0
stimnodes_file =./stim.vtx #file of the nodes where a stimulus is applied (working with case 999)
[discretization]
mesh_dir = # mesh directory
mesh_file = # mesh file name (.mesh)
fibers_dir = # fibers directory
fibers_file = # fibres file name
dist_dir = # distances directory
dist_file = # distance directory file name (for close heterogeneous)
timestep = 0.1 # time step for activation time time stepping
nsub_ode = 1
eikotype = Dijkstra #types: Dijkstra, FastMarching
# timestep_ode deprecated; use nsub_ode to define the nb of ode subiterations
#timestep_ode = 0.1 # time step for the solution of the ionic model
fastSolver = 0 # 0: solve ionic model; 1: fast evaluation (vm=0,0.5,1; w from analytic solution)
[output]
output = 0
binary = 1
writer = vtk #vtk, ensight
verbose = 1
dir = ./output
case = heart