program vtk_fortran_write_vtu
!< VTK_Fortran test: write VTU file.
use penf
use vtk_fortran, only : vtk_file
implicit none
type(vtk_file) :: a_vtk_file !< A VTK file.
integer(I4P), parameter :: np = 27_I4P !< Number of points.
integer(I4P), parameter :: nc = 11_I4P !< Number of cells.
real(R4P), dimension(1:np) :: x = [0,1,2,0,1,2,0,1,2,0,1,2,0,1,2,0,1,2,0,1,2,0,1,2,0,1,2] !< X coordinates of points.
real(R4P), dimension(1:np) :: y = [0,0,0,1,1,1,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1] !< Y coordinates of points.
real(R4P), dimension(1:np) :: z = [0,0,0,0,0,0,1,1,1,1,1,1,2,2,2,3,3,3,4,4,4,5,5,5,6,6,6] !< Z coordinates of points.
integer(I1P), dimension(1:nc) :: cell_type = [12_I1P,12_I1P,10_I1P,10_I1P,&
7_I1P,6_I1P,9_I1P,5_I1P,5_I1P,3_I1P,1_I1P] !< Cells type.
integer(I4P), dimension(1:nc) :: offset = [8_I4P,16_I4P,20_I4P,24_I4P,30_I4P,&
36_I4P,40_I4P,43_I4P,46_I4P,48_I4P,49_I4P] !< Cells offset.
integer(I4P), dimension(1:49) :: connect !< Connectivity.
real(R8P), dimension(1:np) :: v !< One component points-variable.
integer(I4P), dimension(1:np) :: v_u !< X component points-variable.
integer(I4P), dimension(1:np) :: v_v !< Y component points-variable.
integer(I4P), dimension(1:np) :: v_w !< Z component points-variable.
integer(I4P), dimension(1:np) :: v_x !< X component points-variable.
integer(I4P), dimension(1:np) :: v_y !< Y component points-variable.
integer(I4P), dimension(1:np) :: v_z !< Z component points-variable.
integer(I4P) :: error !< Status error.
logical :: test_passed(1) !< List of passed tests.
connect = [0 ,1 ,4 ,3 ,6 ,7 ,10,9 , &
1 ,2 ,5 ,4 ,7 ,8 ,11,10, &
6 ,10,9 ,12, &
5 ,11,10,14, &
15,16,17,14,13,12, &
18,15,19,16,20,17, &
22,23,20,19, &
21,22,18, &
22,19,18, &
26,25, &
24]
v = [0.0,1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,17.0, &
18.0,19.0,20.0,21.0,22.0,23.0,24.0,25.0,26.0]
v_U = [1,1,0,1,1,0,1,1,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
v_V = [0,1,2,0,1,2,0,1,2,0,1,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
v_W = [0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]
v_X = [1,1,0,1,1,0,1,1,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
v_Y = [0,1,2,0,1,2,0,1,2,0,1,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
v_Z = [0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]
! ascii
error = a_vtk_file%initialize(format='ascii', filename='XML_UNST_T-ascii.vtu', mesh_topology='UnstructuredGrid')
call write_data
error = a_vtk_file%finalize()
! raw
error = a_vtk_file%initialize(format='raw', filename='XML_UNST_T-raw.vtu', mesh_topology='UnstructuredGrid')
call write_data
error = a_vtk_file%finalize()
! binary
error = a_vtk_file%initialize(format='binary', filename='XML_UNST_T-binary.vtu', mesh_topology='UnstructuredGrid')
call write_data
error = a_vtk_file%finalize()
test_passed = .true. ! nothing to test yet
print "(A,L1)", new_line('a')//'Are all tests passed? ', all(test_passed)
stop
contains
subroutine write_data
!< Write data.
error = a_vtk_file%xml_writer%write_piece(np=np, nc=nc)
error = a_vtk_file%xml_writer%write_geo(np=np, nc=nc, x=x, y=y, z=z)
error = a_vtk_file%xml_writer%write_connectivity(nc=nc, connectivity=connect, offset=offset, cell_type=cell_type)
error = a_vtk_file%xml_writer%write_dataarray(location='node', action='open')
error = a_vtk_file%xml_writer%write_dataarray(data_name='scalars', x=v)
error = a_vtk_file%xml_writer%write_dataarray(data_name='tensor', u=v_u, v= v_v, w=v_w, x=v_x, y=v_y, z=v_z)
error = a_vtk_file%xml_writer%write_dataarray(location='node', action='close')
error = a_vtk_file%xml_writer%write_piece()
endsubroutine write_data
endprogram vtk_fortran_write_vtu