#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
Go to the source code of this file.
|
| subroutine | s8lagr3 (pm, vol0, rho, eint, mxt, voln, rho0, dvol, vd2, nel) |
◆ s8lagr3()
| subroutine s8lagr3 |
( |
| pm, |
|
|
| vol0, |
|
|
| rho, |
|
|
| eint, |
|
|
integer, dimension(mvsiz) | mxt, |
|
|
| voln, |
|
|
| rho0, |
|
|
| dvol, |
|
|
| vd2, |
|
|
integer, intent(in) | nel ) |
Definition at line 28 of file s8lagr3.F.
32
33
34
35#include "implicit_f.inc"
36
37
38
39#include "mvsiz_p.inc"
40
41
42
43#include "param_c.inc"
44
45
46
47 INTEGER, INTENT(IN) :: NEL
48 my_real pm(npropm,*), vol0(*), rho(*), eint(*)
49 my_real voln(mvsiz), rho0(mvsiz), dvol(mvsiz), vd2(mvsiz)
50 INTEGER MXT(MVSIZ)
51
52 INTEGER I, MX
53
54 mx = mxt(1)
55 DO i=1,nel
56 rho0(i)=pm(1,mx)
57 ENDDO
58
59 DO i=1,nel
60 dvol(i)=voln(i)-(rho0(1)/rho(i))*vol0(i)
61 ENDDO
62
63 DO i=1,nel
64 vd2(i)=zero
65 rho(i)=rho0(1)*(vol0(i)/voln(i))
66 eint(i)=eint(i)*vol0(i)
67 ENDDO
68
69 RETURN