#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
Go to the source code of this file.
|
| subroutine | inimu2 (pm, ix, f, df) |
◆ inimu2()
| subroutine inimu2 |
( |
| pm, |
|
|
integer, dimension(7,*) | ix, |
|
|
| f, |
|
|
| df ) |
Definition at line 28 of file inimu2.F.
29
30
31
32
33#include "implicit_f.inc"
34
35
36
37#include "com01_c.inc"
38#include "com04_c.inc"
39#include "param_c.inc"
40
41
42
43 INTEGER IX(7,*)
45 . pm(npropm,*), f(numnod,*), df(numnod,*)
46
47
48
49 INTEGER :: I, MAT, LAW, M, J, N
51
52
53 DO i=1,numelq
54 mat=iabs(ix(1,i))
55 law=int(pm(19,mat))
56 IF(law/=20)cycle
57 DO m=1,nmult
58 alph=pm(30+m,mat)
59 DO j=2,5
60 n=ix(j,i)
61 f(n,m)=f(n,m)+alph
62 df(n,m)=df(n,m)+one
63 END DO
64 END DO
65 END DO
66
67 RETURN