59
60
61
63 USE sensor_mod
64 use element_mod , only : nixr
65
66
67
68#include "implicit_f.inc"
69
70
71
72#include "com04_c.inc"
73#include "mvsiz_p.inc"
74#include "parit_c.inc"
75
76
77
78 INTEGER, INTENT(IN) :: NFT,NSENSOR
79 INTEGER, INTENT(IN) :: IGRE
80 INTEGER, INTENT(IN) :: JSMS
81 INTEGER JFT,JLT,IOUT,NEL,NUVAR,IPROP,
82 . IPARTR(*),IXR(NIXR,*),IADR(3,*),GRTH(*),IGRTH(*),
83 . NELTST,ITYPTST,IGTYP,NC1(*),NC2(*)
85 . dt1, dt2t, uvar(nuvar,*),stifn(*),stifr(*), off(*), e(*),
86 . x(*), vx(*), vr(*),dx(*),dy(*),dz(*), rx(*),ry(*),rz(*),
87 . f(*) , forx(*), fory(*), forz(*),
88 . ar(*), xmom(*), ymom(*), zmom(*),
89 . fx1(mvsiz),fy1(mvsiz),fz1(mvsiz),
90 . fx2(mvsiz),fy2(mvsiz),fz2(mvsiz),
91 . mx1(mvsiz),my1(mvsiz),mz1(mvsiz),
92 . mx2(mvsiz),my2(mvsiz),mz2(mvsiz),rby(*),
93 . partsav(*),tani(15,*),ms(*),in(*),fsky(*),fskyv(*),gresav(*),
94 . msrt(*), dmelrt(*),x0_err(3,*)
95 DOUBLE PRECISION XDP(3,*)
96 TYPE(H3D_DATABASE) :: H3D_DATA
97 TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) , INTENT(IN) :: SENSOR_TAB
98
99
100
101 INTEGER ISENS
103 . umas(mvsiz),uiner(mvsiz),
104 . rot1(3,mvsiz),rot2(3,mvsiz),sti(3,mvsiz),stir(3,mvsiz),
105 . xkm(mvsiz),xkr(mvsiz),xcm(mvsiz),xcr(mvsiz)
106 DOUBLE PRECISION XL(MVSIZ,3)
107
108C
109 CALL rskew33(jft ,jlt ,ixr ,iout ,iprop,
110 . nuvar ,uvar ,rby ,x ,xl ,
111 . rot1 ,rot2 ,dx ,dy ,dz ,
112 . rx ,ry ,rz ,vr ,igtyp,
113 . nsensor,sensor_tab,isens ,nc1 ,nc2 ,
114 . xdp)
115
116 CALL ruser33(nel ,iout ,iprop ,nuvar ,uvar ,
117 . forx ,fory ,forz ,xmom ,ymom ,
118 . zmom ,xkm ,xkr ,xcm ,xcr ,
119 . xl ,umas ,uiner ,off ,e ,
120 . rot1 ,rot2 ,dx ,dy ,dz ,
121 . rx ,ry ,rz ,igtyp ,isens ,
122 . x0_err)
123
125 1 jft, jlt, dt2t, neltst,
126 2 ityptst, ixr, e, sti,
127 3 stir, off, xkm, xkr,
128 4 xcm, xcr, umas, uiner,
129 5 forx, fory, forz, xmom,
130 6 ymom, zmom, rot1, rot2,
131 7 msrt, dmelrt, nuvar, uvar,
132 8 igtyp, jsms)
133
135 1 jft, jlt, e, partsav,
136 2 umas, vx, ipartr, gresav,
137 3 grth, igrth, nc1, nc2,
138 4 igre)
139 CALL ranim33 (jft,jlt,forx,fory,forz,xmom,ymom,zmom,xl,tani,h3d_data)
140
141 IF(iparit==0)THEN
142 CALL rcum33(jft ,jlt ,xl ,nc1 ,nc2 ,
143 . f ,forx ,fory ,forz ,ar ,
144 . xmom ,ymom ,zmom ,sti ,stir ,
145 . stifn ,stifr ,fx1 ,fx2 ,fy1 ,
146 . fy2 ,fz1 ,fz2 ,mx1 ,mx2 ,
147 . my1 ,my2 ,mz1 ,mz2 ,nuvar ,
148 . uvar)
149 ELSE
151 1 jft, jlt, xl, fsky,
152 2 fskyv, forx, fory, forz,
153 3 xmom, ymom, zmom, sti,
154 4 stir, fx1, fx2, fy1,
155 5 fy2, fz1, fz2, mx1,
156 6 mx2, my1, my2, mz1,
157 7 mz2, iadr, nuvar, uvar,
158 8 nft)
159 ENDIF
160
161 RETURN
subroutine ranim33(jft, jlt, forx, fory, forz, xmom, ymom, zmom, xl, tani, h3d_data)
subroutine rbilan33(jft, jlt, eint, partsav, umas, v, ipartr, gresav, grth, igrth, nc1, nc2, igre)
subroutine rcum33(jft, jlt, xl, nc1, nc2, f, forx, fory, forz, xm, xmom, ymom, zmom, sti, stir, stifn, stifr, fx1, fx2, fy1, fy2, fz1, fz2, mx1, mx2, my1, my2, mz1, mz2, nuvar, uvar)
subroutine rcum33p(jft, jlt, xl, fsky, fskyv, forx, fory, forz, xmom, ymom, zmom, sti, stir, fx1, fx2, fy1, fy2, fz1, fz2, mx1, mx2, my1, my2, mz1, mz2, iadr, nuvar, uvar, nft)
subroutine rdtime33(jft, jlt, dt2t, neltst, ityptst, ixr, eint, sti, stir, off, xkm, xkr, xcm, xcr, umas, uiner, fx, fy, fz, xmom, ymom, zmom, rot1, rot2, msrt, dmelrt, nuvar, uvar, jntyp, jsms)
subroutine rskew33(jft, jlt, ixr, iout, iprop, nuvar, uvar, rby, x, xl, rot1, rot2, dx, dy, dz, rx, ry, rz, vr, igtyp, nsensor, sensor_tab, isens, nc1, nc2, xdp)
subroutine ruser33(nel, iout, iprop, nuvar, uvar, fx, fy, fz, xmom, ymom, zmom, xkm, xkr, xcm, xcr, xl, mass, iner, off, eint, rot1, rot2, dx, dy, dz, rx, ry, rz, igtyp, isens, x0_err)