1*53c9b1fcSKenneth E. Jansen module solvedata 2*53c9b1fcSKenneth E. Jansen 3*53c9b1fcSKenneth E. Jansen integer nsolflow,npermDims, nTmpDims, nPermDimsS, nTmpDimsS 4*53c9b1fcSKenneth E. Jansen 5*53c9b1fcSKenneth E. Jansen real*8, allocatable :: lhsP(:,:), lhsK(:,:), lhsS(:,:) 6*53c9b1fcSKenneth E. Jansen real*8, allocatable :: aperm(:,:), atemp(:,:) 7*53c9b1fcSKenneth E. Jansen real*8, allocatable :: apermS(:,:,:), atempS(:,:) 8*53c9b1fcSKenneth E. Jansen 9*53c9b1fcSKenneth E. Jansen end module 10*53c9b1fcSKenneth E. Jansen 11*53c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 12*53c9b1fcSKenneth E. Jansenc allocate the arrays 13*53c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 14*53c9b1fcSKenneth E. Jansen 15*53c9b1fcSKenneth E. Jansen 16*53c9b1fcSKenneth E. Jansen subroutine aSD 17*53c9b1fcSKenneth E. Jansen 18*53c9b1fcSKenneth E. Jansen use solvedata 19*53c9b1fcSKenneth E. Jansen include "common.h" 20*53c9b1fcSKenneth E. Jansen if(nsolflow.eq.1) then 21*53c9b1fcSKenneth E. Jansen allocate (lhsP(4,nnz_tot)) 22*53c9b1fcSKenneth E. Jansen allocate (lhsK(9,nnz_tot)) 23*53c9b1fcSKenneth E. Jansen if(leslib.eq.1) then 24*53c9b1fcSKenneth E. Jansen allocate (aperm(nshg,nPermDims)) 25*53c9b1fcSKenneth E. Jansen allocate (atemp(nshg,nTmpDims)) 26*53c9b1fcSKenneth E. Jansen endif 27*53c9b1fcSKenneth E. Jansen endif 28*53c9b1fcSKenneth E. Jansen 29*53c9b1fcSKenneth E. Jansen return 30*53c9b1fcSKenneth E. Jansen end 31*53c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 32*53c9b1fcSKenneth E. Jansenc delete the arrays 33*53c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 34*53c9b1fcSKenneth E. Jansen 35*53c9b1fcSKenneth E. Jansen 36*53c9b1fcSKenneth E. Jansen subroutine dSD 37*53c9b1fcSKenneth E. Jansen 38*53c9b1fcSKenneth E. Jansen use solvedata 39*53c9b1fcSKenneth E. Jansen 40*53c9b1fcSKenneth E. Jansen deallocate (lhsP) 41*53c9b1fcSKenneth E. Jansen 42*53c9b1fcSKenneth E. Jansen return 43*53c9b1fcSKenneth E. Jansen end 44