Skip to content

Commit

Permalink
Update export to vtk subroutines
Browse files Browse the repository at this point in the history
  • Loading branch information
@gha3mi
gha3mi committed Apr 1, 2024
1 parent a97525b commit 360fd5d
Show file tree
Hide file tree
Showing 6 changed files with 100 additions and 64 deletions.
26 changes: 16 additions & 10 deletions src/NURBS/forcad_nurbs_curve.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -327,17 +327,20 @@ impure subroutine export_Xc(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xc'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
if (size(this%Xc, 2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), 0.0_rk , i = 1, nc)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*3
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(2, elemConn(i,1),elemConn(i,2), i = 1, size(elemConn,1))
(2, elemConn(i,1)-1,elemConn(i,2)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (3 , i = 1, size(elemConn,1))
Expand Down Expand Up @@ -366,17 +369,20 @@ impure subroutine export_Xg(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xg'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
if (size(this%Xg, 2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), 0.0_rk , i = 1, ng)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*3
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(2, elemConn(i,1),elemConn(i,2), i = 1, size(elemConn,1))
(2, elemConn(i,1)-1,elemConn(i,2)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (3 , i = 1, size(elemConn,1))
Expand Down
26 changes: 16 additions & 10 deletions src/NURBS/forcad_nurbs_surface.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -373,17 +373,20 @@ impure subroutine export_Xc(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xc'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
if (size(this%Xc,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), 0.0_rk , i = 1, nc)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*5
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(4, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3), i = 1, size(elemConn,1))
(4, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (9 , i = 1, size(elemConn,1))
Expand Down Expand Up @@ -412,17 +415,20 @@ impure subroutine export_Xg(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xg'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
if (size(this%Xg,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), 0.0_rk , i = 1, ng)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*5
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(4, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3), i = 1, size(elemConn,1))
(4, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (9 , i = 1, size(elemConn,1))
Expand Down
30 changes: 18 additions & 12 deletions src/NURBS/forcad_nurbs_volume.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -408,18 +408,21 @@ impure subroutine export_Xc(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xc'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
if (size(this%Xc,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), 0.0_rk , i = 1, nc)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*9
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(8, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3),&
elemConn(i,5),elemConn(i,6),elemConn(i,8),elemConn(i,7), i = 1, size(elemConn,1))
(8, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1,&
elemConn(i,5)-1,elemConn(i,6)-1,elemConn(i,8)-1,elemConn(i,7)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (12 , i = 1, size(elemConn,1))
Expand Down Expand Up @@ -448,18 +451,21 @@ impure subroutine export_Xg(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xg'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
if (size(this%Xg,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), 0.0_rk , i = 1, ng)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*9
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(8, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3),&
elemConn(i,5),elemConn(i,6),elemConn(i,8),elemConn(i,7), i = 1, size(elemConn,1))
(8, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1,&
elemConn(i,5)-1,elemConn(i,6)-1,elemConn(i,8)-1,elemConn(i,7)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (12 , i = 1, size(elemConn,1))
Expand Down
26 changes: 16 additions & 10 deletions src/Rational_Bezier/forcad_bezier_curve.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -289,17 +289,20 @@ impure subroutine export_Xc(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xc'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
if (size(this%Xc,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), 0.0_rk , i = 1, nc)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*3
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(2, elemConn(i,1),elemConn(i,2), i = 1, size(elemConn,1))
(2, elemConn(i,1)-1,elemConn(i,2)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (3 , i = 1, size(elemConn,1))
Expand Down Expand Up @@ -329,17 +332,20 @@ impure subroutine export_Xg(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xg'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
if (size(this%Xg,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), 0.0_rk , i = 1, ng)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*3
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(2, elemConn(i,1),elemConn(i,2), i = 1, size(elemConn,1))
(2, elemConn(i,1)-1,elemConn(i,2)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (3 , i = 1, size(elemConn,1))
Expand Down
26 changes: 16 additions & 10 deletions src/Rational_Bezier/forcad_bezier_surface.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -312,17 +312,20 @@ impure subroutine export_Xc(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xc'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
if (size(this%Xc,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), 0.0_rk , i = 1, nc)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*5
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(4, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3), i = 1, size(elemConn,1))
(4, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (9 , i = 1, size(elemConn,1))
Expand Down Expand Up @@ -351,17 +354,20 @@ impure subroutine export_Xg(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xg'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
if (size(this%Xg,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), 0.0_rk , i = 1, ng)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*5
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(4, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3), i = 1, size(elemConn,1))
(4, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (9 , i = 1, size(elemConn,1))
Expand Down
30 changes: 18 additions & 12 deletions src/Rational_Bezier/forcad_bezier_volume.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -332,18 +332,21 @@ impure subroutine export_Xc(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xc'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', nc, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
if (size(this%Xc,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), 0.0_rk , i = 1, nc)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xc(i,1), this%Xc(i,2), this%Xc(i,3) , i = 1, nc)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*9
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(8, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3),&
elemConn(i,5),elemConn(i,6),elemConn(i,8),elemConn(i,7), i = 1, size(elemConn,1))
(8, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1,&
elemConn(i,5)-1,elemConn(i,6)-1,elemConn(i,8)-1,elemConn(i,7)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (12 , i = 1, size(elemConn,1))
Expand Down Expand Up @@ -372,18 +375,21 @@ impure subroutine export_Xg(this, filename)

open(newunit=nunit, file=filename, action='write')
write(nunit,'(a)') '# vtk DataFile Version 2.0'
write(nunit,'(a)') 'Xg'
write(nunit,'(a)') 'Generated by ForCAD'
write(nunit,'(a)') 'ASCII'
write(nunit,'(a)') 'DATASET UNSTRUCTURED_GRID'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng+1, 'double'
write(nunit,'(a," ",g0," ",a)') 'POINTS', ng, 'double'

write(nunit,'(f24.18,f24.18,f24.18)') 0.0_rk, 0.0_rk, 0.0_rk
write(nunit,'(f24.18,f24.18,f24.18)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
if (size(this%Xg,2) == 2) then
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), 0.0_rk , i = 1, ng)
else
write(nunit,'(g0," ",g0," ",g0)') (this%Xg(i,1), this%Xg(i,2), this%Xg(i,3) , i = 1, ng)
end if

write(nunit,'(a," ",g0," ",g0)') 'CELLS', size(elemConn,1), size(elemConn,1)*9
write(nunit,'(g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0," ",g0)')&
(8, elemConn(i,1),elemConn(i,2),elemConn(i,4),elemConn(i,3),&
elemConn(i,5),elemConn(i,6),elemConn(i,8),elemConn(i,7), i = 1, size(elemConn,1))
(8, elemConn(i,1)-1,elemConn(i,2)-1,elemConn(i,4)-1,elemConn(i,3)-1,&
elemConn(i,5)-1,elemConn(i,6)-1,elemConn(i,8)-1,elemConn(i,7)-1, i = 1, size(elemConn,1))

write(nunit,'(a," ",g0)') 'CELL_TYPES', size(elemConn,1)
write(nunit,'(g0)') (12 , i = 1, size(elemConn,1))
Expand Down

0 comments on commit 360fd5d

Please sign in to comment.