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)
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 TYPE(timer_), INTENT(inout) :: TIMERS
79 INTEGER ICODR(*),ISKEW(*),WEIGHT(*),LPBY(*),NPBY(NNPBY,*),
80 . KIND(NRBYKIN),FR_RBY6(*),IAD_RBY(*),IRBKIN_L(*),
81 . NRBYKIN_L, NATIV_SMS(*),DIMFB,STABSEN,TABSENSOR(*),
82 . NODREAC(*),CPTREAC,NDAMP_VREL_RBYG,SIZE_RBY6_C
83 my_real RBY(NRBY,*) ,A(3,*) ,AR(3,*) ,X(3,*) ,VR(3,*),FSAV(NTHVKI,*),
84 . IN(*) ,STIFN(*),STIFR(*),FANI(3,*),MS(*),V(3,*) ,
85 . BUFSF(*), SKEW(LSKEW,*),FTHREAC(6,*)
86 DOUBLE PRECISION RBY6(8,6,NRBYKIN),RBY6_C(2,6,SIZE_RBY6_C)
87 DOUBLE PRECISION FBSAV6(12,6,DIMFB),RBID(12,6),WFEXT
88 TYPE (SURF_), DIMENSION(NSURF) :: IGRSURF
89 INTEGER ,INTENT(IN) :: NDAMP_VREL,IPARIT,SDAMP
90 INTEGER ,INTENT(IN) :: ID_DAMP_VREL(NDAMP_VREL),TAGSLV_RBY(NUMNOD)
91 my_real, INTENT(INOUT) :: DAMPR(NRDAMP,NDAMP),DAMP(SDAMP)
92 TYPE(GROUP_), INTENT(IN) :: IGRNOD(NGRNOD)
93C-----------------------------------------------
94C L o c a l V a r i a b l e s
95C-----------------------------------------------
96 INTEGER J,K,N,KK,N2,ILAGM,
97 . isu,adrsrf,im, l, lsend, lrecv, nrbdim,iparsens,isect,
98 . dim
99C-------------------------------------
100C INIT RBY6 to ZERO for deactivated Rigid Bodies (non multi-thread)
101C-------------------------------------
102!$OMP SINGLE
103 rbid = zero
104 DO kk=1,nrbykin_l
105 n=irbkin_l(kk)
106 IF(npby(7,n)<=0) THEN
107 DO k = 1, 6
108 rby6(1,k,n) = zero
109 rby6(2,k,n) = zero
110 rby6(3,k,n) = zero
111 rby6(4,k,n) = zero
112 rby6(5,k,n) = zero
113 rby6(6,k,n) = zero
114 rby6(7,k,n) = zero
115 rby6(8,k,n) = zero
116 ENDDO
117 ENDIF
118 ENDDO
119!$OMP END SINGLE
120C-------------------------------------
121C CALCUL FORCE SUPER RIGID BODIES (non multi-thread)
122C-------------------------------------
123
124!$OMP SINGLE
125
126 DO kk=1,nrbykin_l
127 n=irbkin_l(kk)
128 k = kind(n)
129 IF (npby(4,n)/=0) THEN
130 IF(npby(7,n)>0) THEN
131 n2 = ninter+nrwall+n
132C Cumul des Forces / Moments sur la surface :
133C A chaque instant (et a TT=0. dans le Starter),
134C le Rigid Body transmet la position du noeud main
135C a la surface
136C ET
137C Dans l'Interface TYPE 14, Les forces, moments (et rigidites)
138C sont calcules en ce point !!!
139 isu=npby(8,n)
140 IF(ispmd==0.AND.isu/=0) THEN
141 im=npby(1,n)
142 adrsrf=igrsurf(isu)%IAD_BUFR
143 a(1,im) =a(1,im) +bufsf(adrsrf+25)
144 a(2,im) =a(2,im) +bufsf(adrsrf+26)
145 a(3,im) =a(3,im) +bufsf(adrsrf+27)
146 ar(1,im)=ar(1,im)+bufsf(adrsrf+28)
147 ar(2,im)=ar(2,im)+bufsf(adrsrf+29)
148 ar(3,im)=ar(3,im)+bufsf(adrsrf+30)
149 stifn(im)=stifn(im)+bufsf(adrsrf+31)
150 stifr(im)=stifr(im)+bufsf(adrsrf+32)
151C Sortie sur TH des forces d'interface.
152 fsav(10,n2)=fsav(10,n2)+bufsf(adrsrf+25)*dt1
153 fsav(11,n2)=fsav(11,n2)+bufsf(adrsrf+26)*dt1
154 fsav(12,n2)=fsav(12,n2)+bufsf(adrsrf+27)*dt1
155 fsav(13,n2)=fsav(13,n2)+bufsf(adrsrf+28)*dt1
156 fsav(14,n2)=fsav(14,n2)+bufsf(adrsrf+29)*dt1
157 fsav(15,n2)=fsav(15,n2)+bufsf(adrsrf+30)*dt1
158 END IF
159c
160 iparsens=0
161 isect=0
162 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
163 . tabsensor(n+nsect+nintsub+ninter+nrwall)
164 IF(isect/=0) THEN
165 iparsens=1
166 CALL rgbodfp(
167 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
168 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
169 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
170 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
171 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
172 6 nodreac,fthreac,cptreac,npby(5,n))
173 ELSE
174 CALL rgbodfp(
175 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
176 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
177 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
178 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
179 5 npby(2,n),nativ_sms,rbid , iparsens,
180 6 nodreac,fthreac,cptreac,npby(5,n))
181 ENDIF
182 ENDIF
183 ENDIF
184 ENDDO
185
186!$OMP END SINGLE
187
188C-------------------------------------
189C CALCUL FORCE RIGID BODIES CLASSIQUES (multi-thread)
190C-------------------------------------
191
192C-------------------------------------
193C Traitement rigid bodies (sans packet), avec scheduling dynamique
194C-------------------------------------
195
196!$OMP DO SCHEDULE(DYNAMIC,1)
197
198 DO kk=1,nrbykin_l
199 n = irbkin_l(kk)
200 k = kind(n)
201 IF (npby(4,n)==0) THEN
202 IF(npby(7,n)>0) THEN
203 n2 = ninter+nrwall+n
204C Cumul des Forces / Moments sur la surface :
205C A chaque instant (et a TT=0. dans le Starter),
206C le Rigid Body transmet la position du noeud main
207C a la surface
208C ET
209C Dans l'Interface TYPE 14, Les forces, moments (et rigidites)
210C sont calcules en ce point !!!
211 isu=npby(8,n)
212 IF(ispmd==0.AND.isu/=0) THEN
213 im=npby(1,n)
214 adrsrf=igrsurf(isu)%IAD_BUFR
215 a(1,im) =a(1,im) +bufsf(adrsrf+25)
216 a(2,im) =a(2,im) +bufsf(adrsrf+26)
217 a(3,im) =a(3,im) +bufsf(adrsrf+27)
218 ar(1,im)=ar(1,im)+bufsf(adrsrf+28)
219 ar(2,im)=ar(2,im)+bufsf(adrsrf+29)
220 ar(3,im)=ar(3,im)+bufsf(adrsrf+30)
221 stifn(im)=stifn(im)+bufsf(adrsrf+31)
222 stifr(im)=stifr(im)+bufsf(adrsrf+32)
223C Sortie sur TH des forces d'interface.
224 fsav(10,n2)=fsav(10,n2)+bufsf(adrsrf+25)*dt1
225 fsav(11,n2)=fsav(11,n2)+bufsf(adrsrf+26)*dt1
226 fsav(12,n2)=fsav(12,n2)+bufsf(adrsrf+27)*dt1
227 fsav(13,n2)=fsav(13,n2)+bufsf(adrsrf+28)*dt1
228 fsav(14,n2)=fsav(14,n2)+bufsf(adrsrf+29)*dt1
229 fsav(15,n2)=fsav(15,n2)+bufsf(adrsrf+30)*dt1
230 END IF
231c
232 iparsens=0
233 isect=0
234 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
235 . tabsensor(n+nsect+nintsub+ninter+nrwall)
236 IF(isect/=0) THEN
237 iparsens=1
238 CALL rgbodfp(
239 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
240 2 lpby(k) ,npby(1,n) ,in ,vr ,stifn ,
241 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
242 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
243 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
244 6 nodreac,fthreac,cptreac,npby(5,n))
245 ELSE
246 CALL rgbodfp(
247 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
248 2 lpby(k) ,npby(1,n) ,in ,vr ,stifn ,
249 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
250 4 1 ,icodr ,iskew ,skew ,rby6(1,1,n),
251 5 npby(2,n),nativ_sms ,rbid ,iparsens ,
252 6 nodreac,fthreac,cptreac,npby(5,n))
253 ENDIF
254 ENDIF
255 ENDIF
256 ENDDO
257
258!$OMP END DO
259
260C-------------------------------------
261C Contribution of damping forces from /DAMP/VREL
262C-------------------------------------
263 IF ((ndamp_vrel_rbyg > 0).AND.(nrbykin_l > 0)) THEN
264!$OMP SINGLE
265 dim = 3+3*iroddl
266 CALL damping_vref_rby(igrnod ,ngrnod ,v ,vr ,a ,
267 . x ,ms ,dampr ,nrdamp ,ndamp ,
268 . ndamp_vrel,iparit ,numnod ,dt1 ,id_damp_vrel,
269 . tt ,nnpby ,nrbykin ,npby ,rby6 ,
270 . rby6_c ,tagslv_rby,weight ,lskew ,numskw ,
271 . dim ,damp ,skew ,wfext ,size_rby6_c)
272!$OMP END SINGLE
273 ENDIF
274
275 IF (nspmd > 1) THEN
276
277!$OMP SINGLE
278
279 IF (imon>0) CALL startime(timers,11)
280C
281C Comm non multi-thread
282C
283C flag to enable/disable all to all communication for RBY.
284 IF(nintstamp == 0) THEN
285 IF (ndamp_vrel_rbyg == 0) THEN
286C routine A2A for Parith/ON computation of A AR STIFN STIFR on main nodes of RB
287 nrbdim=8
288 CALL spmd_exch_a_rb6(nrbdim,iad_rby,fr_rby6,iad_rby(nspmd+1),rby6)
289 ELSE
290C routine A2A for Parith/ON computation of A AR STIFN STIFR C CR on main nodes of RB for /DAMP/VREL
291 nrbdim=10
292 CALL spmd_exch_a_rb6_vrel(nrbdim,iad_rby,fr_rby6,iad_rby(nspmd+1),rby6,
293 . rby6_c,size_rby6_c)
294 ENDIF
295C routine G/S de calcul Parith/ON de A AR STIFN STIFR noeuds main RB
296 ELSE
297 CALL spmd_exch_a_rb6g(npby,rby6)
298 END IF
299
300 IF (imon>0) CALL stoptime(timers,11)
301
302!$OMP END SINGLE
303
304 END IF
305C
306C Traitement fin parith/on multi-thread
307C
308!$OMP DO SCHEDULE(DYNAMIC,1)
309
310 DO kk=1,nrbykin_l
311 n = irbkin_l(kk)
312 k = kind(n)
313 n2 = ninter+nrwall+n
314 IF(npby(7,n)>0) THEN
315c
316 iparsens=0
317 isect=0
318 IF(stabsen/=0) isect=tabsensor(n+nsect+nintsub+ninter+nrwall+1)-
319 . tabsensor(n+nsect+nintsub+ninter+nrwall)
320 IF(isect/=0) THEN
321 iparsens=1
322 CALL rgbodfp(
323 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
324 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
325 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
326 4 2 ,icodr ,iskew ,skew, rby6(1,1,n),
327 5 npby(2,n),nativ_sms,fbsav6(1,1,isect) ,iparsens,
328 6 nodreac,fthreac,cptreac,npby(5,n))
329 ELSE
330 CALL rgbodfp(
331 1 a ,ar ,x ,fsav(1,n2),rby(1,n),
332 2 lpby(k) ,npby(1,n),in ,vr ,stifn ,
333 3 stifr ,fani(1,1+2*(nsect+n-1)),weight,ms ,v ,
334 4 2 ,icodr ,iskew ,skew, rby6(1,1,n),
335 5 npby(2,n),nativ_sms,rbid ,iparsens,
336 6 nodreac,fthreac,cptreac,npby(5,n))
337 ENDIF
338 ENDIF
339 ENDDO
340
341!$OMP END DO
342C
343C-------------------------------------
344C Computation of stability of main node for DAMP/VREL
345C-------------------------------------
346 IF (ndamp_vrel_rbyg > 0) THEN
347!$OMP SINGLE
348 call damping_vref_rby_stiff(numnod,nnpby,nrbykin,nrbykin_l,npby,
349 . rby6_c,ms,in,stifn,stifr,size_rby6_c,
350 . irbkin_l )
351!$OMP END SINGLE
352 ENDIF
353C
354 RETURN
355 END
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)
Definition rbyfor.F:53
subroutine resol(timers, element, nodes, coupling, af, iaf, iskwn, neth, ipart, nom_opt, kxx, ixx, ixtg, ixs, ixq, ixt, ixp, ixr, ifill, mat_elem, ims, npc, ibcl, ibfv, idum, las, laccelm, nnlink, lnlink, iparg, dd_iad, igrv, iexlnk, kinet, ipari, nprw, iconx, npby, lpby, lrivet, nstrf, ljoint, nodpor, monvol, ilink, llink, linale, neflsw, nnflsw, icut, cluster, itask, inoise, thke, damp, pm, skews, geo, eani, bufmat, bufgeo, bufsf, w, veul, fill, dfill, alph, wb, dsave, asave, msnf, tf, forc, vel, fsav, fzero, xlas, accelm, agrv, fr_wave, failwave, parts0, elbuf, rwbuf, sensors, rwsav, rby, rivet, secbuf, volmon, lambda, wa, fv, partsav, uwa, val2, phi, segvar, r, crflsw, flsw, fani, xcut, anin, tani, secfcum, bufnois, idata, rdata, iframe, kxsp, ixsp, nod2sp, ispsym, ispcond, xframe, spbuf, xspsym, vspsym, pv, fsavd, ibvel, lbvel, wasph, w16, isphio, lprtsph, lonfsph, vsphio, fbvel, lagbuf, ibcslag, iactiv, dampr, gjbufi, gjbufr, rbmpc, ibmpc, sphveln, nbrcvois, nbsdvois, lnrcvois, lnsdvois, nercvois, nesdvois, lercvois, lesdvois, npsegcom, lsegcom, nporgeo, ixtg1, npbyl, lpbyl, rbyl, igeo, ipm, madprt, madsh4, madsh3, madsol, madnod, madfail, iad_rby, fr_rby, fr_wall, iad_rby2, fr_rby2, iad_i2m, fr_i2m, addcni2, procni2, iadi2, fr_mv, iadmv2, fr_ll, fr_rl, iadcj, fr_cj, fr_sec, iad_sec, iad_cut, fr_cut, rg_cut, newfront, fr_mad, fxbipm, fxbrpm, fxbnod, fxbmod, fxbglm, fxbcpm, fxbcps, fxblm, fxbfls, fxbdls, fxbdep, fxbvit, fxbacc, fxbelm, fxbsig, fxbgrvi, fxbgrvr, eigipm, eigibuf, eigrpm, lnodpor, fr_i18, graphe, iflow, rflow, lgrav, dd_r2r, fasolfr, fr_lagf, llagf, lprw, icontact, rcontact, sh4tree, sh3tree, ipadmesh, padmesh, msc, mstg, inc, intg, ptg, iskwp, nskwp, isensp, nsensp, iaccp, naccp, ipart_state, acontact, pcontact, factiv, sh4trim, sh3trim, mscnd, incnd, ibfflux, fbfflux, rbym, irbym, lnrbym, icodrbym, ibcv, fconv, ibftemp, fbftemp, iad_rbym, fr_rbym, weight_rm, ms_ply, zi_ply, inod_pxfem, iel_pxfem, iadc_pxfem, adsky_pxfem, icode_ply, icodt_ply, iskew_ply, admsms, madclnod, nom_sect, mcpc, mcptg, dmelc, dmeltg, mssa, dmels, mstr, dmeltr, msp, dmelp, msrt, dmelrt, ibcr, fradia, res_sms, table, irbe2, lrbe2, iad_rbe2, fr_rbe2, phie, msf, procne_pxfem, iadsdp_pxfem, iadrcp_pxfem, icfield, lcfield, cfield, msz2, diag_sms, iloadp, lloadp, loadp, inod_crk, iel_crk, iadc_crk, adsky_crk, cne_crk, procne_crk, iadsdp_crk, iadrcp_crk, ibufssg_io, ibc_ply, dmint2, ibordnode, elbuf_tab, por, nodedge, iad_edge, fr_edge, fr_nbedge, crknodiad, lgauge, gauge, igaup, ngaup, nodlevxf, dd_r2r_elem, nodglobxfe, sph2sol, sol2sph, irst, dmsph, wagap, xfem_tab, elcutc, nodenr, kxfenod2elc, enrtag, rthbu f, kxig3d, ixig3d, knot, wige, wsmcomp, stack, cputime_mp_glob, cputime_mp, tab_ump, poin_ump, sol2sph_typ, irunn_bis, addcsrect, iad_frnor, fr_nor, procnor, iad_fredg, fr_edg, drape_sh4n, drape_sh3n, tab_mat, nativ0_sms, multi_fvm, segquadfr, ms_2d, h3d_data, subsets, igrnod, igrbric, igrquad, igrsh4n, igrsh3n, igrtruss, igrbeam, igrspring, igrpart, igrsurf, forneqs, nloc_dmg, iskwp_l, knotlocpc, knotlocel, pinch_data, tag_skins6, ale_connectivity, xcell, xface, ne_nercvois, ne_nesdvois, ne_lercvois, ne_lesdvois, ibcscyc, lbcscyc, t_monvol, id_global_vois, face_vois, dynain_data, fcont_max, ebcs_tab, diffusion, kloadpinter, loadpinter, dgaploadint, drapeg, user_windows, output, interfaces, dt, loads, python, dpl0cld, vel0cld, ndamp_vrel, id_damp_vrel, fr_damp_vrel, ndamp_vrel_rbyg, names_and_titles, unitab, liflow, lrflow, glob_therm, pblast, rbe3)
Definition resol.F:633
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
subroutine spmd_exch_a_rb6_vrel(nrbdim, iad_rby, fr_rby6, icsize, rby6, rby6_c, size_rby6_c)
subroutine spmd_exch_a_rb6(nrbdim, iad_rby, fr_rby6, icsize, rbf6)
subroutine spmd_exch_a_rb6g(npby, rbf6)
subroutine startime(event, itask)
Definition timer.F:93
subroutine stoptime(event, itask)
Definition timer.F:135