29
30
31
32#include "implicit_f.inc"
33
34
35
36#include "com01_c.inc"
37#include "flowcom.inc"
38
39
40
41
42 INTEGER IFUNC, IFLOW(*)
44 . rflow(*)
45 REAL WA4(*)
46
47
48
49 INTEGER IADI, IADR, I, ITYP, NINOUT, NNO, NEL, II1, II2,
50 . IR1, IR2, J, JJ, NNO_L, NNI_L, II3, II4, JJJ, NNI
51
52 iadi=0
53 iadr=0
54 DO i=1,nflow
55 ityp=iflow(iadi+2)
56 IF (ityp==1) THEN
57 ninout=iflow(iadi+4)
58 nno=iflow(iadi+5)
59 nel=iflow(iadi+6)
60 nni=iflow(iadi+7)
61 nno_l=iflow(iadi+16)
62 nni_l=iflow(iadi+22)
63 ii1=iadi+niflow
64 ii2=ii1+nno+3*nel+ninout*niioflow
65 IF(nspmd == 1) THEN
66 ii3=ii2+nni+nel
67 ii4=ii3+nno
68 ELSE
69 ii3=ii2+nni+nel+2*nno
70 ii4=ii3+2*nno
71 ENDIF
72 ir1=iadr+nrflow
73 ir2=ir1+nno+nni
74 IF (ifunc==3) THEN
75 DO j=1,nno_l
76 jj=iflow(ii3+j)
77 jjj=iflow(ii1+jj)
78 wa4(jjj)=rflow(ir2+jj)
79 ENDDO
80 DO j=1,nni_l
81 jj=iflow(ii4+j)
82 jjj=iflow(ii2+jj)
83 wa4(jjj)=rflow(ir2+nno+jj)
84 ENDDO
85 ELSEIF (ifunc==13) THEN
86 DO j=1,nno_l
87 jj=iflow(ii3+j)
88 jjj=iflow(ii1+jj)
89 wa4(jjj)=rflow(ir1+jj)
90 ENDDO
91 DO j=1,nni_l
92 jj=iflow(ii4+j)
93 jjj=iflow(ii2+jj)
94 wa4(jjj)=rflow(ir1+nno+jj)
95 ENDDO
96 ENDIF
97 ENDIF
98 iadr=iadr+iflow(iadi+15)
99 iadi=iadi+iflow(iadi+14)
100 ENDDO
101
102 RETURN