// genesis
// Setting the Resistance of the Gap Junction (1/conductance)
float Rgap = 4.0e9 // Ohms (1/Mhos) from Traub 2005 paper
//float P6RSb_SEPX = 25e-6
//float P6RSb_SEPY = 25e-6
//float P6RSb_NX = 5 ;defined in netparams.g
//float P6RSb_NY = 5 ;defined in netparams.g
int i,j,dexN,dexE,dexS,dexW
int conn1dex, conn2dex
float test1, test2
str conn1loc="", conn2loc=""
float randconn
int jnodedex,inodedex,jnodecalc
if ({columntype == 1})
float globalconnprob=0.0156 //after Traub 2005
end
if ({columntype == 0})
float globalconnprob=0.0078 //after Traub 2005
end
echo P6RSb - P6RSb Gap Junctions
for (i=0;i<=(P6RSb_NY-1);i=i+1)
for (j=0;j<=(P6RSb_NX-1);j=j+1)
dexN=i-1
dexE=j+1
dexS=i+1
dexW=j-1
// N
randconn = { rand 0 1.0 }
if ({dexN>=0} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=(i-1)*P6RSb_NX+j
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
if ({dexN==-1} && {randconn<{globalconnprob}} && { {mynode} >= {sqrtNnodes} }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(P6RSb_NY-1)*P6RSb_NX+j
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
// NE
randconn = { rand 0 1.0 }
if ({dexN>=0} && {dexE<P6RSb_NX} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=(i-1)*P6RSb_NX+j+1
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
float jnodecalcarg = {mynode}/{sqrtNnodes}
jnodecalc = {trunc {jnodecalcarg} }
jnodedex = {mynode} - {sqrtNnodes}*jnodecalc
if ({dexN==-1} && {randconn<{globalconnprob}} && {dexE<P6RSb_NX} && { {mynode} >= {sqrtNnodes} }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(P6RSb_NY-1)*P6RSb_NX+j+1
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
if ({dexN==-1} && {randconn<{globalconnprob}} && {dexE==P6RSb_NX} && { {mynode} >= {sqrtNnodes} } && {jnodedex} < {sqrtNnodes-1}) //for nodal edge corner
conn1dex=i*P6RSb_NX+j
conn2dex=(P6RSb_NY-1)*P6RSb_NX
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - {sqrtNnodes} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - {sqrtNnodes} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
if ( {dexE==P6RSb_NX} && {dexN>-1} && {randconn<{globalconnprob}} && { {jnodedex} < {sqrtNnodes-1} } ) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(i-1)*P6RSb_NX
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
// E
randconn = { rand 0 1.0 }
if ({dexE<P6RSb_NX} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=i*P6RSb_NX+j+1
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
float jnodecalcarg = {mynode}/{sqrtNnodes}
jnodecalc = {trunc {jnodecalcarg} }
jnodedex = {mynode} - {sqrtNnodes}*jnodecalc
if ({dexE==P6RSb_NX} && {randconn<{globalconnprob}} && { {jnodedex} < {sqrtNnodes-1} }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=i*P6RSb_NX
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
// SE
randconn = { rand 0 1.0 }
//pk if ({dexE<P6RSb_NX} && {dexS>P6RSb_NY} && {randconn<{globalconnprob}})
if ({dexE<P6RSb_NX} && {dexS<P6RSb_NY} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=(i+1)*P6RSb_NX+j+1
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
float jnodecalcarg = {mynode}/{sqrtNnodes}
jnodecalc = {trunc {jnodecalcarg} }
jnodedex = {mynode} - {sqrtNnodes}*jnodecalc
if ({dexE==P6RSb_NX} && {randconn<{globalconnprob}} && {dexS<P6RSb_NY} && { {jnodedex} < {sqrtNnodes-1} }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(i+1)*P6RSb_NX
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
//pk if ({dexS==P6RSb_NY} && {randconn<{globalconnprob}} && {dexE==P6RSb_NX} && { {jnodedex} < {{sqrtNnodes}-1} } ) //for nodal edge corner
if ({dexS==P6RSb_NY} && {randconn<{globalconnprob}} && {dexE==P6RSb_NX} && { {jnodedex} < {sqrtNnodes-1} } && {{mynode} < {{sqrtNnodes}*{sqrtNnodes-1}}}) //for nodal edge corner
conn1dex=i*P6RSb_NX+j
conn2dex=0
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + {sqrtNnodes} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + {sqrtNnodes} + 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
if ( {dexE<P6RSb_NX} && {dexS==P6RSb_NY} && {randconn<{globalconnprob}} && { {mynode} < {{sqrtNnodes}*{sqrtNnodes-1}} } ) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(j+1)
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
// S
randconn = { rand 0 1.0 }
if ({dexS<P6RSb_NY} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=(i+1)*P6RSb_NX+j
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
//pk if ({dexS==P6RSb_NY} && {randconn<{globalconnprob}} && { {mynode} >= {sqrtNnodes}*{sqrtNnodes-1} }) //for nodal edge
if ({dexS==P6RSb_NY} && {randconn<{globalconnprob}} && { {mynode} < {sqrtNnodes}*{sqrtNnodes-1} }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=j
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
// SW
randconn = { rand 0 1.0 }
if ({dexS<P6RSb_NY} && {dexW>=0} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=(i+1)*P6RSb_NX+j-1
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
float jnodecalcarg = {mynode}/{sqrtNnodes}
jnodecalc = {trunc {jnodecalcarg} }
jnodedex = {mynode} - {sqrtNnodes}*jnodecalc
if ({dexE==1} && {randconn<{globalconnprob}} && {dexS<P6RSb_NY} && { {jnodedex} > 0 }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(i+2)*P6RSb_NX-1
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
if ({dexS==P6RSb_NY} && {randconn<{globalconnprob}} && {dexE==1} && { {jnodedex} > 0 } && { {mynode} < {{sqrtNnodes}*{sqrtNnodes-1}}}) //for nodal edge corner
conn1dex=i*P6RSb_NX+j
conn2dex=P6RSb_NX-1
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + {sqrtNnodes} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + {sqrtNnodes} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
if ( {dexE>1} && {dexS==P6RSb_NY} && {randconn<{globalconnprob}} &&{ {mynode} < {{sqrtNnodes}*{sqrtNnodes-1}} } ) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(j-1)
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} + {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
raddmsg@{{mynode} + {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
// W
randconn = { rand 0 1.0 }
if ({dexW>=0} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=i*P6RSb_NX+j-1
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
float jnodecalcarg = {mynode}/{sqrtNnodes}
jnodecalc = {trunc {jnodecalcarg} }
jnodedex = {mynode} - {sqrtNnodes}*jnodecalc
if ({dexW==-1} && {randconn<{globalconnprob}} && { {jnodedex} > 0 }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(i+1)*P6RSb_NX-1
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
// NW
randconn = { rand 0 1.0 }
if ({dexW>=0} && {dexN>=0} && {randconn<{globalconnprob}})
conn1dex=i*P6RSb_NX+j
conn2dex=(i-1)*P6RSb_NX+j-1
conn1loc="axonLa"
conn2loc="axonRa"
addmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
addmsg /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm
end
float jnodecalcarg = {mynode}/{sqrtNnodes}
jnodecalc = {trunc {jnodecalcarg} }
jnodedex = {mynode} - {sqrtNnodes}*jnodecalc
if ({dexN==-1} && {randconn<{globalconnprob}} && {dexE>1} && { {mynode} >= {sqrtNnodes} }) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=(P6RSb_NY-1)*P6RSb_NX+j-1
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - {sqrtNnodes}} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
if ({dexN==-1} && {randconn<{globalconnprob}} && {dexE==1} && { {mynode} >= {sqrtNnodes} } && { {jnodedex} > 0 }) //for nodal edge corner
conn1dex=i*P6RSb_NX+j
conn2dex=(P6RSb_NY*P6RSb_NX)-1
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - {sqrtNnodes} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - {sqrtNnodes} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
if ( {dexE==1} && {dexN>-1} && {randconn<{globalconnprob}} && { {jnodedex} > 0 } ) //for nodal edge
conn1dex=i*P6RSb_NX+j
conn2dex=i*P6RSb_NX-1
conn1loc="axonLa"
conn2loc="axonRa"
raddmsg /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
//pk raddmsg@{{mynode} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
raddmsg@{{mynode} - 1} /P6RSbnet/P6RSb[{conn2dex}]/{conn2loc} /P6RSbnet/P6RSb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
end
end
end
| 
