32
33
34
35
36
37
38#include "implicit_f.inc"
39
40
41
42#include "mvsiz_p.inc"
43
44
45
46#include "vect01_c.inc"
47#include "param_c.inc"
48
49
50
52 . pm(npropm,*),geo(npropg,*), deltax(*), aire(*), vol(*), dtx(*)
53 INTEGER PID(*),MAT(*)
54
55
56
57 INTEGER I
58
60 . ssp(mvsiz), rho0(mvsiz), vis(mvsiz)
61
62
63
64
65
66 DO i=lft,llt
67 ssp(i)=pm(38,mat(i))
68 rho0(i) =pm( 1,mat(i))
69 vis(i) = zero
70 END DO
71
72 IF(jsph==0)THEN
73 CALL dtel(ssp,pm,geo,pid,mat, rho0, vis, deltax, aire, vol, dtx)
74 ELSE
75 CALL dtsph(ssp,pm,geo,pid,mat, rho0, vis, deltax, vol, dtx)
76 ENDIF
77
78 RETURN
subroutine dtel(ssp, pm, geo, pid, mat, rho0, vis, deltax, aire, vol, dtx)
subroutine dtsph(ssp, pm, geo, pid, mat, rho0, vis, deltax, vol, dtx)