153c9b1fcSKenneth E. Jansen module solvedata 253c9b1fcSKenneth E. Jansen 353c9b1fcSKenneth E. Jansen integer nsolflow,npermDims, nTmpDims, nPermDimsS, nTmpDimsS 453c9b1fcSKenneth E. Jansen 553c9b1fcSKenneth E. Jansen real*8, allocatable :: lhsP(:,:), lhsK(:,:), lhsS(:,:) 653c9b1fcSKenneth E. Jansen real*8, allocatable :: aperm(:,:), atemp(:,:) 753c9b1fcSKenneth E. Jansen real*8, allocatable :: apermS(:,:,:), atempS(:,:) 853c9b1fcSKenneth E. Jansen 953c9b1fcSKenneth E. Jansen end module 1053c9b1fcSKenneth E. Jansen 1153c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 1253c9b1fcSKenneth E. Jansenc allocate the arrays 1353c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 14*da6af917SKenneth E. Jansen subroutine aSDf 1553c9b1fcSKenneth E. Jansen use solvedata 1653c9b1fcSKenneth E. Jansen include "common.h" 1753c9b1fcSKenneth E. Jansen allocate (lhsP(4,nnz_tot)) 1853c9b1fcSKenneth E. Jansen allocate (lhsK(9,nnz_tot)) 1953c9b1fcSKenneth E. Jansen if(leslib.eq.1) then 2053c9b1fcSKenneth E. Jansen allocate (aperm(nshg,nPermDims)) 2153c9b1fcSKenneth E. Jansen allocate (atemp(nshg,nTmpDims)) 2253c9b1fcSKenneth E. Jansen endif 23*da6af917SKenneth E. Jansen return 24*da6af917SKenneth E. Jansen end 2553c9b1fcSKenneth E. Jansen 26*da6af917SKenneth E. Jansen subroutine aSDs(nsclrsol) 27*da6af917SKenneth E. Jansen use solvedata 28*da6af917SKenneth E. Jansen include "common.h" 29*da6af917SKenneth E. Jansen allocate (lhsS(nnz_tot,nsclrsol)) 30*da6af917SKenneth E. Jansen if(leslib.eq.1) then 31*da6af917SKenneth E. Jansen allocate (apermS(nshg,nPermDimsS,nsclrsol)) 32*da6af917SKenneth E. Jansen allocate (atempS(nshg,nTmpDimsS)) 33*da6af917SKenneth E. Jansen endif 3453c9b1fcSKenneth E. Jansen return 3553c9b1fcSKenneth E. Jansen end 3653c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 3753c9b1fcSKenneth E. Jansenc delete the arrays 3853c9b1fcSKenneth E. Jansenc----------------------------------------------------------------------- 39*da6af917SKenneth E. Jansen subroutine dSDf 4053c9b1fcSKenneth E. Jansen use solvedata 41*da6af917SKenneth E. Jansen include "common.h" 4253c9b1fcSKenneth E. Jansen deallocate (lhsP) 43*da6af917SKenneth E. Jansen deallocate (lhsK) 44*da6af917SKenneth E. Jansen if(leslib.eq.1) then 45*da6af917SKenneth E. Jansen deallocate (aperm) 46*da6af917SKenneth E. Jansen deallocate (atemp) 47*da6af917SKenneth E. Jansen endif 48*da6af917SKenneth E. Jansen return 49*da6af917SKenneth E. Jansen end 5053c9b1fcSKenneth E. Jansen 51*da6af917SKenneth E. Jansen subroutine dSDs 52*da6af917SKenneth E. Jansen use solvedata 53*da6af917SKenneth E. Jansen include "common.h" 54*da6af917SKenneth E. Jansen deallocate (lhsS) 55*da6af917SKenneth E. Jansen if(leslib.eq.1) then 56*da6af917SKenneth E. Jansen deallocate (apermS) 57*da6af917SKenneth E. Jansen deallocate (atempS) 58*da6af917SKenneth E. Jansen endif 5953c9b1fcSKenneth E. Jansen return 6053c9b1fcSKenneth E. Jansen end 61