1c----------------------------------------------------------------------- 2c 3c module for time averaged statistics (conservative projection). 4c 5c----------------------------------------------------------------------- 6 module stats 7 8 integer nResDims, nSolDims, nLhsDims, nTimeStep, stsResFlg 9 integer stsCompFreq, stsWriteFreq, stsResetFreq, step1, 10 & stsType 11 12 real*8, allocatable :: stsVec(:,:) 13 14 real*8, allocatable :: stsReg(:) 15 real*8, allocatable :: stsMInv(:,:) 16 real*8, allocatable :: stsB(:,:) 17 real*8, allocatable :: stsDInv(:,:) 18 real*8, allocatable :: stsCInv(:,:) 19 20 real*8, allocatable :: stsPres(:), stsPresSqr(:), stsVel(:,:), 21 & stsVelSqr(:,:), stsVelReg(:,:), 22 & stsStress(:,:) 23 24 end module 25 26c----------------------------------------------------------------------- 27c create the new statistics arrays 28c----------------------------------------------------------------------- 29 subroutine initStats(x, iBC, iper, ilwork) 30 31 use stats 32 include "common.h" 33 34 real*8 x(numnp,3) 35 integer ilwork(nlwork), iper(nshg), iBC(nshg) 36 37 write(*,*) 'Stats are not developed for compressible code' 38 return 39 end 40