xref: /phasta/phSolver/incompressible/solvedata.f (revision 53c9b1fc264f732de8bfb0deaf531857e2919987)
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