I was wondering if there is a way to use epicardial potentials calculated with a bidomain model as a kind of dirichlet boundary condition to compute potentials at the torso surface? I have calculated the potentials with the torso as bath surrounding, however I am willing to compute it in a uncoupled way, which means that there is a no flux boundary condition on the extracellular heart part as well instead of the perfect transmission between torso and heart.
I have tried to do it with "stimulus.stimtype = 10", however when using this stimtype it seems that I cannot define the specific nodes to prescribe the potential on by using the vtx_file.
I also tried "stimulus.stimtype =2" and using "stimulus.vtx_file", however when using this in combination with the La Place solver it had to be calculated with a vertex adjustment file per time step. Thereby, I thought of using the stimtype 2 with a predefined pulse file (trc), however the pulse shape is different for each node on the epicardial surface.
Is there another way to compute the uncoupled problem, either by assigning the stimulus differently than I described above, or by defining a no-flux boundary on heart-torso surface instead of the automatic assigned perfect transmission?