Learning Forum

[David Weed]

Hi,

The following code will give you the element volume fraction (volume per integration point):

DOUBLE PRECISION vect2(1)
INTEGER v_num

v_num = 1

 

print*,'element num',elemId

print*,'int point num',kDomIntPt

call get_ElmData ('EJAC',elemId, kDomIntPt,v_num,vect2)

print*,'vect2 value',vect2