44
45
46
47 USE intbufdef_mod
48
49
50
51#include "implicit_f.inc"
52
53
54
55#include "mvsiz_p.inc"
56
57
58
59#include "param_c.inc"
60
61
62
63 INTEGER (*), ISKEW(*), (*), ITAB(*)
64 my_real x(*), v(*), a(*), skew(*), fn(*), msmn(*), ft(*), msmt(*), ms(*) ,fsav(nthvki,*)
65
66 TYPE(INTBUF_STRUCT_) INTBUF_TAB
67
68
69
70 INTEGER NG, IGIMP, NGROUS
71 INTEGER, DIMENSION(MVSIZ) :: IX1,IX2,IX3,IX4
72 my_real,
DIMENSION(MVSIZ) :: x1,x2,x3,x4,xi
73 my_real,
DIMENSION(MVSIZ) :: y1,y2,y3,y4,yi
74 my_real,
DIMENSION(MVSIZ) :: z1,z2,z3,z4,zi
75 my_real,
DIMENSION(MVSIZ) :: xface,n1,n2,n3
77 my_real,
DIMENSION(MVSIZ) :: x0,y0,z0,ans
78 my_real,
DIMENSION(MVSIZ) :: xx1,xx2,xx3,xx4
79 my_real,
DIMENSION(MVSIZ) :: yy1,yy2,yy3,yy4
80 my_real,
DIMENSION(MVSIZ) :: zz1,zz2,zz3,zz4
81 my_real,
DIMENSION(MVSIZ) :: xi1,xi2,xi3,xi4
82 my_real,
DIMENSION(MVSIZ) :: yi1,yi2,yi3,yi4
83 my_real,
DIMENSION(MVSIZ) :: zi1,zi2,zi3,zi4
84 my_real,
DIMENSION(MVSIZ) :: xn1,xn2,xn3,xn4
85 my_real,
DIMENSION(MVSIZ) :: yn1,yn2,yn3,yn4
86 my_real,
DIMENSION(MVSIZ) :: zn1,zn2,zn3,zn4
87 my_real,
DIMENSION(MVSIZ) :: xp,yp,zp
88 my_real,
DIMENSION(MVSIZ) :: h1,h2,h3,h4
89 my_real startt, fric, gap, stopt,dist
90 INTEGER :: NMN,NSN,NTY
91 INTEGER :: LFT,LLT,NFT
92
93
94
95 nmn =ipari(6)
96 IF(nmn == 0)RETURN
97
98
99
100 nsn =ipari(5)
101 nty =ipari(7)
102
103
104
105 fric =intbuf_tab%VARIABLES(1)
106 gap =intbuf_tab%VARIABLES(2)
107 ngrous=1+(nsn-1)/nvsiz
108
109
110
111 DO 150 ng=1,ngrous
112 nft=(ng-1)*nvsiz
113 lft=1
114 llt=
min(nvsiz,nsn-nft)
116 1 x, intbuf_tab%IRECTM,intbuf_tab%LMSR, intbuf_tab%MSR,
117 2 intbuf_tab%NSV, intbuf_tab%ILOCS, intbuf_tab%NSEGM, xi,
118 3 yi, zi, xface, lft,
119 4 llt, nft)
121 1 x, intbuf_tab%IRECTM,intbuf_tab%LMSR, intbuf_tab%MSR,
122 2 intbuf_tab%NSV, intbuf_tab%ILOCS, intbuf_tab%IRTLM, intbuf_tab%NSEGM,
123 3 xface, lft, llt, nft)
125 1 x, v, a
126 2 intbuf_tab%MSR, intbuf_tab%NSV, intbuf_tab%IRTLM, ix1,
127 3 ix2, ix3, ix4, x1,
128 4 x2, x3, x4, xi,
129 5 y1, y2, y3, y4,
130 6 yi, z1, z2, z3,
131 7 z4, zi, lft, llt,
132 8 nft)
134 1 x1, x2, x3, x4,
135 2 xi, y1, y2, y3,
136 3 y4, yi, z1, z2,
137 4 z3, z4, zi, xface,
138 5 n1, n2, n3, ssc,
139 6 ttc, x0, y0, z0,
140 7 xx1, xx2, xx3, xx4,
141 8 yy1, yy2, yy3, yy4,
142 9 zz1, zz2, zz3, zz4,
143 a xi1, xi2, xi3, xi4,
144 b yi1, yi2, yi3, yi4,
145 c zi1, zi2, zi3, zi4,
146 d xn1, xn2, xn3, xn4,
147 e yn1, yn2, yn3, yn4,
148 f zn1, zn2, zn3, zn4,
151 1 gap,
area, thk, alp,
152 2 lft, llt)
154 1 igimp, nty, dist, x1,
155 2 x2, x3, x4, xi,
156 3 y1, y2, y3, y4,
157 4 yi, z1, z2, z3,
158 5 z4, zi, xface, n1,
159 6 n2, n3, ssc, ttc,
160 7 alp, ans, xp, yp,
161 8 zp, h1, h2, h3,
162 9 h4, lft, llt)
164 1 intbuf_tab%ILOCS,intbuf_tab%CSTS, intbuf_tab%N, xface,
165 2 n1, n2, n3, ssc,
166 3 ttc, thk, ans, lft,
167 4 llt, nft)
168 150 CONTINUE
169
170
171
173 1 (x ,skew ,a ,fn ,ft ,
174 2 intbuf_tab%IRECTM,intbuf_tab%LMSR,intbuf_tab%CSTS,intbuf_tab%MSR,
175 . intbuf_tab%NSV,
176 3 intbuf_tab%ILOCS,intbuf_tab%IRTLM,ms ,intbuf_tab%N,lcod ,
177 4 iskew ,fric ,msmn ,msmt , nmn, nsn)
178
179
180
182 1 (x ,skew ,a ,fn ,ft ,
183 2 intbuf_tab%IRECTM,intbuf_tab%LMSR,intbuf_tab%CSTS,intbuf_tab%MSR,
184 . intbuf_tab%NSV,
185 3 intbuf_tab%ILOCS,intbuf_tab%IRTLM,lcod ,v ,iskew ,
186 4 intbuf_tab%N,ms ,intbuf_tab%VARIABLES(1),msmn ,msmt ,
187 5 intbuf_tab%STFNS,itab ,fsav, nmn, nsn )
188
189 RETURN
subroutine area(d1, x, x2, y, y2, eint, stif0)
subroutine i3cst3(x1, x2, x3, x4, xi, y1, y2, y3, y4, yi, z1, z2, z3, z4, zi, xface, n1, n2, n3, ssc, ttc, x0, y0, z0, xx1, xx2, xx3, xx4, yy1, yy2, yy3, yy4, zz1, zz2, zz3, zz4, xi1, xi2, xi3, xi4, yi1, yi2, yi3, yi4, zi1, zi2, zi3, zi4, xn1, xn2, xn3, xn4, yn1, yn2, yn3, yn4, zn1, zn2, zn3, zn4, area, lft, llt)
subroutine i3dis3(igimp, nty, dist, x1, x2, x3, x4, xi, y1, y2, y3, y4, yi, z1, z2, z3, z4, zi, xface, n1, n2, n3, ssc, ttc, alp, ans, xp, yp, zp, h1, h2, h3, h4, lft, llt)
subroutine i3gap3(gap, area, thk, alp, lft, llt)
subroutine i3loc3(x, irect, lmsr, msr, nsv, iloc, nseg, xi, yi, zi, xface, lft, llt, nft)
subroutine i3msr3(x, irect, lmsr, msr, nsv, iloc, irtl, nseg, xface, lft, llt, nft)
subroutine i9avel(x, skew, a, fn, ft, irect, lmsr, crst, msr, nsv, iloc, irtl, lcode, v, iskew, nor, ms, frigap, msmn, msmt, efric, itab, fsav, nmn, nsn)
subroutine i9cor3(x, v, a, irect, msr, nsv, irtl, ix1, ix2, ix3, ix4, x1, x2, x3, x4, xi, y1, y2, y3, y4, yi, z1, z2, z3, z4, zi, lft, llt, nft)
subroutine i9frms(x, skew, a, fn, ft, irect, lmsr, crst, msr, nsv, iloc, irtl, ms, nor, lcode, iskew, fric, msmn, msmt, nmn, nsn)
subroutine i9imp3(iloc, cst, nor, xface, n1, n2, n3, ssc, ttc, thk, ans, lft, llt, nft)