OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
rbyfor.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| rbyfor ../engine/source/constraints/general/rbody/rbyfor.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| damping_vref_rby ../engine/source/assembly/damping_vref_rby.F90
29!|| damping_vref_rby_stiff ../engine/source/assembly/damping_vref_rby_stiff.F90
30!|| rgbodfp ../engine/source/constraints/general/rbody/rgbodfp.F
31!|| spmd_exch_a_rb6 ../engine/source/mpi/kinematic_conditions/spmd_exch_a_rb6.F
32!|| spmd_exch_a_rb6_vrel ../engine/source/mpi/kinematic_conditions/spmd_exch_a_rb6.F
33!|| spmd_exch_a_rb6g ../engine/source/mpi/kinematic_conditions/spmd_exch_a_rb6g.F
34!|| startime ../engine/source/system/timer_mod.F90
35!|| stoptime ../engine/source/system/timer_mod.F90
36!||--- uses -----------------------------------------------------
37!|| damping_vref_rby_mod ../engine/source/assembly/damping_vref_rby.F90
38!|| damping_vref_rby_stiff_mod ../engine/source/assembly/damping_vref_rby_stiff.F90
39!|| groupdef_mod ../common_source/modules/groupdef_mod.F
40!|| timer_mod ../engine/source/system/timer_mod.F90
41!||====================================================================
42 SUBROUTINE rbyfor(TIMERS,
43 1 RBY ,A ,AR ,X ,VR ,
44 2 FSAV ,IN ,STIFN ,STIFR ,FANI ,
45 3 LPBY ,NPBY ,WEIGHT ,MS ,V ,
46 4 IGRSURF ,BUFSF ,ICODR ,ISKEW ,SKEW ,
47 5 KIND ,IAD_RBY,FR_RBY6,RBY6 ,IRBKIN_L,
48 6 NRBYKIN_L,NATIV_SMS,DIMFB,FBSAV6 ,STABSEN ,
49 7 TABSENSOR,NODREAC,FTHREAC,CPTREAC,DAMPR ,
50 8 SDAMP ,DAMP ,NDAMP_VREL,ID_DAMP_VREL,IGRNOD,
51 9 TAGSLV_RBY,IPARIT,WFEXT,NDAMP_VREL_RBYG,SIZE_RBY6_C,
52 A RBY6_C ,NHIER_RBY)
53C-----------------------------------------------
54C M o d u l e s
55C-----------------------------------------------
56 USE timer_mod
57 USE groupdef_mod
58 USE damping_vref_rby_mod
59 USE damping_vref_rby_stiff_mod
60C----6---------------------------------------------------------------7---------8
61C I m p l i c i t T y p e s
62C-----------------------------------------------
63#include "implicit_f.inc"
64#include "comlock.inc"
65C-----------------------------------------------
66C C o m m o n B l o c k s
67C-----------------------------------------------
68#include "com01_c.inc"
69#include "com04_c.inc"
70#include "com08_c.inc"
71#include "param_c.inc"
72#include "task_c.inc"
73#include "timeri_c.inc"
74#include "intstamp_c.inc"
75C-----------------------------------------------------------------
76C D u m m y A r g u m e n t s
77C-----------------------------------------------
78 INTEGER, INTENT(IN) :: NHIER_RBY
79 TYPE(TIMER_), INTENT(inout) :: TIMERS
80 INTEGER ICODR(*),ISKEW(*),WEIGHT(*),LPBY(*),NPBY(NNPBY,*),
81 . KIND(NRBYKIN),FR_RBY6(*),IAD_RBY(*),IRBKIN_L(*),
82 . NRBYKIN_L, NATIV_SMS(*),DIMFB,STABSEN,TABSENSOR(*),
83 . NODREAC(*),CPTREAC,NDAMP_VREL_RBYG,SIZE_RBY6_C
84 my_real RBY(NRBY,*) ,A(3,*) ,AR(3,*) ,X(3,*) ,VR(3,*),FSAV(NTHVKI,*),
85 . IN(*) ,STIFN(*),STIFR(*),FANI(3,*),MS(*),V(3,*) ,
86 . BUFSF(*), SKEW(LSKEW,*),FTHREAC(6,*)
87 DOUBLE PRECISION RBY6(8,6,NRBYKIN),RBY6_C(2,6,SIZE_RBY6_C)
88 DOUBLE PRECISION FBSAV6(12,6,DIMFB),RBID(12,6),WFEXT
89 TYPE (SURF_), DIMENSION(NSURF) :: IGRSURF
90 INTEGER ,INTENT(IN) :: NDAMP_VREL,IPARIT,SDAMP
91 INTEGER ,INTENT(IN) :: ID_DAMP_VREL(NDAMP_VREL),TAGSLV_RBY(NUMNOD)
92 my_real, INTENT(INOUT) :: DAMPR(NRDAMP,NDAMP),DAMP(SDAMP)
93 TYPE(GROUP_), INTENT(IN) :: IGRNOD(NGRNOD)
94C-----------------------------------------------
95C L o c a l V a r i a b l e s
96C-----------------------------------------------
97 INTEGER J,K,N,KK,N2,
98 . isu,adrsrf,im, nrbdim,iparsens,isect,
99 . dim,nhi
100C-------------------------------------
101C INIT RBY6 to ZERO for deactivated Rigid Bodies (non multi-thread)
102C-------------------------------------
103!$OMP SINGLE
104 rbid = zero
105 DO kk=1,nrbykin_l
106 n=irbkin_l(kk)
107 IF(npby(7,n)<=0) THEN
108 DO k = 1, 6
109 rby6(1,k,n) = zero
110 rby6(2,k,n) = zero
111 rby6(3,k,n) = zero
112 rby6(4,k,n) = zero
113 rby6(5,k,n) = zero
114 rby6(6,k,n) = zero
115 rby6(7,k,n) = zero
116 rby6(8,k,n) = zero
117 ENDDO
118 ENDIF
119 ENDDO
120!$OMP END SINGLE
121 IF (nhier_rby ==0) THEN
122C-------------------------------------
123C FORCE COMPUTATION FOR SUPER RIGID BODIES (non multi-thread)
124C-------------------------------------
125
126!$OMP SINGLE
127
128 DO kk=1,nrbykin_l
129 n=irbkin_l(kk)
130 k = kind(n)
131 IF (npby(4,n)/=0) THEN
132 IF(npby(7,n)>0) THEN
133 n2 = ninter+nrwall+n
134C sum of forces / moments on the surface:
135C at each time step (and at TT=0. in the starter),
136C the rigid body transmits the position of the main node
137C to the surface
138C ET
139C in interface type 14, the forces, moments (and stiffnesses)
140C are computed at this point!!!
141 isu=npby(8,n)
142 IF(ispmd==0.AND.isu/=0) THEN
143 im=npby(1,n)
144 adrsrf=igrsurf(isu)%IAD_BUFR
145 a(1,im) =a(1,im) +bufsf(adrsrf+25)
146 a(2,im) =a(2,im) +bufsf(adrsrf+26)
147 a(3,im) =a(3,im) +bufsf(adrsrf+27)
148 ar(1,im)=ar(1,im)+bufsf(adrsrf+28)
149 ar(2,im)=ar(2,im)+bufsf(adrsrf+29)
150 ar(3,im)=ar(3,im)+bufsf(adrsrf+30)
151 stifn(im)=stifn(im)+bufsf(adrsrf+31)
152 stifr(im)=stifr(im)+bufsf(adrsrf+32)
153C output to th of interface forces.
154 fsav(10,n2)=fsav(10,n2)+bufsf(adrsrf+25)*dt1
155 fsav(11,n2)=fsav(11,n2)+bufsf(adrsrf+26)*dt1
156 fsav(12,n2)=fsav(12,n2)+bufsf(adrsrf+27)*dt1
157 fsav(13,n2)=fsav(13,n2)+bufsf(adrsrf+28)*dt1
158 fsav(14,n2)=fsav(14,n2)+bufsf(adrsrf+29)*dt1
159 fsav(15,n2)=fsav(15,n2)+bufsf(adrsrf+30)*dt1
160 END IF
161c
162 iparsens=0
163 isect=0
164 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
165 . tabsensor(n+nsect+nintsub+ninter+nrwall)
166 IF(isect/=0) THEN
167 iparsens=1
168 CALL rgbodfp(
169 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
170 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
171 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
172 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
173 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
174 6 nodreac,fthreac,cptreac,npby(5,n))
175 ELSE
176 CALL rgbodfp(
177 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
178 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
179 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
180 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
181 5 npby(2,n),nativ_sms,rbid , iparsens,
182 6 nodreac,fthreac,cptreac,npby(5,n))
183 ENDIF
184 ENDIF
185 ENDIF
186 ENDDO
187
188!$OMP END SINGLE
189
190C-------------------------------------
191C FORCE COMPUTATION FOR CLASSICAL RIGID BODIES (multi-thread)
192C-------------------------------------
193
194C-------------------------------------
195C Processing rigid bodies (without packet), with dynamic scheduling
196C-------------------------------------
197
198!$OMP DO SCHEDULE(DYNAMIC,1)
199
200 DO kk=1,nrbykin_l
201 n = irbkin_l(kk)
202 k = kind(n)
203 IF (npby(4,n)==0) THEN
204 IF(npby(7,n)>0) THEN
205 n2 = ninter+nrwall+n
206C sum of forces / moments on the surface:
207C at each time step (and at tt=0. in the starter),
208C the rigid body transmits the position of the main node
209C to the surface
210C ET
211C in interface type 14, the forces, moments (and stiffnesses)
212C are computed at this point!!!
213 isu=npby(8,n)
214 IF(ispmd==0.AND.isu/=0) THEN
215 im=npby(1,n)
216 adrsrf=igrsurf(isu)%IAD_BUFR
217 a(1,im) =a(1,im) +bufsf(adrsrf+25)
218 a(2,im) =a(2,im) +bufsf(adrsrf+26)
219 a(3,im) =a(3,im) +bufsf(adrsrf+27)
220 ar(1,im)=ar(1,im)+bufsf(adrsrf+28)
221 ar(2,im)=ar(2,im)+bufsf(adrsrf+29)
222 ar(3,im)=ar(3,im)+bufsf(adrsrf+30)
223 stifn(im)=stifn(im)+bufsf(adrsrf+31)
224 stifr(im)=stifr(im)+bufsf(adrsrf+32)
225C output to th of interface forces.
226 fsav(10,n2)=fsav(10,n2)+bufsf(adrsrf+25)*dt1
227 fsav(11,n2)=fsav(11,n2)+bufsf(adrsrf+26)*dt1
228 fsav(12,n2)=fsav(12,n2)+bufsf(adrsrf+27)*dt1
229 fsav(13,n2)=fsav(13,n2)+bufsf(adrsrf+28)*dt1
230 fsav(14,n2)=fsav(14,n2)+bufsf(adrsrf+29)*dt1
231 fsav(15,n2)=fsav(15,n2)+bufsf(adrsrf+30)*dt1
232 END IF
233c
234 iparsens=0
235 isect=0
236 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
237 . tabsensor(n+nsect+nintsub+ninter+nrwall)
238 IF(isect/=0) THEN
239 iparsens=1
240 CALL rgbodfp(
241 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
242 2 lpby(k) ,npby(1,n) ,in ,vr ,stifn ,
243 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
244 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
245 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
246 6 nodreac,fthreac,cptreac,npby(5,n))
247 ELSE
248 CALL rgbodfp(
249 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
250 2 lpby(k) ,npby(1,n) ,in ,vr ,stifn ,
251 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
252 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
253 5 npby(2,n),nativ_sms ,rbid ,iparsens ,
254 6 nodreac,fthreac,cptreac,npby(5,n))
255 ENDIF
256 ENDIF
257 ENDIF
258 ENDDO
259
260!$omp END DO
261
262C-------------------------------------
263C Contribution of damping forces from /DAMP/VREL
264C-------------------------------------
265 IF ((ndamp_vrel_rbyg > 0).AND.(nrbykin_l > 0)) THEN
266!$OMP SINGLE
267 dim = 3+3*iroddl
268 CALL damping_vref_rby(igrnod ,ngrnod ,v ,vr ,a ,
269 . x ,ms ,dampr ,nrdamp ,ndamp ,
270 . ndamp_vrel,iparit ,numnod ,dt1 ,id_damp_vrel,
271 . tt ,nnpby ,nrbykin ,npby ,rby6 ,
272 . rby6_c ,tagslv_rby,weight ,lskew ,numskw ,
273 . dim ,damp ,skew ,wfext ,size_rby6_c,
274 . nhier_rby )
275!$OMP END SINGLE
276 ENDIF
277
278 IF (nspmd > 1) THEN
279
280!$OMP SINGLE
281
282 IF (imon>0) CALL startime(timers,11)
283C
284C Comm non multi-thread
285C
286C flag to enable/disable all to all communication for RBY.
287 IF(nintstamp == 0) THEN
288 IF (ndamp_vrel_rbyg == 0) THEN
289C routine A2A for Parith/ON computation of A AR STIFN STIFR on main nodes of RB
290 nrbdim=8
291 CALL spmd_exch_a_rb6(nrbdim,iad_rby,fr_rby6,iad_rby(nspmd+1),rby6)
292 ELSE
293C routine A2A for Parith/ON computation of A AR STIFN STIFR C CR on main nodes of RB for /DAMP/VREL
294 nrbdim=10
295 CALL spmd_exch_a_rb6_vrel(nrbdim,iad_rby,fr_rby6,iad_rby(nspmd+1),rby6,
296 . rby6_c,size_rby6_c)
297 ENDIF
298C routine G/S for Parith/ON computation of A AR STIFN STIFR main nodes RB
299 ELSE
300 CALL spmd_exch_a_rb6g(npby,rby6)
301 END IF
302
303 IF (imon>0) CALL stoptime(timers,11)
304
305!$OMP END SINGLE
306
307 END IF
308C
309C End of parith/on multi-thread processing
310C
311!$OMP DO SCHEDULE(DYNAMIC,1)
312
313 DO kk=1,nrbykin_l
314 n = irbkin_l(kk)
315 k = kind(n)
316 n2 = ninter+nrwall+n
317 IF(npby(7,n)>0) THEN
318c
319 iparsens=0
320 isect=0
321 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
322 . tabsensor(n+nsect+nintsub+ninter+nrwall)
323 IF(isect/=0) THEN
324 iparsens=1
325 CALL rgbodfp(
326 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
327 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
328 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
329 4 2 ,icodr ,iskew ,skew, rby6(1,1,n),
330 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
331 6 nodreac,fthreac,cptreac,npby(5,n))
332 ELSE
333 CALL rgbodfp(
334 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
335 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
336 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
337 4 2 ,icodr ,iskew ,skew, rby6(1,1,n),
338 5 npby(2,n),nativ_sms,rbid ,iparsens,
339 6 nodreac,fthreac,cptreac,npby(5,n))
340 ENDIF
341 ENDIF
342 ENDDO
343
344!$OMP END DO
345C-------------------------------------
346C Computation of stability of main node for DAMP/VREL
347C-------------------------------------
348 IF (ndamp_vrel_rbyg > 0) THEN
349!$OMP SINGLE
350 call damping_vref_rby_stiff(numnod,nnpby,nrbykin,nrbykin_l,npby,
351 . rby6_c,ms,in,stifn,stifr,size_rby6_c,
352 . irbkin_l ,nhier_rby)
353!$OMP END SINGLE
354 ENDIF
355C
356 ELSE ! NHIER_RBY>0
357!$OMP SINGLE
358 DO nhi = 0 ,nhier_rby
359 DO n=1,nrbykin
360 im = npby(1,n)
361 IF(npby(7,n)==0 .OR. im<=0 .OR. npby(20,n)/=nhi) cycle
362 k = kind(n)
363 n2 = ninter+nrwall+n
364 isu=npby(8,n)
365 IF(ispmd==0.AND.isu/=0) THEN
366 im=npby(1,n)
367 adrsrf=igrsurf(isu)%IAD_BUFR
368 a(1,im) =a(1,im) +bufsf(adrsrf+25)
369 a(2,im) =a(2,im) +bufsf(adrsrf+26)
370 a(3,im) =a(3,im) +bufsf(adrsrf+27)
371 ar(1,im)=ar(1,im)+bufsf(adrsrf+28)
372 ar(2,im)=ar(2,im)+bufsf(adrsrf+29)
373 ar(3,im)=ar(3,im)+bufsf(adrsrf+30)
374 stifn(im)=stifn(im)+bufsf(adrsrf+31)
375 stifr(im)=stifr(im)+bufsf(adrsrf+32)
376C output to th of interface forces.
377 fsav(10,n2)=fsav(10,n2)+bufsf(adrsrf+25)*dt1
378 fsav(11,n2)=fsav(11,n2)+bufsf(adrsrf+26)*dt1
379 fsav(12,n2)=fsav(12,n2)+bufsf(adrsrf+27)*dt1
380 fsav(13,n2)=fsav(13,n2)+bufsf(adrsrf+28)*dt1
381 fsav(14,n2)=fsav(14,n2)+bufsf(adrsrf+29)*dt1
382 fsav(15,n2)=fsav(15,n2)+bufsf(adrsrf+30)*dt1
383 END IF
384c
385 iparsens=0
386 isect=0
387 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
388 . tabsensor(n+nsect+nintsub+ninter+nrwall)
389 IF(isect/=0) THEN
390 iparsens=1
391 CALL rgbodfp(
392 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
393 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
394 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
395 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
396 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
397 6 nodreac,fthreac,cptreac,npby(5,n))
398 ELSE
399 CALL rgbodfp(
400 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
401 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
402 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
403 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
404 5 npby(2,n),nativ_sms,rbid , iparsens,
405 6 nodreac,fthreac,cptreac,npby(5,n))
406 END IF
407 END DO !N=1,NRBYKIN
408C-------------------------------------
409C Contribution of damping forces from /DAMP/VREL
410C-------------------------------------
411 IF ((ndamp_vrel_rbyg > 0).AND.(nrbykin_l > 0)) THEN
412 dim = 3+3*iroddl
413 CALL damping_vref_rby(igrnod ,ngrnod ,v ,vr ,a ,
414 . x ,ms ,dampr ,nrdamp ,ndamp ,
415 . ndamp_vrel,iparit ,numnod ,dt1 ,id_damp_vrel,
416 . tt ,nnpby ,nrbykin ,npby ,rby6 ,
417 . rby6_c ,tagslv_rby,weight ,lskew ,numskw ,
418 . dim ,damp ,skew ,wfext ,size_rby6_c,
419 . nhi )
420 ENDIF
421
422 IF (nspmd > 1) THEN
423 IF (imon>0) CALL startime(timers,11)
424 IF(nintstamp == 0) THEN
425 IF (ndamp_vrel_rbyg == 0) THEN
426 nrbdim=8
427 CALL spmd_exch_a_rb6(nrbdim,iad_rby,fr_rby6,iad_rby(nspmd+1),rby6)
428 ELSE
429 nrbdim=10
430 CALL spmd_exch_a_rb6_vrel(nrbdim,iad_rby,fr_rby6,iad_rby(nspmd+1),rby6,
431 . rby6_c,size_rby6_c)
432 ENDIF
433 ELSE
434 CALL spmd_exch_a_rb6g(npby,rby6)
435 END IF
436 IF (imon>0) CALL stoptime(timers,11)
437 END IF
438! back to main node : IFLAG=2
439 DO n=1,nrbykin
440 im = npby(1,n)
441 k = kind(n)
442 n2 = ninter+nrwall+n
443 IF(npby(7,n)==0 .OR. im<=0 .OR. npby(20,n)/=nhi) cycle
444!
445 iparsens=0
446 isect=0
447 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
448 . tabsensor(n+nsect+nintsub+ninter+nrwall)
449 IF(isect/=0) THEN
450 iparsens=1
451 CALL rgbodfp(
452 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
453 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
454 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
455 4 2 ,icodr ,iskew ,skew, rby6(1,1,n),
456 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
457 6 nodreac,fthreac,cptreac,npby(5,n))
458 ELSE
459 CALL rgbodfp(
460 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
461 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
462 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
463 4 2 ,icodr ,iskew ,skew, rby6(1,1,n),
464 5 npby(2,n),nativ_sms,rbid ,iparsens,
465 6 nodreac,fthreac,cptreac,npby(5,n))
466 END IF
467 END DO
468C-----------------------------------
469C Computation of stability of main node for DAMP/VREL
470C-----------------------------------
471 IF (ndamp_vrel_rbyg > 0) THEN
472 call damping_vref_rby_stiff(numnod,nnpby,nrbykin,nrbykin_l,npby,
473 . rby6_c,ms,in,stifn,stifr,size_rby6_c,
474 . irbkin_l ,nhi )
475 ENDIF
476 END DO !NHI = ,NHIER
477!$OMP END SINGLE
478 END IF !IF (NHIER ==0)
479C
480 RETURN
481 END
groupdef_mod
Definition groupdef_mod.F:668
rbyfor
subroutine rbyfor(timers, rby, a, ar, x, vr, fsav, in, stifn, stifr, fani, lpby, npby, weight, ms, v, igrsurf, bufsf, icodr, iskew, skew, kind, iad_rby, fr_rby6, rby6, irbkin_l, nrbykin_l, nativ_sms, dimfb, fbsav6, stabsen, tabsensor, nodreac, fthreac, cptreac, dampr, sdamp, damp, ndamp_vrel, id_damp_vrel, igrnod, tagslv_rby, iparit, wfext, ndamp_vrel_rbyg, size_rby6_c, rby6_c, nhier_rby)
Definition rbyfor.F:53
rgbodfp
subroutine rgbodfp(af, am, x, fs, rby, nod, m, in, vr, stifn, stifr, fopta, weight, ms, v, iflag, icodr, iskew, skew, rbf6, nsn, nativ_sms, fbsav6, iparsens, nodreac, fthreac, cptreac, ispher)
Definition rgbodfp.F:39
spmd_exch_a_rb6_vrel
subroutine spmd_exch_a_rb6_vrel(nrbdim, iad_rby, fr_rby6, icsize, rby6, rby6_c, size_rby6_c)
Definition spmd_exch_a_rb6.F:170
spmd_exch_a_rb6
subroutine spmd_exch_a_rb6(nrbdim, iad_rby, fr_rby6, icsize, rbf6)
Definition spmd_exch_a_rb6.F:37
spmd_exch_a_rb6g
subroutine spmd_exch_a_rb6g(npby, rbf6)
Definition spmd_exch_a_rb6g.F:34
startime
subroutine startime(event, itask)
Definition timer.F:93
stoptime
subroutine stoptime(event, itask)
Definition timer.F:135