Skip to content

Commit b25177f

Browse files
gsjaardemagsjaardema
authored
Merge pull request #644 from sandialabs/NUMBERS_support_tet_meshes
Numbers support tet meshes
2 parents fde7419 + d45faab commit b25177f

28 files changed

+605
-98
lines changed

‎packages/seacas/applications/numbers/nu_cavity.f‎

Lines changed: 19 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -1,13 +1,15 @@
1-
C Copyright(C) 1999-2020 National Technology & Engineering Solutions
1+
C Copyright(C) 1999-2020, 2025 National Technology & Engineering Solutions
22
C of Sandia, LLC (NTESS). Under the terms of Contract DE-NA0003525 with
33
C NTESS, the U.S. Government retains certain rights in this software.
44
C
55
C See packages/seacas/LICENSE for details
66

77
SUBROUTINE CAVITY (A, CRD, IDESS, NEESS, NNESS, IPEESS, IPNESS,
88
* LTEESS, LTNESS, FACESS, DISP, NUMNP, NDIM, NUMESS,
9-
* TIME, ITMSEL, TITLE, CENT, CENTER)
9+
* TIME, ITMSEL, TITLE, CENT, CENTER, NNODES)
1010

11+
C .. nnodes is number of nodes per element
12+
1113
include 'nu_io.blk'
1214
DIMENSION A(*), CRD(NUMNP,NDIM), IDESS(*), NEESS(*),
1315
* NNESS(*), IPEESS(*), IPNESS(*), LTEESS(*), LTNESS(*),
@@ -27,8 +29,14 @@ SUBROUTINE CAVITY (A, CRD, IDESS, NEESS, NNESS, IPEESS, IPNESS,
2729
IFLG = IFND(NCAV)
2830
IPTR = IPNESS(IFLG)
2931
IF (NDIM .EQ. 3) THEN
30-
CALL VOL3D( CRD, LTNESS(IPTR), NEESS(IFLG), VOLUME,
31-
* NDIM, NUMESS, CENT, NUMNP, CENTER)
32+
if (nnodes .eq. 8) then
33+
CALL VOL3D( CRD, LTNESS(IPTR), NEESS(IFLG), VOLUME,
34+
* NDIM, NUMESS, CENT, NUMNP, CENTER)
35+
endif
36+
if (nnodes .eq. 4) then
37+
CALL TVOL3D( CRD, LTNESS(IPTR), NEESS(IFLG), VOLUME,
38+
* NDIM, NUMESS, CENT, NUMNP, CENTER)
39+
endif
3240
ELSE
3341
CALL VOL2D( CRD, LTNESS(IPTR), NEESS(IFLG), VOLUME,
3442
* NDIM, NUMESS, AXI, CENT, NUMNP, CENTER)
@@ -79,8 +87,13 @@ SUBROUTINE CAVITY (A, CRD, IDESS, NEESS, NNESS, IPEESS, IPNESS,
7987
C NOTE: Positive delcav = shrink in cavity volume
8088

8189
IF (NDIM .EQ. 3) THEN
82-
CALL DVOL3D(CRD, DISP, LTNESS(IPTR),
83-
* NEESS(IFLG), DELCAV, NDIM, NUMNP)
90+
if (nnodes .eq. 8) then
91+
CALL DVOL3D(CRD, DISP, LTNESS(IPTR),
92+
* NEESS(IFLG), DELCAV, NDIM, NUMNP)
93+
else
94+
CALL DTVOL3D(CRD, DISP, LTNESS(IPTR),
95+
* NEESS(IFLG), DELCAV, NDIM, NUMNP)
96+
end if
8497
ELSE
8598
CALL DVOL2D(CRD, DISP, LTNESS(IPTR),
8699
* NEESS(IFLG), DELCAV, NDIM, AXI, NUMNP)

‎packages/seacas/applications/numbers/nu_cgcal2.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,4 @@
1-
C Copyright(C) 1999-2020 National Technology & Engineering Solutions
1+
C Copyright(C) 1999-2020, 2025 National Technology & Engineering Solutions
22
C of Sandia, LLC (NTESS). Under the terms of Contract DE-NA0003525 with
33
C NTESS, the U.S. Government retains certain rights in this software.
44
C
@@ -8,7 +8,7 @@ SUBROUTINE CGCAL2(CRD,IX,MAT,MASS,VOL,DENS,VOLM,CG,ZITOT,XXX,
88
* XG,XI,XINI,AJ,NNODES,NDIM,NQUAD,VOLMN,IELM,NELBLK,
99
* AXI,NUMNP)
1010

11-
DIMENSION CRD(NUMNP,*), IX(NNODES,*), MAT(6,*), MASS(*),
11+
DIMENSION CRD(NUMNP,*), IX(NNODES,*), MAT(7,*), MASS(*),
1212
* DENS(*), VOLM(*), CG(*), ZITOT(*),VOLMN(4,*),IELM(4,*),
1313
* XXX(NDIM+1,NQUAD,*),XG(NDIM,*), XI(NDIM,*), XINI(*),
1414
* AJ(2,*)

‎packages/seacas/applications/numbers/nu_cgcal3.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,4 @@
1-
C Copyright(C) 1999-2020 National Technology & Engineering Solutions
1+
C Copyright(C) 1999-2020, 2025 National Technology & Engineering Solutions
22
C of Sandia, LLC (NTESS). Under the terms of Contract DE-NA0003525 with
33
C NTESS, the U.S. Government retains certain rights in this software.
44
C
@@ -9,7 +9,7 @@ SUBROUTINE CGCAL3(CRD,IX,MAT,MASS,VOL,DENS,VOLM,CG,ZITOT,XXX,
99

1010
C ... CALCULATE PROPERTIES FOR THREE-DIMENSIONAL MESH --- BRICKS ONLY
1111

12-
DIMENSION CRD(NUMNP,*), IX(NNODES,*), MAT(6,*), MASS(*),
12+
DIMENSION CRD(NUMNP,*), IX(NNODES,*), MAT(7,*), MASS(*),
1313
* DENS(*), VOLM(*), CG(*), ZITOT(*),VOLMN(4,*),IELM(4,*),
1414
* XXX(3+1,NQUAD,*),XG(3,*), XI(3,*), XINI(*),
1515
* AJ(3,*)

‎packages/seacas/applications/numbers/nu_cgcal3t.f‎

Lines changed: 164 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,164 @@
1+
C Copyright(C) 2025 National Technology & Engineering Solutions
2+
C of Sandia, LLC (NTESS). Under the terms of Contract DE-NA0003525 with
3+
C NTESS, the U.S. Government retains certain rights in this software.
4+
C
5+
C See packages/seacas/LICENSE for details
6+
7+
SUBROUTINE CGCAL3T(CRD,IX,MAT,MASS,VOL,DENS,VOLM,CG,ZITOT,XXX,
8+
* XG,XI,XINI,AJ,NNODES,NDIM,NQUAD,VOLMN,IELM,NELBLK,NUMNP)
9+
10+
C ... CALCULATE PROPERTIES FOR THREE-DIMENSIONAL TET MESH
11+
12+
DIMENSION CRD(NUMNP,*), IX(4,*), MAT(7,*), MASS(*),
13+
* DENS(*), VOLM(*), CG(*), ZITOT(*),VOLMN(4,*),IELM(4,*),
14+
* XXX(3+1,NQUAD,*),XG(3,*), XI(3,*), XINI(*),
15+
* AJ(3,*)
16+
DIMENSION ZI(6), ZMOM(3)
17+
DIMENSION CCC(3,8)
18+
19+
REAL MASS, MASSE
20+
REAL M(3,3)
21+
22+
C ... VOLMN(1,*) = MINIMUM VOLUME (AREAS FOR 2-D)
23+
C VOLMN(2,*) = MAXIMUM VOLUME
24+
C VOLMN(3,*) = TOTAL VOLUME
25+
C VOLMN(4,*) = MAXIMUM SUM OF SQUARES OF INVERSE LENGTHS
26+
27+
DO I=1, NELBLK
28+
VOLMN(1,I) = 1.0E30
29+
VOLMN(2,I) = 0.0
30+
VOLMN(3,I) = 0.0
31+
IELM (3,I) = 0
32+
MASS(I) = 0.0
33+
VOLM(I) = 0.0
34+
END DO
35+
36+
DO I=1,3
37+
ZITOT(I) = 0.0
38+
ZITOT(I+3) = 0.0
39+
ZMOM(I) = 0.0
40+
END DO
41+
42+
ZMAS = 0.0
43+
VOL = 0.0
44+
45+
C ... GET QUADRATURE POINT LOCATIONS, EVALUATE SHAPE FUNCTIONS
46+
47+
DO IBLK = 1, NELBLK
48+
IF (MAT(5,IBLK) .NE. 1) CONTINUE
49+
IELBEG = MAT(3,IBLK)
50+
IELEND = MAT(4,IBLK)
51+
MIEL = IBLK
52+
DO IEL = IELBEG, IELEND
53+
54+
C ... CALCULATE AREA, VOLUME, AND MOMENTS OF INERTIA OF ELEMENT
55+
56+
DO I=1,3
57+
ZI(I) = 0.0
58+
ZI(I+3) = 0.0
59+
CG(I) = 0.0
60+
END DO
61+
VOLUME = 0.0
62+
63+
X1 = CRD(IX(1,IEL),1)
64+
X2 = CRD(IX(2,IEL),1)
65+
X3 = CRD(IX(3,IEL),1)
66+
X4 = CRD(IX(4,IEL),1)
67+
68+
Y1 = CRD(IX(1,IEL),2)
69+
Y2 = CRD(IX(2,IEL),2)
70+
Y3 = CRD(IX(3,IEL),2)
71+
Y4 = CRD(IX(4,IEL),2)
72+
73+
Z1 = CRD(IX(1,IEL),3)
74+
Z2 = CRD(IX(2,IEL),3)
75+
Z3 = CRD(IX(3,IEL),3)
76+
Z4 = CRD(IX(4,IEL),3)
77+
78+
xini(1) = (x1 + x2 + x3 + x4) / 4.0
79+
xini(2) = (y1 + y2 + y3 + y4) / 4.0
80+
xini(3) = (z1 + z2 + z3 + z4) / 4.0
81+
82+
m(1,1) = x1 - x2
83+
m(1,2) = y1 - y2
84+
m(1,3) = z1 - z2
85+
86+
m(2,1) = x2 - x3
87+
m(2,2) = y2 - y3
88+
m(2,3) = z2 - z3
89+
90+
m(3,1) = x3 - x4
91+
m(3,2) = y3 - y4
92+
m(3,3) = z3 - z4
93+
94+
det = (m(1,1) * (m(2,2) * m(3,3) - m(2,3) * m(3,2)) -
95+
$ m(2,1) * (m(1,2) * m(3,3) - m(1,3) * m(3,2)) +
96+
$ m(3,1) * (m(1,2) * m(2,3) - m(1,3) * m(2,2)))
97+
98+
det = -det / 6.0
99+
100+
DETW = DET * DENS(MIEL)
101+
102+
CG(1) = CG(1) + DETW * XINI(1)
103+
CG(2) = CG(2) + DETW * XINI(2)
104+
CG(3) = CG(3) + DETW * XINI(3)
105+
106+
ZI(1) = ZI(1) + DETW * XINI(1)**2
107+
ZI(2) = ZI(2) + DETW * XINI(2)**2
108+
ZI(3) = ZI(3) + DETW * XINI(3)**2
109+
110+
ZI(4) = ZI(4) + DETW * XINI(1) * XINI(2)
111+
ZI(5) = ZI(5) + DETW * XINI(1) * XINI(3)
112+
ZI(6) = ZI(6) + DETW * XINI(2) * XINI(3)
113+
114+
VOLUME = VOLUME + DET
115+
116+
117+
C ... DETERMINE MIN/MAX ELEMENT VOLUMES FOR EACH MATERIAL AND
118+
C COUNT NUMBER OF ELEMENTS FOR EACH MATERIAL
119+
120+
IELM(3,MIEL) = IELM(3,MIEL) + 1
121+
VOLMN(3,MIEL) = VOLMN(3,MIEL) + VOLUME
122+
IF (VOLUME .LT. VOLMN(1,MIEL)) THEN
123+
VOLMN(1,MIEL) = VOLUME
124+
IELM(1,MIEL) = IEL
125+
END IF
126+
C ... Changed from else if to if so 1 element and equal size blocks get correct volume
127+
IF (VOLUME .GT. VOLMN(2,MIEL)) THEN
128+
VOLMN(2,MIEL) = VOLUME
129+
IELM(2,MIEL) = IEL
130+
END IF
131+
132+
if (volume .le. 0.0) then
133+
write (*,*) 'Zero or negative volume at element',
134+
& iel
135+
end if
136+
ZITOT(1) = ZITOT(1) + ZI(2) + ZI(3)
137+
ZITOT(2) = ZITOT(2) + ZI(1) + ZI(3)
138+
ZITOT(3) = ZITOT(3) + ZI(1) + ZI(2)
139+
ZITOT(4) = ZITOT(4) + ZI(4)
140+
ZITOT(5) = ZITOT(5) + ZI(5)
141+
ZITOT(6) = ZITOT(6) + ZI(6)
142+
MASSE = VOLUME * DENS(MIEL)
143+
MASS(MIEL)= MASS(MIEL) + MASSE
144+
VOLM(MIEL)= VOLM(MIEL) + VOLUME
145+
VOL = VOL + VOLUME
146+
ZMAS = ZMAS + MASSE
147+
ZMOM(1) = ZMOM(1) + CG(1)
148+
ZMOM(2) = ZMOM(2) + CG(2)
149+
ZMOM(3) = ZMOM(3) + CG(3)
150+
END DO
151+
END DO
152+
FIX = SIGN(0.5, ZMAS) + SIGN(0.5, -ZMAS)
153+
DO I=1,3
154+
CG(I) = ZMOM(I) / (ZMAS + FIX)
155+
END DO
156+
ZITOT(1) = ZITOT(1) - ZMAS * (CG(2)**2 + CG(3)**2)
157+
ZITOT(2) = ZITOT(2) - ZMAS * (CG(1)**2 + CG(3)**2)
158+
ZITOT(3) = ZITOT(3) - ZMAS * (CG(1)**2 + CG(2)**2)
159+
ZITOT(4) = ZITOT(4) - ZMAS * CG(1) * CG(2)
160+
ZITOT(5) = ZITOT(5) - ZMAS * CG(1) * CG(3)
161+
ZITOT(6) = ZITOT(6) - ZMAS * CG(2) * CG(3)
162+
163+
RETURN
164+
END

‎packages/seacas/applications/numbers/nu_command.f‎

Lines changed: 32 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,4 @@
1-
C Copyright(C) 1999-2020, 2022 National Technology & Engineering Solutions
1+
C Copyright(C) 1999-2020, 2022, 2025 National Technology & Engineering Solutions
22
C of Sandia, LLC (NTESS). Under the terms of Contract DE-NA0003525 with
33
C NTESS, the U.S. Government retains certain rights in this software.
44
C
@@ -14,7 +14,7 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
1414
include 'exodusII.inc'
1515
CHARACTER*80 TITLE, COMMENT
1616
INTEGER IA(*)
17-
DIMENSION A(*), TIME(*), MAT(6,*), DISP(NUMNP,*), CRD(NUMNP,*),
17+
DIMENSION A(*), TIME(*), MAT(7,*), DISP(NUMNP,*), CRD(NUMNP,*),
1818
* LINK(*), DENSTY(*), WAVE(*)
1919
INTEGER ITMSEL(*), ISEVOK(*)
2020
CHARACTER*(MXSTLN) NAMEGL(*), NAMENV(*), NAMEEL(*)
@@ -102,8 +102,8 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
102102

103103
C ... GET SOME SCRATCH SPACE
104104

105-
CALL MDRSRV ('SCRTCH', ISCR, NDIM*NUMNP)
106-
CALL MDRSRV ('SCRTC2', ISCR2, NDIM*NUMNP)
105+
CALL MDRSRV ('SCRTCH', ISCR, NDIM*MAX(NUMEL, NUMNP))
106+
CALL MDRSRV ('SCRTC2', ISCR2, NDIM*MAX(NUMEL, NUMNP))
107107
CALL MDRSRV ('SORTMP', ISMP, MAX(NUMEL, NUMNP) )
108108
CALL MDRSRV ('NODSEL', INDSEL, NUMNP)
109109
CALL MDRSRV ('ELMSEL', IELSEL, NUMEL)
@@ -247,6 +247,11 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
247247
C-----------------------------------------------------------------------
248248
ELSE IF (NAME .EQ. 'TIMESTEP') THEN
249249

250+
if (mat(7,1) .eq. 4) then
251+
CALL PRTERR ('CMDERR',
252+
* 'TIMESTEP command not supported for TET meshes')
253+
GO TO 20
254+
end if
250255
CALL FFREAL (IFLD, KV, RV,
251256
* 'common material wavespeed', 0.0, CWAVE, *20)
252257
IF (CWAVE .GT. 0.) THEN
@@ -310,7 +315,7 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
310315
END IF
311316
CALL CAVITY (A, CRD, A(IBC1), A(IBC2), A(IBC3), A(IBC4),
312317
* A(IBC5),A(IBC6), A(IBC7), A(IBC8), DISP, NUMNP, NDIM,
313-
* NUMESS, TIME, ITMSEL, TITLE, CENT, CENTER)
318+
* NUMESS, TIME, ITMSEL, TITLE, CENT, CENTER, MAT(7,1))
314319

315320
ELSE IF (NAME .EQ. 'MCAVITY') THEN
316321
NUMCAV = 0
@@ -346,10 +351,15 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
346351
END IF
347352
CALL MULTI_CAVITY (A, CRD, A(IBC1), A(IBC2), A(IBC3), A(IBC4),
348353
* A(IBC5),A(IBC6), A(IBC7), A(IBC8), DISP, NUMNP, NDIM,
349-
* NUMESS, TIME, ITMSEL, TITLE, CENT, CENTER)
354+
* NUMESS, TIME, ITMSEL, TITLE, CENT, CENTER, MAT(7,1))
350355

351356
C-----------------------------------------------------------------------
352357
ELSE IF (NAME .EQ. 'GAP') THEN
358+
if (mat(7,1) .eq. 4) then
359+
CALL PRTERR ('CMDERR',
360+
* 'GAP command not supported for TET meshes')
361+
GO TO 20
362+
end if
353363
CALL FFINTG (IFLD, KV, IVAL,
354364
* 'gap master surface', 0, IMAS, *20)
355365
CALL FFINTG (IFLD, KV, IVAL,
@@ -373,6 +383,11 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
373383

374384
C-----------------------------------------------------------------------
375385
ELSE IF (NAME .EQ. 'OVERLAP') THEN
386+
if (mat(7,1) .eq. 4) then
387+
CALL PRTERR ('CMDERR',
388+
* 'OVERLAP command not supported for TET meshes')
389+
GO TO 20
390+
end if
376391
CALL FFINTG (IFLD, KV, IVAL,
377392
* 'overlap master surface', 0, IMAS, *20)
378393
CALL FFINTG (IFLD, KV, IVAL,
@@ -529,12 +544,12 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
529544
END IF
530545

531546
ELSE IF (LISTYP .EQ. 'VOLUME') THEN
532-
CALL PRVOL (ndim, CRD, link, numnp, numel, 8,
547+
CALL PRVOL (ndim, CRD, link, numnp, numel, mat(7,1),
533548
& a(ismp), IHARD)
534549
CALL PRTERR ('CMDSPEC',
535550
* 'Element Volumes were written to the list file')
536551
ELSE IF (LISTYP .EQ. 'NODALVOL') THEN
537-
CALL PRNVOL (ndim, CRD, link, numnp, numel, 8,
552+
CALL PRNVOL (ndim, CRD, link, numnp, numel, mat(7,1),
538553
& a(ismp), IHARD)
539554
CALL PRTERR ('CMDSPEC',
540555
* 'Nodal Volumes were written to the list file')
@@ -739,11 +754,11 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
739754
END IF
740755
IF (ALLTIM) THEN
741756
CALL LIMITS (A(IXYZMN), A(IXYZMX), CRD, LINK, MAT, NDIM,
742-
* NELBLK, 2**NDIM, ALLTIM, TIME, ITMSEL,
757+
* NELBLK, MAT(7,1), ALLTIM, TIME, ITMSEL,
743758
* DISP, NUMNP)
744759
ELSE
745760
CALL LIMITS (A(IXYZMN), A(IXYZMX), CRD, LINK, MAT, NDIM,
746-
* NELBLK, 2**NDIM, ALLTIM, TIME, ITMSEL,
761+
* NELBLK, MAT(7,1), ALLTIM, TIME, ITMSEL,
747762
* CRD, NUMNP)
748763
END IF
749764
CALL MDDEL ('XYZMIN')
@@ -838,6 +853,11 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
838853
END IF
839854
C ----------------------------------------
840855
ELSE IF (NAME .EQ. 'CONDITIO' ) THEN
856+
if (mat(7,1) .eq. 4) then
857+
CALL PRTERR ('CMDERR',
858+
* 'CONDITION command not supported for TET meshes')
859+
GO TO 20
860+
end if
841861
LTMP = .FALSE.
842862
CALL FFCHAR (IFLD, KV, CV, 'NODEBUG', CTMP)
843863
IF (MATSTR(CTMP, 'DEBUG', 1)) LTMP = .TRUE.
@@ -861,13 +881,12 @@ SUBROUTINE COMMAND (A, IA, TITLE, TIME, ITMSEL, MAT, DISP,
861881
CALL MEMERR
862882
STOP 'MEMORY'
863883
END IF
864-
NNODES = 2**NDIM
865884
IF (NDIM .EQ. 2) THEN
866-
CALL CON2D(CRD, NDIM, NUMNP, LINK, NNODES, NUMEL, MAT,
885+
CALL CON2D(CRD, NDIM, NUMNP, LINK, 2**NDIM, NUMEL, MAT,
867886
* NELBLK, A(IELSEL), A(IASPEC), A(ISKEW), A(ITAPER),
868887
* A(IAREA), A(ISUMR), A(IISUM), LTMP)
869888
ELSE
870-
CALL CON3D(CRD, NDIM, NUMNP, LINK, NNODES, NUMEL, MAT,
889+
CALL CON3D(CRD, NDIM, NUMNP, LINK, 2**NDIM, NUMEL, MAT,
871890
* NELBLK, A(IELSEL), A(IASPEC), A(ISKEW), A(ITAPER),
872891
* A(IAREA), A(ISUMR), A(IISUM), A(ISKX), A(ISKY), A(ISKZ),
873892
* A(ITPX), A(ITPY), A(ITPZ), A(IJAC), LTMP)

‎packages/seacas/applications/numbers/nu_con2d.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,4 +1,4 @@
1-
C Copyright(C) 1999-2020 National Technology & Engineering Solutions
1+
C Copyright(C) 1999-2020, 2025 National Technology & Engineering Solutions
22
C of Sandia, LLC (NTESS). Under the terms of Contract DE-NA0003525 with
33
C NTESS, the U.S. Government retains certain rights in this software.
44
C
@@ -47,7 +47,7 @@ SUBROUTINE CON2D (CRD, NDIM, NUMNP, IX, NNODES, NUMEL, MAT,
4747
C=======================================================================
4848

4949
REAL CRD(NUMNP, NDIM), ASPECT(*), SKEW(*), TAPER(*), AREA(*)
50-
INTEGER IX(NNODES, NUMEL), MAT(6, NELBLK), ISUMRY(2,4,NELBLK)
50+
INTEGER IX(NNODES, NUMEL), MAT(7, NELBLK), ISUMRY(2,4,NELBLK)
5151
REAL SUMRY(4,4,NELBLK)
5252
LOGICAL SELECT(*), DEBUG, ISABRT
5353
include 'nu_io.blk'

0 commit comments

Comments
 (0)