Type debug to Turn On Debugging: debug Specify Debugging Destination, screen or file: screen List Functions to Debug - or all - end ends list Names following all will not be debugged - end ends : new_bdry : end Enter NX: 60 Enter NY: 60 Enter NZ: 60 Enter dx: 0.2 Enter dy: 0.2 Enter dz: 0.2 Enter dt: 0.01 Enter final time: 640 Enter the frequency of states dump for restart: 10000 Enter the frequency of plotting: 1000 Enter intracellular conductivity, (si_x,si_y,si_z): 1 0.081 0.081 Enter extracellular conductivity, (se_x,se_y,se_z): 1 12.4 12.4 Enter the fiber direction of the first layer in degree: 60 Enter the rate of the angle change between layers: -1.666667 Two ionic models are available FN (f) FHN-OMER (o) 3V-SIM (3) Enter a model (f(dflt),3): o Enter the initial value for v: 0 Three problems are available Monodomain (m) 1d Bidomain (1) 2d or 3d Bidomain (b) Enter a problem (m,1,b(dflt)): b Enter the resistance: 1 Enter the surface volume ratio: 1 Four solution methods are available. Crank_Nicolson (c) Semi_Implicit (i) Split_Dense (s) Split_ADI (t) Enter a method (c,i,s or t): t Enter y if you want to turn off bidomain: n Choose an initial condition. Plane Wave starting from left (p) Radial wave starting from a corner (c) Radial wave starting from the center (s) Read from a file (f) Enter the choice (p,c,s(default), or f): f Enter the file name from which you want to read: out/omer2d_plot/states.ts0004000 : n When the dimension of the init_file is 2, and the simulation dimension is 3, it needs to know how to fill the 3d volume with the initial wave. Stack along z axis (z) Stack along y axis (y) Dimension of init wave is 3. (3) Enter the choice (z, y, or 3(default)): z Choose a method for stimulation. Cutting half (h) Stimulting the center (c) No stimulus (n) Enter a method for stimulation (h, c(default), or n): n Enter the time for the stimulus: 20 Enter the strength of the stimulus: 0.8 Enter the radius of the stimulus: 2