Mesh regions, which are taged other than 1, are not present in vm.igb results file
Dear OpenCARP Team,
I am trying to simulate simple tissue mesh with two different regions (1 and 2) tagged in block.elem file. But every time I use any elements tagged with region 2, resulting vm.igb file doesn't contain elements from region other than 1.
e.g. if I tag half of mesh elements to region 2, then resulting vm.igb file is proportionally smaller. I did quite a bit digging into documentation and code without any luck.
Questions: how to get vm of all tagged mesh regions in results file? is such behavior expected? should I use dynamic re-tagging?
Here is my config of regions (for testing purposes, parameters for each region are identical):
# Set monodomain conductivities
cmd += ['-num_gregions', 2,
'-gregion[0].name', "myocardium",
'-gregion[0].num_IDs', 1, # one tag will be given in the next line
'-gregion[0].ID[0]', 1, # use these settings for the mesh region with tag 1
# '-gregion[0].ID[1]', 2, # use these settings for the mesh region with tag 1
# mondomain conductivites will be calculated as half of the harmonic mean of intracellular
# and extracellular conductivities
'-gregion[0].g_el', 0.625, # extracellular conductivity in longitudinal direction
'-gregion[0].g_et', 0.236, # extracellular conductivity in transverse direction
'-gregion[0].g_en', 0.236, # extracellular conductivity in sheet direction
'-gregion[0].g_il', 0.174, # intracellular conductivity in longitudinal direction
'-gregion[0].g_it', 0.019, # intracellular conductivity in transverse direction
'-gregion[0].g_in', 0.019, # intracellular conductivity in sheet direction
'-gregion[0].g_mult', 0.5, # scale all conducitivites by a factor (to reduce conduction velocity)
'-gregion[1].name', "fibro",
'-gregion[1].num_IDs', 1, # one tag will be given in the next line
'-gregion[1].ID[0]', 2, # use these settings for the mesh region with tag 1
# mondomain conductivites will be calculated as half of the harmonic mean of intracellular
# and extracellular conductivities
'-gregion[1].g_el', 0.625, # extracellular conductivity in longitudinal direction
'-gregion[1].g_et', 0.236, # extracellular conductivity in transverse direction
'-gregion[1].g_en', 0.236, # extracellular conductivity in sheet direction
'-gregion[1].g_il', 0.174, # intracellular conductivity in longitudinal direction
'-gregion[1].g_it', 0.019, # intracellular conductivity in transverse direction
'-gregion[1].g_in', 0.019, # intracellular conductivity in sheet direction
'-gregion[1].g_mult', 0.5 # scale all conducitivites by a factor (to reduce conduction velocity)
]
# Define the ionic model to use
cmd += ['num_imp_regions', 2,
'imp_region[0].im', 'Courtemanche'
'imp_region[0].num_IDs', 1,
'imp_region[0].ID[0]', 1, # use this setting for the mesh region with tag 1
# 'imp_region[0].ID[1]', 2
'imp_region[1].im', 'Courtemanche'
'imp_region[1].num_IDs', 1,
'imp_region[1].ID[0]', 1, # use this setting for the mesh region with tag 1
# 'imp_region[1].ID[1]', 2
]