Mesh Nodes Numbering and xml
Hello all,
I wrote a post-processing subroutine (Finite Element Method) in HDf5 format. I wrote also xml file to read the data in Paraview:
if (rank.eq.0) then
print*,'| ...WRITING: '// ou //'/sol' &
//trim(adjustl(fnum))//'.xmf'
print*,'|'
mesh_filename= 'Mesh.h5'
sol_filename= 'sol'//trim(adjustl(fnum))//'.h5'
GridName = 'sol'//trim(adjustl(fnum))
xml_filename = 'paraview'// ou //'/sol'//trim(adjustl(fnum))//'.xmf'
i1 = index(GridName,' ')-1
ixmf=987
open(ixmf,file=xml_filename)
write(ixmf,10) GridName(1:i1)
write(ixmf,40) num
write(ixmf,20) neles,neles,9, trim(mesh_filename)//':/Elements'
write(ixmf,30) nodes,3, trim(mesh_filename)//':/Nodes'
write(ixmf,50) 'Phase-field', 'Scalar', 'Node', nodes, trim(sol_filename)//':/Mu'
write(ixmf,50) 'Pressure', 'Scalar', 'Node', nodes, trim(sol_filename)//':/P'
write(ixmf,50) 'Displacement', 'Vector', 'Node', nodes*3, trim(sol_filename)//':/U'
write(ixmf,50) 'Velocity', 'Vector', 'Node', nodes*3, trim(sol_filename)//':/V'
write(ixmf,50) 'GradMu', 'Vector', 'Node', nodes*3, trim(sol_filename)//':/GradMu'
write(ixmf,50) 'Stress', 'Tensor', 'Node', 9*nodes, trim(sol_filename)//':/Stress'
write(ixmf,50) 'Strain', 'Tensor', 'Node', 9*nodes, trim(sol_filename)//':/Strain'
write(ixmf,60)
close(ixmf)
10 format( '',/,& '',/,& '',/,& ' ',/,& ' ')
20 format( ' ',/,& ' ', A, '',/,& ' ')
30 format( ' ',/,& ' ', A, '',/,& ' ')
40 format( ' ')
50 format( ' ',/,& ' ', A, '',/,& ' ')
60 format( ' ',/,& ' ',/,& '')
When I display the data in Paraview, I got this:
I think I have issue either the nodes numbering in each element or the xml format. Could anyone help?
Thank you, Talal
Edited by Talal Alotaibi