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
59 . ssp(mvsiz), rho0(mvsiz), vis(mvsiz)
60
61
62
63
64
65 DO i=lft,llt
66 ssp(i)=pm(38,mat(i))
67 rho0(i) =pm( 1,mat(i))
68 vis(i) = zero
69 END DO
70
71 IF(jsph==0)THEN
72 CALL dtel(ssp,pm,geo,pid,mat, rho0, vis, deltax, aire, vol, dtx)
73 ELSE
74 CALL dtsph(ssp,pm,geo,pid,mat, rho0, vis, deltax, vol, dtx)
75 ENDIF
76
77 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)