35
36
37
39
40
41
42#include "implicit_f.inc"
43#include "comlock.inc"
44#include "task_c.inc"
45
46
47
48#include "param_c.inc"
49
50
51
52 INTEGER, INTENT(IN) :: NEL
53 INTEGER, INTENT(IN) :: NFT
54 INTEGER, INTENT(IN) :: JALE
55 INTEGER, INTENT(IN) :: JEUL
56 INTEGER, INTENT(IN) :: NUMMAT
57 INTEGER, INTENT(IN) ::
58 my_real pm(npropm,nummat),flux(6,*),bufmat(*)
59 INTEGER, INTENT(IN) :: IXS(NIXS,S_IXS),IPM(NPROPMI,NUMMAT)
60 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
61
62
63
64 INTEGER I, J,JJ,II,K,IFLG,IADBUF,IAD2
65
66
67
68 IF(jeul+jale==1)THEN
69 DO i=1,nel
70 iad2 = ale_connect%ee_connect%iad_connect(i + nft)
71 DO k=1,6
72 jj = ale_connect%ee_connect%connected(iad2 + k - 1)
73 IF(jj>0)THEN
74 j = ixs(1,jj)
75 IF(nint(pm(50,j))==8)THEN
76 !test law11
77 IF(nint(pm(19,j))==11)THEN
78 flux(k,i) =
max(flux(k,i),zero)
79
80 ELSEIF(nint(pm(19,j))==51)THEN
81 iadbuf = ipm(27,j)
82 iflg = nint(bufmat(31+iadbuf-1))
83 IF (iflg == 6)THEN
84 flux(k,i) =
max(flux(k,i),zero)
85 ENDIF
86 endif!law test
87 endif
88 endif
89 ENDDO
90 ENDDO
91 ENDIF
92
93 RETURN