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