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#include "vect01_c.inc"
45
46
47
49 . pm(npropm,*),geo(npropg,*), deltax(*), aire(*), vol(*), dtx(*)
50 INTEGER PID(*),MAT(*)
51
52
53
54 INTEGER I, MX
55
57 . ssp(mvsiz), rho0(mvsiz), vis(mvsiz)
58
59
60
61 DO i=lft,llt
62 mx =mat(i)
63 ssp(i) =pm(49,mx)
64 rho0(i) =pm( 1,mx)
65 vis(i) =zero
66 END DO
67
68
69
70 IF(jsph==0)THEN
71 CALL dtel(ssp,pm,geo,pid,mat, rho0, vis, deltax, aire, vol, dtx)
72 ELSE
73 CALL dtsph(ssp,pm,geo,pid,mat, rho0, vis, deltax, vol, dtx)
74 ENDIF
75
76 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)