From f8507cfefd4a0940b08ee5c51be727681fdf1115 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Fri, 3 Apr 2009 09:46:10 +0000 Subject: [PATCH] suppressed output of subnode coordinates to *.out file, but still available in 'mesh_tell_statistics' --- trunk/mesh.f90 | 27 +++++++++++++++++++-------- 1 file changed, 19 insertions(+), 8 deletions(-) diff --git a/trunk/mesh.f90 b/trunk/mesh.f90 index 69f340ab9..f8112ea9c 100644 --- a/trunk/mesh.f90 +++ b/trunk/mesh.f90 @@ -2837,7 +2837,6 @@ SUBROUTINE mesh_get_nodeElemDimensions (unit) real(pReal), dimension(3) :: centerOfGravity allocate(mesh_ipVolume(mesh_maxNips,mesh_NcpElems)) ; mesh_ipVolume = 0.0_pReal - !write(6,'(a10,x,a20,x,a20,x,a25,3(x,a6))') 'FE_Nips','FE_NipNeighbors','FE_NipFaceNodes','FE_subNodeOnIPFace','x','y','z' do e = 1,mesh_NcpElems ! loop over cpElems t = mesh_element(2,e) ! get elemType do i = 1,FE_Nips(t) ! loop over IPs of elem @@ -2847,10 +2846,6 @@ SUBROUTINE mesh_get_nodeElemDimensions (unit) do n = 1,FE_NipFaceNodes ! loop over nodes on interface gravityNode(FE_subNodeOnIPFace(n,f,i,t)) = 1 gravityNodePos(:,FE_subNodeOnIPFace(n,f,i,t)) = mesh_subNodeCoord(:,FE_subNodeOnIPFace(n,f,i,t),e) - write(6,'(i10,x,i20,x,i20,x,i25,3(x,f6.3))') i,f,n,FE_subNodeOnIPFace(n,f,i,t),& - mesh_subNodeCoord(1,FE_subNodeOnIPFace(n,f,i,t),e),& - mesh_subNodeCoord(2,FE_subNodeOnIPFace(n,f,i,t),e),& - mesh_subNodeCoord(3,FE_subNodeOnIPFace(n,f,i,t),e) end do end do @@ -2978,9 +2973,25 @@ SUBROUTINE mesh_get_nodeElemDimensions (unit) write (6,"(i5,x,i5,x,e15.8)") e,i,mesh_IPvolume(i,e) do f = 1,FE_NipNeighbors(mesh_element(2,e)) write (6,"(i33,x,e15.8,x,3(f6.3,x))") f,mesh_ipArea(f,i,e),mesh_ipAreaNormal(:,f,i,e) - enddo - enddo - enddo + end do + end do + end do + write (6,*) + !write(6,'(a5,x,a5,x,a15,x,a15,x,a20,3(x,a8))') 'elem','IP','IP neighbor','IPFaceNodes','subNodeOnIPFace','x','y','z' + !do e = 1,mesh_NcpElems ! loop over cpElems + ! t = mesh_element(2,e) ! get elemType + ! do i = 1,FE_Nips(t) ! loop over IPs of elem + ! do f = 1,FE_NipNeighbors(t) ! loop over interfaces of IP + ! do n = 1,FE_NipFaceNodes ! loop over nodes on interface + ! write(6,'(i5,x,i5,x,i15,x,i15,x,i20,3(x,f8.3))') e,i,f,n,FE_subNodeOnIPFace(n,f,i,t),& + ! mesh_subNodeCoord(1,FE_subNodeOnIPFace(n,f,i,t),e),& + ! mesh_subNodeCoord(2,FE_subNodeOnIPFace(n,f,i,t),e),& + ! mesh_subNodeCoord(3,FE_subNodeOnIPFace(n,f,i,t),e) + ! end do + ! end do + ! end do + !end do + write (6,*) write (6,*) write (6,*) "Input Parser: STATISTICS" write (6,*)