Commit 8a46059f authored by Stefan's avatar Stefan
Browse files

print correct dofs

parent 4dde21af
......@@ -1283,36 +1283,13 @@ c-----------------------------------------------------------------------
subroutine dofcnt
include 'SIZE'
include 'TOTAL'
COMMON /SCRNS/ WORK(LCTMP1)
integer*8 i8glsum
integer*8 ntot,ntotp,ntotv, nn
nxyz = lx1*ly1*lz1
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
integer*8 ntotv
ntot1=lx1*ly1*lz1*nelv
ntot2=lx2*ly2*lz2*nelv
ntotv = glsc2(tmult,tmask,ntot1)
nn = ntot2
ntotp = i8glsum(nn,1)
c if (ifflow) ntotv = glsc2(vmult,v1mask,ntot1)
c if (ifsplit) ntotp = glsc2(vmult,pmask ,ntot1)
c if (nio.eq.0) write(6,*) ' dofs:',ntotv,ntotp
if (ifflow) ntotv = glsc2(vmult,v1mask,ntot1)
if (nio.eq.0) write(6,*) 'dofs: ',ntotv
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