fixed integer division problem in python code snippet

for j1 in range(0,twojmax+1):

for j1 in range(0,twojmax+1):

    if(diagonal==2): 

    if(diagonal==2): 

        print j1/2,j1/2,j1/2

        print j1/2.,j1/2.,j1/2.

    elif(diagonal==1):

    elif(diagonal==1):

        for j in range(0,min(twojmax,2*j1)+1,2): 

        for j in range(0,min(twojmax,2*j1)+1,2): 

            print j1/2,j1/2,j/2

            print j1/2.,j1/2.,j/2.

    elif(diagonal==0):

    elif(diagonal==0):

        for j2 in range(0,j1+1):

        for j2 in range(0,j1+1):

            for j in range(j1-j2,min(twojmax,j1+j2)+1,2): 

            for j in range(j1-j2,min(twojmax,j1+j2)+1,2): 

                print j1/2,j2/2,j/2

                print j1/2.,j2/2.,j/2.

    elif(diagonal==3):

    elif(diagonal==3):

        for j2 in range(0,j1+1):

        for j2 in range(0,j1+1):

            for j in range(j1-j2,min(twojmax,j1+j2)+1,2): 

            for j in range(j1-j2,min(twojmax,j1+j2)+1,2): 

                if (j>=j1): print j1/2,j2/2,j/2 :pre

                if (j>=j1): print j1/2.,j2/2.,j/2. :pre

Compute {snad/atom} evaluates a per-atom array. The columns are

Compute {snad/atom} evaluates a per-atom array. The columns are

arranged into {ntypes} blocks, listed in order of atom type {I}.  Each

arranged into {ntypes} blocks, listed in order of atom type {I}.  Each