#IN CASE GLOBAL NODES ARE STORED IN BLOCK
vrowinc=0
#Assembly Matrix
for vrow in range(numvars): #Variables
#print("vrow",vrow)
ir=0
imax=int(var_dim[vrow])
for n in range (4): #Nodes
for i in range(imax):
d=elnodes.astype(int)[e][n]
#print("ir glob",ir, vrowinc+var_dim[vrow]*d+i)
vnrow[ir]=vrowinc+var_dim[vrow]*d+i
ir=ir+1
# print("vnrow",vnrow.astype(int))
if vrow == 0:
for row in range(4*imax):
UV[row,0]=Uglob[int(vnrow[row])]
elif vrow == 1:
for row in range(4*imax):
UF[row,0]=Uglob[int(vnrow[row])]
vrowinc+=numnodes*var_dim[vrow]