volumetric weighted average of SDV1 using python

I've been trying to obtain volumetric weighted average of SDV1 in my model by
using element volume EVOL. I wrote the following python code to do it for
me, but it gives an error on multiplication part (el_VOL*VR) "incompatible
data locations". I know SDV1 is stored at integration point and EVOL at
whole element level and that is probably why the error. I probably should
use IVOL; however, I already finished simulations and didn't store IVOL as I
thought EVOL is enough, especially since the elements are C3D8R therefore
EVOL and IVOL should be the same numbers. Does anyone know how to manipulate
this code to perform the multiplications?

odbPath = "C:\\Temp\\cap03.odb"
odb = session.openOdb(name=odbPath,readOnly=FALSE)
allSteps = session.odbData[odbPath].steps.keys()
step = odb.steps[allSteps[1]]
allFrames = session.odbData[odbPath].steps[allSteps[1]].frames.keys()
myInstance = odb.rootAssembly.instances["PART-1-1"]
for j in range(len(allFrames)):
	frame = step.frames[j]
	VR = frame.fieldOutputs['SDV1'].getSubset(region=myInstance)
	el_VOL = frame.fieldOutputs['EVOL'].getSubset(region=myInstance)
	e_times_V = el_VOL*VR
	newField = frame.FieldOutput(name='e_times_V', description='NA',
field=e_times_V)
	print 'stepName = ', allSteps[1], ' frameNumber = ', allFrames[j]
odb.save()
odb.close()



--
Sent from: http://abaqus-users.1086179.n5.nabble.com/