1         subroutine vortGLB  (y, x, shp, shgl, ilwork, vortG)  argument
7         real*8  y(nshg,ndof),  x(numnp,nsd)                local
48         call vortELE(y,x,tmpshp,tmpshgl,mien(iblk)%p,vortIntG,lpmassG)
78       subroutine vortELE(y,x,shp1,shgl1,ien,vortIntG,lpmassG)  argument
83        real*8 y(nshg,ndof), x(numnp,nsd)  local
108        call localy (y, yl, ien, ndof,  'gather  ')
192        subroutine PRcalcDuct  (y)  argument
197        real*8 y(nshg,ndof)    local
207        Mach(:)=sqrt(y(:,1)**2+y(:,2)**2+y(:,3)**2)/sqrt(1.4*287*y(:,5)) 
208        PresRec(:)=y(:,4)*(1+0.2*Mach(:)**2)**3.5/Ptotal_contra 
209        PresCoeff(:)=(y(:,4)-Pstatic_throat)/(Ptotal_contra-Pstatic_throat) 
210        velocity(:,1:3)=y(:,1:3)