// genesis

// Setting the Resistance of the Gap Junction (1/conductance)
float Rgap = 4.0e9	// Ohms (1/Mhos) from Traub 2005 paper

//float	P5IBb_SEPX = 25e-6
//float	P5IBb_SEPY = 25e-6
//float	P5IBb_NX = 5
//float	P5IBb_NY = 5

int i,j,dexN,dexE,dexS,dexW
int conn1dex, conn2dex
float test1, test2
str conn1loc="", conn2loc=""
float globalconnprob=0.0034 //after Traub 2005
float randconn
int jnodedex,inodedex,jnodecalc

echo P5IBb - P5IBc Gap Junctions

for (i=0;i<=(P5IBb_NY-1);i=i+1)
    for (j=0;j<=(P5IBc_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*P5IBb_NX+j
            conn2dex=(i-1)*P5IBc_NX+j

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}  RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm

        end

        if ({dexN==-1} && {randconn<{globalconnprob}} && { {mynode} >= {sqrtNnodes} }) //for nodal edge
            conn1dex=i*P5IBb_NX+j
            conn2dex=(P5IBc_NY-1)*P5IBc_NX+j

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm 

        end

        // NE
        randconn = { rand 0 1.0 }

        if ({dexN>=0} && {dexE<P5IBb_NX} && {randconn<{globalconnprob}})
            conn1dex=i*P5IBb_NX+j
            conn2dex=(i-1)*P5IBc_NX+j+1

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm

        end

        float jnodecalcarg = {mynode}/{sqrtNnodes}

        jnodecalc = {trunc {jnodecalcarg} }

        jnodedex = {mynode} - {sqrtNnodes}*jnodecalc

        if ({dexN==-1} && {randconn<{globalconnprob}} && {dexE<P5IBb_NX} && { {mynode} >= {sqrtNnodes} }) //for nodal edge
            conn1dex=i*P5IBb_NX+j
            conn2dex=(P5IBc_NY-1)*P5IBc_NX+j+1

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        if ({dexN==-1} && {randconn<{globalconnprob}} && {dexE==P5IBb_NX} && { {mynode} >= {sqrtNnodes} } && {jnodedex} < {sqrtNnodes-1}) //for nodal edge corner

            conn1dex=i*P5IBb_NX+j
            conn2dex=(P5IBc_NY-1)*P5IBc_NX

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - {sqrtNnodes} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - {sqrtNnodes} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
        end

        if ( {dexE==P5IBb_NX} && {dexN>-1} && {randconn<{globalconnprob}} && { {jnodedex} < {sqrtNnodes-1} } ) //for nodal edge

            conn1dex=i*P5IBb_NX+j
            conn2dex=(i-1)*P5IBc_NX

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
        end

        // E
        randconn = { rand 0 1.0 }

        if ({dexE<P5IBb_NX} && {randconn<{globalconnprob}})
            conn1dex=i*P5IBb_NX+j
            conn2dex=i*P5IBc_NX+j+1

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm

        end

        float jnodecalcarg = {mynode}/{sqrtNnodes}

        jnodecalc = {trunc {jnodecalcarg} }

        jnodedex = {mynode} - {sqrtNnodes}*jnodecalc

        if ({dexE==P5IBb_NX} && {randconn<{globalconnprob}} && { {jnodedex} < {sqrtNnodes-1} }) //for nodal edge

            conn1dex=i*P5IBb_NX+j
            conn2dex=i*P5IBc_NX

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        // SE
        randconn = { rand 0 1.0 }

        //pk if ({dexE<P5IBb_NX} && {dexS>P5IBb_NY} && {randconn<{globalconnprob}})
        if ({dexE<P5IBb_NX} && {dexS<P5IBb_NY} && {randconn<{globalconnprob}})

            conn1dex=i*P5IBb_NX+j
            conn2dex=(i+1)*P5IBc_NX+j+1

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm

        end

        float jnodecalcarg = {mynode}/{sqrtNnodes}

        jnodecalc = {trunc {jnodecalcarg} }

        jnodedex = {mynode} - {sqrtNnodes}*jnodecalc

        if ({dexE==P5IBb_NX} && {randconn<{globalconnprob}} && {dexS<P5IBb_NY} && { {jnodedex} < {sqrtNnodes-1} }) //for nodal edge

            conn1dex=i*P5IBb_NX+j
            conn2dex=(i+1)*P5IBc_NX

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end
        //pk if ({dexS==P5IBb_NY} && {randconn<{globalconnprob}} && {dexE==P5IBb_NX} && { {jnodedex} < {sqrtNnodes-1} }) //for nodal edge corner
        if ({dexS==P5IBb_NY} && {randconn<{globalconnprob}} && {dexE==P5IBb_NX} && { {jnodedex} < {sqrtNnodes-1} } && {{mynode} < {{sqrtNnodes}*{sqrtNnodes-1}}}) //for nodal edge corner

            conn1dex=i*P5IBb_NX+j
            conn2dex=0

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
           //pk raddmsg@{{mynode} + {sqrtNnodes} +1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + {sqrtNnodes} +1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        if ( {dexE<P5IBb_NX} && {dexS==P5IBb_NY} && {randconn<{globalconnprob}} && { {mynode} < {{sqrtNnodes}*{sqrtNnodes-1}} } ) //for nodal edge

            conn1dex=i*P5IBb_NX+j
            conn2dex=(j+1)

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} + {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
        end

        // S
        randconn = { rand 0 1.0 }

        if ({dexS<P5IBb_NY} && {randconn<{globalconnprob}})
            conn1dex=i*P5IBb_NX+j
            conn2dex=(i+1)*P5IBc_NX+j

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm

        end

        //pk if ({dexS==P5IBb_NY} && {randconn<{globalconnprob}} && { {mynode} >= {sqrtNnodes}*{sqrtNnodes-1} }) //for nodal edge
        if ({dexS==P5IBb_NY} && {randconn<{globalconnprob}} && { {mynode} < {sqrtNnodes}*{sqrtNnodes-1} }) //for nodal edge    
            conn1dex=i*P5IBb_NX+j
            conn2dex=j

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} + {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
        end

        // SW
        randconn = { rand 0 1.0 }

        if ({dexS<P5IBb_NY} && {dexW>=0} && {randconn<{globalconnprob}})
            conn1dex=i*P5IBb_NX+j
            conn2dex=(i+1)*P5IBc_NX+j-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} RAXIAL {Rgap} Vm

        end

        float jnodecalcarg = {mynode}/{sqrtNnodes}

        jnodecalc = {trunc {jnodecalcarg} }

        jnodedex = {mynode} - {sqrtNnodes}*jnodecalc

        if ({dexE==1} && {randconn<{globalconnprob}} && {dexS<P5IBb_NY} && { {jnodedex} > 0 }) //for nodal edge

            conn1dex=i*P5IBb_NX+j
            conn2dex=(i+2)*P5IBc_NX-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        if ({dexS==P5IBb_NY} && {randconn<{globalconnprob}} && {dexE==1} && { {jnodedex} > 0 } && { {mynode} < {{sqrtNnodes}*{sqrtNnodes-1}}}) //for nodal edge corner

            conn1dex=i*P5IBb_NX+j
            conn2dex=P5IBc_NX-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} + {sqrtNnodes} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + {sqrtNnodes} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        if ( {dexE>1} && {dexS==P5IBb_NY} && {randconn<{globalconnprob}} &&{ {mynode} < {{sqrtNnodes}*{sqrtNnodes-1}} } ) //for nodal edge

            conn1dex=i*P5IBb_NX+j
            conn2dex=(j-1)

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} + {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} + {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        // W
        randconn = { rand 0 1.0 }

        if ({dexW>=0} && {randconn<{globalconnprob}})
            conn1dex=i*P5IBb_NX+j
            conn2dex=i*P5IBc_NX+j-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{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*P5IBb_NX+j
            conn2dex=(i+1)*P5IBc_NX-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm
        end

        // NW
        randconn = { rand 0 1.0 }

        if ({dexW>=0} && {dexN>=0} && {randconn<{globalconnprob}})
            conn1dex=i*P5IBb_NX+j
            conn2dex=(i-1)*P5IBc_NX+j-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            addmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} RAXIAL {Rgap} Vm
            addmsg /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{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*P5IBb_NX+j
            conn2dex=(P5IBc_NY-1)*P5IBc_NX+j-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes}} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes}} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - {sqrtNnodes}} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        if ({dexN==-1} && {randconn<{globalconnprob}} && {dexE==1} && { {mynode} >= {sqrtNnodes} } && { {jnodedex} > 0 }) //for nodal edge corner

            conn1dex=i*P5IBb_NX+j
            conn2dex=(P5IBc_NY*P5IBc_NX)-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - {sqrtNnodes} - 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - {sqrtNnodes} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + {sqrtNnodes} + 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - {sqrtNnodes} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

        if ( {dexE==1} && {dexN>-1} && {randconn<{globalconnprob}} && { {jnodedex} > 0 } ) //for nodal edge

            conn1dex=i*P5IBb_NX+j
            conn2dex=i*P5IBc_NX-1

            conn1loc="axonLa"
            conn2loc="axonRa"

            raddmsg /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc}@{{mynode} - 1} RAXIAL {Rgap} Vm
            //pk raddmsg@{{mynode} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{{mynode} + 1} RAXIAL {Rgap} Vm
            raddmsg@{{mynode} - 1} /P5IBcnet/P5IBc[{conn2dex}]/{conn2loc} /P5IBbnet/P5IBb[{conn1dex}]/{conn1loc}@{mynode} RAXIAL {Rgap} Vm

        end

    end
end