Verify operation of readGmsh
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed arrows point from an interface to procedures which implement that interface. This could include the module procedures in a generic interface or the implementation in a submodule of an interface in a parent module.
subroutine testWriteVTK
!! Verify operation of readGmsh
type(mesh_t)::m
real(wp),dimension(:),allocatable::s
real(wp),dimension(:,:),allocatable::v
integer::N,k
call execute_command_line('gmsh -2 ./input/square.geo -o square.msh')
call m%readGmsh('square.msh')
N = size(m%nodes)
allocate( s(N) , v(N,3) )
v(:,3) = 0.0_wp
do k=1,N
s(k) = norm2( m%nodes(k)%x )
v(k,1:2) = m%nodes(k)%x
end do
call m%writeVTK('square.vtk')
call m%appendScalarVTK('square.vtk',s,'s')
call m%appendVectorVTK('square.vtk',v,'v')
end subroutine testWriteVTK