Commit 98e96957 authored by Stefan's avatar Stefan
Browse files

print dof correctly

parent e2e94af5
......@@ -1285,13 +1285,35 @@ c-----------------------------------------------------------------------
subroutine dofcnt
include 'SIZE'
include 'TOTAL'
COMMON /SCRNS/ WORK(LCTMP1)
integer*8 ntotv
integer*8 ntot,ntotp,ntotv
nxyz = nx1*ny1*nz1
nel = nelv
! unique points on v-mesh
vpts = glsum(vmult,nel*nxyz) + .1
nvtot=vpts
! unique points on pressure mesh
work(1)=nel*nxyz
ppts = glsum(work,1) + .1
ntot=ppts
C
if (nio.eq.0) write(6,'(A,2i13)')
& 'gridpoints unique/tot: ',nvtot,ntot
ntot1=nx1*ny1*nz1*nelv
ntot2=nx2*ny2*nz2*nelv
ntot1=lx1*ly1*lz1*nelv
ntotv = glsc2(tmult,tmask,ntot1)
if (ifflow) ntotv = glsc2(vmult,v1mask,ntot1)
if (nio.eq.0) write(6,*) 'dofs: ',ntotv
ntotp = i8glsum(ntot2,1)
if (ifflow) ntotv = glsc2(vmult,v1mask,ntot1)
if (ifsplit) ntotp = glsc2(vmult,pmask ,ntot1)
if (nio.eq.0) write(6,'(A,2i13)')
$ 'dofs vel/pr: ',ntotv,ntotp
return
end
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment