Sorry, I don't follow what you're saying. I have a 3D rectangular mesh with
dimensions (JD,KD,LD), so there are JD*KD*LD points and (JD-1)*(KD-1)*(LD-1)
cells. I have arrays X,Y,Z dimensioned (JD,KD,LD). a cell-centered array
WD dimensioned (JD-1,KD-1,LD-1), and the connectivity array ELT2NODE dimensioned
(8,JD-1,KD-1,LD-1). I have NCELLS set to (JD-1)*(KD-1)*(LD-1) and NPTS
set to JD*KD*LD. I set up the arrays X,Y,Z,WD, ELT2NODE and have the calls
I = TECINI112( 'Centroid Z coordinate'//C_NULL_CHAR,...
I = TECZNE112 ( 'Simple brick mesh'//C_NULL_CHAR,...
I = TECDAT112 ( NPTS,X,ISDOUBLE )
I = TECDAT112 ( NPTS,Y,ISDOUBLE )
I = TECDAT112 ( NPTS,Z,ISDOUBLE )
I = TECDAT112 ( NCELLS,WD,ISDOUBLE )
I = TECNODE112 ( NCELLS,ELT2NODE )
I = TECEND112 ( )
I'm not sure where a multiplication by 8 is supposed to occur.