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