#include "implicit_f.inc"
#include "mvsiz_p.inc"
Go to the source code of this file.
|
| subroutine | sesa30 (sig, eint, qold, vol, qnew, vnew, dv2, nel, p01) |
◆ sesa30()
| subroutine sesa30 |
( |
| sig, |
|
|
| eint, |
|
|
| qold, |
|
|
| vol, |
|
|
| qnew, |
|
|
| vnew, |
|
|
| dv2, |
|
|
integer | nel, |
|
|
| p01 ) |
Definition at line 28 of file sesa30.F.
30
31
32
33#include "implicit_f.inc"
34
35
36
37#include "mvsiz_p.inc"
38
39
40
41
42
43
44 INTEGER NEL
46 . sig(nel,6), eint(*), qold(*), vol(*) ,qnew(*), vnew(*), dv2(*), p01(*)
47
48
49
50 INTEGER I
52 . p(mvsiz)
53
54 DO i=1,nel
55 p(i)=p01(i)
56 eint(i)=
max(zero,eint(i)-dv2(i)*(p(i)+qnew(i)))
57 eint(i)=eint(i)/
max(em30,vol(i))
58 ENDDO
59
60 DO i=1,nel
61 sig(i,1)=sig(i,1)-p(i)
62 sig(i,2)=sig(i,2)-p(i)
63 sig(i,3)=sig(i,3)-p(i)
64 ENDDO
65
66 RETURN