Test writing VTK files
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.
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
!! Test writing VTK files
real(wp),dimension(:),allocatable::x,y
real(wp),dimension(:,:),allocatable::XX,YY
real(wp),dimension(:,:),allocatable::F,U,V
integer::N,M,iou
N = 100
M = 101
x = linspace(0.0_wp,1.0_wp,N)
y = linspace(0.0_wp,1.0_wp,M)
XX = MeshGridX(x,y)
YY = MeshGridY(x,y)
F = XX*YY
U = -YY
V = XX
open(file='data.vtk',newunit=iou)
call writeHeaderVTK(iou,'Test File')
call writeGridVTK(iou,XX,YY)
call writeScalarVTK(iou,'F',F)
call writeVectorVTK(iou,'U',U,V)
close(iou)
end subroutine testWriteVTK